Compile (Fortran 77) and run phy4231orbits.f.

Plot with phy4241orbits.plt.