Quadratic term in the mean anomaly.

miroslavbroz · Nov 12, 2024 · c28e09e · c28e09e
1 parent c648a28
commit c28e09e
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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)