Optimized solution of Kepler's equation (1972) [pdf] (ntrs.nasa.gov)

This document presents a detailed description of KEPLER, an IBM 360 computer program used for the solution of Kepler’s equation for eccentric anomaly. KEPLER employs a second-order Newton-Raphson differential correction process, and it is faster than previously developed programs by an order of magnitude.