Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

13-May-2025 16:29:20

pendulum_nonlinear_ode_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Solve pendulum_nonlinear_ode().

  parameters:
    g     = 9.81
    l     = 2
    m     = 1
    t0    = 0
    y0    = (3.11018,0)
    tstop = 50

  Initial energy             = 9.80758
  Energy needed to flip      = 9.81

pendulum_nonlinear_ode23():
  Test pendulum_nonlinear_ode() using ode23().
  Number of adaptive steps = 413
  Graphics saved in "pendulum_nonlinear_ode23_theta.png"
  Graphics saved in "pendulum_nonlinear_ode23_thetadot.png"
  Graphics saved in "pendulum_nonlinear_ode23_phase.png"
  Graphics saved in "pendulum_nonlinear_ode23_energy.png"

pendulum_nonlinear_midpoint():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Solve pendulum_nonlinear_ode() using midpoint().

  Number of equal steps to use = 1000
  Graphics saved in "pendulum_nonlinear_midpoint_theta.png"
  Graphics saved in "pendulum_nonlinear_midpoint_thetadot.png"
  Graphics saved in "pendulum_nonlinear_midpoint_phase.png"
  Graphics saved in "pendulum_nonlinear_midpoint_energy.png"

pendulum_nonlinear_period_test():
  Compute the period as a function of initial energy.

  parameters:
    g     = 9.81
    l     = 2
    m     = 1
    t0    = 0
    y0    = (3.11018,0)
    tstop = 50
    Energy needed to flip      = 9.81
  Graphics saved in "pendulum_nonlinear_period.png"

pendulum_nonlinear_ode_test():
  Normal end of execution.

13-May-2025 16:29:31