function u=dFu(t,y)
%
global k c;
%k=1;c=1/2;
u(1,1)=0;u(1,2)=1;
u(2,1)=-k*cos(y(1));u(2,2)=-c;
%u=u';

