function v=lin_rw(t,x)
v=[x(2),-x(1)+0.01*x(1)^2];