function dy=LIV1_eq1(t,y,g,z)  

%---2006.5.18-----
%%%   System
%==========================================================
%
%                       g
%                dd_x= --- x
%                       z
%
%           d_y1=y2
%           d_y2= g/z * y1
%===========================================================
%  y(1,1)•ÏˆÊ
%  y(2,1)‘¬“x
%  dy(1,1)‘¬“x
%  dy(2,1)‰Á‘¬“x

dy=zeros(2,1);

dy(1,1)=y(2,1); 

dy(2,1)= g/z * y(1,1) ;








