The earlier put up very briefly stated that the integral illustration for Bessel features was motived by fixing Kepler’s equation. This put up will go into extra element.
Kepler’s equation
There are a number of methods to explain the place of a planet in an elliptical orbit round a star. For historic causes, these descriptions have arcane names comparable to imply anomaly, true anomaly, and eccentric anomaly. This put up explains how these three are associated.
For this put up, it is sufficient to say that always you recognize imply anomaly M and need to know eccentric anomaly E. These are associated by way of Kepler’s equation
the place e is the eccentricity of the orbit. You’d like to resolve for E as a perform of M and e, however there’s no elementary approach to try this.
One option to resolve Kepler’s equation is to take a guess at E and plug it into the suitable hand aspect of
to get a brand new E, and hold iterating till the 2 sides are nearer collectively. I write extra about this right here.
One other strategy to fixing Kepler’s equation is to make use of Newton’s technique. I write extra about that right here.
Nonetheless one other strategy is to broaden E in a sine sequence and discover the sequence coefficients. A bonus to this strategy is that after you have the coefficients, you have got an expression for E as a perform of M, and you’ll plug in additional values of M with out having to resolve Kepler’s equation for every worth of M individually.
Sine sequence coefficients
Kepler’s equation is straightforward to resolve at E = 0 and at E = π. In each circumstances, E = M. So the perform E − M is zero at each ends of [0, π], which suggests we attempt to broaden E − M in a sine sequence
We then calculate the Fourier coefficients an as standard.
The second line makes use of integration by elements. The third line makes use of Kepler’s equation. The final line makes use of the definition of the Bessel features Jn given within the earlier put up.
