function f = fpend(y,t)
ome=1.;
f(1)=y(2);
f(2)=-ome*ome*sin(y(1));
end