%
%   fixed_fixed_beam_modes.m  ver 1.0
%
function[YY,YYd,YYdd,YYddd]=fixed_fixed_beam_modes(n,x,C,beta,sq_mass)
%
      YY=zeros(n,1);
     YYd=zeros(n,1);   
    YYdd=zeros(n,1);   
   YYddd=zeros(n,1);     
%
   for i=1:n
        arg=beta(i)*x;
          YY(i)=             (cosh(arg)-cos(arg))-C(i)*(sinh(arg)-sin(arg));
         YYd(i)=(beta(i))^1*((sinh(arg)+sin(arg))-C(i)*(cosh(arg)-cos(arg)));       
        YYdd(i)=(beta(i))^2*((cosh(arg)+cos(arg))-C(i)*(sinh(arg)+sin(arg)));  
       YYddd(i)=(beta(i))^3*((sinh(arg)-sin(arg))-C(i)*(cosh(arg)+cos(arg)));           
   end
%   
      YY=  YY/sq_mass;
    YYdd=YYdd/sq_mass;
   YYddd=YYddd/sq_mass;
