pendulum_midpoint:
  Test pendulum_ode using midpoint().
  Graphics saved in "pendulum_midpoint_theta.png"
  Graphics saved in "pendulum_midpoint_thetadot.png"
  Graphics saved in "pendulum_midpoint_energy.png"
  h(1) = 5.37345, h(n+1) = 5.37345