function value = predator_prime ( t, y )

%*****************************************************************************80
%
%% predator_prime returns the right hand side of the predator test problem.
%
  r = y(1);
  f = y(2);

  drdt =   2.0 * r - 0.001 * r * f;
  dfdt = -10.0 * f + 0.002 * r * f;

  value = [ drdt; dfdt ];

  return
end

