function yprime=f(t,y);
   yprime(1)=y(2);
   yprime(2)=-sin(y(1));
   yprime=yprime';
