Module containing the functions used to compute lpalace coefficients and their derivatives.
Computes the Laplace coefficient b_s^(j)(alpha) and its derivatives using the recursive relations given in Murray & Dermott (1999), including the case where alpha is close to 1 where the G(y) series solution is used.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=DP), | intent(in) | :: | alpha |
alpha is the ratio of the inner to outer semimajor axis (a_inner/a_outer) for the coefficient being computed |
||
| integer(kind=I4B), | intent(in) | :: | j |
j is the index of the Laplace coefficient being computed. |
||
| real(kind=DP), | intent(in) | :: | s |
s is the order of the Laplace coefficient being computed. |
||
| integer(kind=I4B), | intent(in) | :: | n |
n is the order of the derivative. n=0 gives the coefficient itself, n=1 gives the first derivative, etc. |
The value of the Laplace coefficient or its derivative, depending on n.