

function df=feuler_shape(f);

global FD
%f=[theta;dtheta;x;y];

df=zeros(4,1);

df(1)=f(2);
df(2)=-sin(f(1));
df(3)=FD*sin(f(1));
df(4)=FD*cos(f(1));


