function v=harm_osc(t,x)
v=[x(2),-sin(x(1))];