function dydt=damped(t,y)
dydt=[y(2); -3*sin(y(1))-0.2*y(2)];