function yprime=pendde(t,y,flag,s)
yprime=s*[y(2); -sin(y(1))];
