Method is modification of generalized Newton-Ralphson algorithm for analyzing two-point boundary problems. It constructs sequence of solutions that converge to precise dynamic solution in the sequence limit. Program calculates periodic orbits in either circular or elliptical restricted three-body problems
