

function df=dpendulum(f);

df=zeros(2,1);

df(1)=f(2);
df(2)=-sin(f(1));

