function dy=furiko3_eq(t,y,m,c,L)  

%---2011.11.18-----
%%%  nonlinear System
%  y(1,1)Šp“x
%  y(2,1)Šp‘¬“x
%  dy(1,1)Šp‘¬“x
%  dy(2,1)Šp‰Á‘¬“x

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

dy(2,1)=-c/(m*L^2) * y(2,1) -9.8/L * sin( y(1,1) );








