A simple script for demonstrating zero velocity curves in the restricted, circular three-body problem. Intended to allow students to visualize a contour plot of the Jacobi Constant for a given pair of masses. See the script itself for details.
Just type
python threebody.py
Edit the script to change different setup parameters.
Requires pylab. Written for python2.7.