%
  function [f_x, f_u ]   = rhsj(t,x,u)
%
% Evaluate the Jacobian for right-hand side 
%
  f_x = [ 0    .5;
          0    0 ];
%
  f_u = [ -sin(u); cos(u)]; 
%
% return
%
