Skip to content

Commit

Permalink
Correct implementation of Laguerre's method (#27)
Browse files Browse the repository at this point in the history
* Correct implementation of Laguerre's method

* Update time-since-periapsis-and-keplers-equation/universal-variables.md

---------

Co-authored-by: Bryan Weber <[email protected]>
  • Loading branch information
rootmos and bryanwweber committed Jun 16, 2024
1 parent 8ea9ffc commit 3e51bd7
Showing 1 changed file with 1 addition and 1 deletion.
Original file line number Diff line number Diff line change
Expand Up @@ -312,7 +312,7 @@ The Laguerre algorithm can be implemented as:

:::{math}
:label:
\chi_{i + 1} = \chi_{i} - \frac{n f(\chi_i)}{f'(\chi_i) \pm \left[\left(n - 1\right)^2 \left(f'(\chi_i\right)^2 - n\left(n - 1\right) f(\chi_i)f''(\chi_i)\right]}
\chi_{i + 1} = \chi_{i} - \frac{n f(\chi_i)}{f'(\chi_i) \pm \left[\left(n - 1\right)^2 \left(f'(\chi_i)\right)^2 - n\left(n - 1\right) f(\chi_i)f''(\chi_i)\right]^{1/2}}
:::

The sign ambiguity in the denominator is determined by taking the sign of the numerical value of $f'(\chi_i)$. In addition, the solution is relatively insensitive to the choice of the value of $n$, which is an integer constant. It seems as though $n = 5$ is a reasonable value. Choosing $n = 1$ gives the standard Newton's algorithm.
Expand Down

0 comments on commit 3e51bd7

Please sign in to comment.