%
%  ss_beam_modes.m  ver 1.0
%
function[YY,YYd,YYdd,YYddd]=ss_beam_modes(n,x,L,AC)
%
   YY=zeros(n,1);
 YYdd=zeros(n,1);
YYddd=zeros(n,1);
%   
for i=1:n
    arg=i*pi*x/L;
    c=(i*pi/L);
%
       AC_sin_arg=AC*sin(arg);
       AC_cos_arg=AC*cos(arg);
%
       YY(i)=     AC_sin_arg;
      YYd(i)=   c*AC_cos_arg;
     YYdd(i)=-c^2*AC_sin_arg;
    YYddd(i)=-c^3*AC_cos_arg;
end
