function dy=rigid(t,y)
dy = -0.5*y^2;



