predator_rk1  Solve predator-prey equation with the rk1 method
  Plot the solution, phase plane, conservation.

  Input value N not supplied.  Using N = 100000
  Graphics saved as "predator_rk1_100000.png"
  Graphics saved as "predator_rk1_phase_100000.png"