function u=Fu(t,y)
%
global k c;
%k=1;c=1/2;
u(1)=y(2);
u(2)=-c*y(2)-k*sin(y(1));
u=u';

