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