function dy = heiluri(t,y)
dy=[y(2);-sin(y(1))];