function [Cbeam] = CmatrixBeamGeneral ( k_w1, k_w2, k_theta1, k_theta2, beta);

%  characteristic matrix for beam vibration problem;
%  the two end supports are constrained by two transverse spring
%  as well as two rotational springs;


%input beta = beta*L;

%%%%%%%%%%%%%%%%%%%
%  data preparation
%%%%%%%%%%%%%%%%%%%

cc = cos (beta);
ss = sin (beta);

ch = cosh(beta);
sh = sinh(beta);

Cbeam = [     
 beta^3                         -k_w1              -beta^3                                 -k_w1; 
 -k_theta1                     -beta               -k_theta1                                beta;			 
-(beta^3*cc+k_w2*ss)      (beta^3*ss-k_w2*cc)     (beta^3*ch-k_w2*sh)       (beta^3*sh-k_w2*ch);
 (beta*ss - k_theta2*cc)  (beta*cc + k_theta2*ss)  -(beta*sh +k_theta2*ch)   -(beta*ch +k_theta2*sh)];
return;