function V=lin_rw_mat(t,A,x) J=[[0, 1];[ -1, 0]]; V=J*A;