function xout = rhs(x)
%
% 
xout    = x;   % make sure shapes agree
xout(1) = x(2);
xout(2) = x(3);
xout(3) = -x(1);