Computer program uses an iterative method to construct precisely periodic orbits which dynamically approximate solutions that converge to precise dynamical solutions in the limit of the sequence. The method used is a modification of the generalized Newton-Raphson algorithm used in analyzing two point boundary problems