From c28e09e0545a1357f9d9dcb798ab7ca6c000378d Mon Sep 17 00:00:00 2001 From: Miroslav Broz Date: Tue, 12 Nov 2024 08:38:19 +0100 Subject: [PATCH] Quadratic term in the mean anomaly. --- phoebe/dynamics/keplerian.py | 1 + 1 file changed, 1 insertion(+) diff --git a/phoebe/dynamics/keplerian.py b/phoebe/dynamics/keplerian.py index 2921e73c9..3d5e7fff7 100644 --- a/phoebe/dynamics/keplerian.py +++ b/phoebe/dynamics/keplerian.py @@ -141,6 +141,7 @@ def _keplerian(times, component=None): if dpdt != 0.0: P_ = P + dpdt*(t-t0) a *= (P_/P)**2 + M -= n*dpdt/(2.*P)*(t-t0)**2 if dperdt != 0.0: omega += dperdt*(t-t0)