The Lambert <Emphasis Type="Italic">W</Emphasis> function and solutions to Kepler’s equation

This paper presents a method for the truncation of infinite Fourier–Bessel representations for functions requiring a solution to Kepler’s equation. Use is made of the Lambert W function to solve for the desired index that bounds the remainder terms of the series, within the prescribed tolerance. The enforcement of a maximum on the number of Bessel functions is also useful in truncating the Bessel functions themselves, resulting in an analytical representation of the solution to a desired tolerance, without the use of infinite series.