function yprime = sinfunc(t,y)
yprime = -sin(y);
