%
  function x_dot   = rhs(t,x,FLAG,u)
%
% Evaluate the right-hand side of the
% differential equation for the state
%
  x_dot = ...
          [ cos(u(1)) + .5*x(2);
                sin(u(1))      ];
%
% return
%
