08-Jan-2022 08:59:35

predator_prey_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the predator prey ODE.

  Parameters:
    alpha = 2
    beta  = 0.001
    gamma = 10
    delta = 0.002
    t0    = 0
    y0    = (5000,100)
    tstop = 5

predator_prey_equilibrium_test()
  Verify that dy/dt=[0;0] and H=0 for any parameter values.

  (alpha,beta,gamma,delta)  dydt  H

  (0.814724,0.905792,0.126987,0.913376): (0,0)  0
  (0.0975404,0.278498,0.546882,0.957507): (0,0)  0
  (0.157613,0.970593,0.957167,0.485376): (0,0)  0
  (0.141886,0.421761,0.915736,0.792207): (0,0)  0
  (0.655741,0.0357117,0.849129,0.933993): (0,0)  0

predator_prey_euler():
  Use euler() to solve the predator prey ODE.

  Number of equal steps = 200
  Graphics saved as "predator_prey_euler_plot.png"
  Graphics saved as "predator_prey_euler_phase.png"
  Graphics saved as "predator_prey_euler_conservation.png"

predator_prey_midpoint()
  Use midpoint to solve the predator prey ODE.

  Number of equal steps = 200
  Graphics saved as "predator_prey_midpoint_plot.png"
  Graphics saved as "predator_prey_midpoint_phase.png"
  Graphics saved as "predator_prey_midpoint_conservation.png"

predator_prey_ode23():
  Use ode23() to solve the predator prey ODE.

  Number of variable steps = 126
  Graphics saved as "predator_prey_ode23_plot.png"
  Graphics saved as "predator_prey_ode23_phase.png"
  Graphics saved as "predator_prey_ode23_conservation.png"

predator_prey_ode45()
  Use ode45() to solve the predator prey ODE.

  Number of variable steps = 153
  Graphics saved as "predator_prey_ode45_plot.png"
  Graphics saved as "predator_prey_ode45_phase.png"
  Graphics saved as "predator_prey_ode45_conservation.png"

predator_prey_ode_test():
  Normal end of execution.

08-Jan-2022 08:59:43