@interstellar-tools/equations 0.60.2 → 0.61.1

package/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 Alexander Vassbotn Røyne-Helgesen
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

package/dist/categories/angle/compute-angle.d.ts

@@ -1,62 +1,51 @@
-import type { CelestialBodyType, Radians, TimeStepInterface } from '@interstellar-tools/types';
+import type { OrbitInterface, Radians, TimeStepInterface } from '@interstellar-tools/types';
 /**
- * Computes the orbital angle (true anomaly, $ν$) of a celestial body for a given date and time step.
+ * Advances a body along its orbit by one time step and returns its new
+ * **true anomaly** ($\nu$).
  *
- * **Mathematical Explanation:**
+ * The orbit's `angle` is the starting true anomaly. There is no date input:
+ * call this once per frame or step, and feed the result back in as the next
+ * `angle` to move the body over time.
  *
- * This function determines the **true anomaly** ($ν$) based on the time elapsed since the
- * **J2000 epoch** (January 1, 2000) and the provided time step. The calculation follows these steps:
+ * **Step 1: Mean anomaly after the time step ($M$)**
  *
- * **Step 1: Compute Time Elapsed Since J2000**
- *
- * The number of days since **J2000** is computed as:
- * $$ \Delta T = \frac{t - t_{J2000}}{\text{MILLISECONDS\_PER\_DAY}} $$
+ * Done by {@link advanceMeanAnomaly}:
+ * $$ M = M_0 + n \cdot \Delta t, \qquad n = \frac{2\pi}{|P|} $$
  * where:
- * - $t$ is the current date in milliseconds.
- * - $t_{J2000}$ is **J2000** (2000-01-01T00:00:00Z).
- * - **MILLISECONDS_PER_DAY** is the number of milliseconds in one day.
- *
- * **Step 2: Compute the Mean Anomaly ($M$)**
+ * - $M_0$ is the mean anomaly at the start, converted from the true anomaly
+ *   `angle` (equal to it when $e = 0$).
+ * - $n$ is the mean motion, from the period $P$ converted to days. The sign of
+ *   $P$ is ignored, so the angle always increases.
+ * - $\Delta t$ is `timeStep.value` in days, clamped to at least $10^{-5}$ days
+ *   and at most $10\,|P|$ ($|P|$ when $e \geq 0.9$).
  *
- * The **mean anomaly** is calculated as:
- * $$ M = M_0 + n \cdot (\Delta T + \text{timeStep}) $$
- * where:
- * - $M_0$ is the initial mean anomaly at **J2000**.
- * - $n$ is the mean motion (orbital angular velocity).
- * - $\text{timeStep}$ is the time step in days.
+ * The result is wrapped to $[0, 2\pi)$.
  *
- * **Step 3: Solve Kepler’s Equation for Eccentric Anomaly ($E$)**
+ * **Step 2: Solve Kepler's equation for the eccentric anomaly ($E$)**
  *
- * Kepler’s equation:
  * $$ M = E - e \sin(E) $$
- * is solved numerically to obtain the **eccentric anomaly** ($E$).
- *
- * **Step 4: Convert Eccentric Anomaly ($E$) to True Anomaly ($ν$)**
+ * is solved numerically by {@link solveKepler}.
  *
- * Using the relation:
- * $$ ν = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
+ * **Step 3: Convert the eccentric anomaly to the true anomaly ($\nu$)**
  *
- * **Step 5: Adjust the True Anomaly for Certain Bodies**
- * - For **comets**, the final **true anomaly** is adjusted by subtracting the **eccentricity** ($e$).
- * - For **stars and other celestial objects**, the computed **true anomaly** is used directly.
+ * Done by {@link eccentricToTrueAnomaly}:
+ * $$ \nu = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the true anomaly.
- * @param {TimeStepInterface} timeStep - The time step in days since the last frame.
- * @returns {Radians} The computed true anomaly in radians.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TimeStepInterface} timeStep - How far to advance, in days.
+ * @returns {Radians} The true anomaly after the time step, in radians.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
  *
  * @example
  * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
  * import { computeAngle } from '@interstellar-tools/equations';
  *
- * const earth: CelestialBodyType = {
- *   name: 'Earth',
- *   e: 0.0167, // Orbital eccentricity
- *   angle: 0,  // Initial angle
- * };
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ * const timeStep: TimeStepInterface = { value: 1, unit: 'd' };
  *
- * const timeStep: TimeStepInterface = { days: 1 }; // Time step of 1 day
- * const angle: Radians = computeAngle(earth, timeStep);
- * console.log(angle); // Output: Computed true anomaly in radians
+ * // Earth's true anomaly one day after its dataset position
+ * const angle: Radians = computeAngle(earth.orbit!, timeStep);
  * ```
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
@@ -65,5 +54,5 @@ import type { CelestialBodyType, Radians, TimeStepInterface } from '@interstella
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Angles
  */
-export declare const computeAngle: (body: CelestialBodyType, timeStep: TimeStepInterface) => Radians;
+export declare const computeAngle: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, timeStep: TimeStepInterface) => Radians;
 //# sourceMappingURL=compute-angle.d.ts.map

package/dist/categories/angle/compute-angle.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"compute-angle.d.ts","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EACV,iBAAiB,EACjB,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAMnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiEG;AACH,eAAO,MAAM,YAAY,GACvB,MAAM,iBAAiB,EACvB,UAAU,iBAAiB,KAC1B,OASF,CAAC"}
+{"version":3,"file":"compute-angle.d.ts","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EACV,cAAc,EACd,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAMnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsDG;AACH,eAAO,MAAM,YAAY,GACvB,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,UAAU,iBAAiB,KAC1B,OASF,CAAC"}

package/dist/categories/angle/compute-angle.js

@@ -1,64 +1,53 @@
+import { advanceMeanAnomaly } from '../anomalies/advance-mean-anomaly.js';
 import { eccentricToTrueAnomaly } from '../anomalies/eccentric-to-true-anomaly.js';
-import { meanToEccentricAnomaly } from '../anomalies/mean-to-eccentric-anomaly.js';
 import { solveKepler } from '../kepler/solve-kepler.js';
 /**
- * Computes the orbital angle (true anomaly, $ν$) of a celestial body for a given date and time step.
+ * Advances a body along its orbit by one time step and returns its new
+ * **true anomaly** ($\nu$).
  *
- * **Mathematical Explanation:**
+ * The orbit's `angle` is the starting true anomaly. There is no date input:
+ * call this once per frame or step, and feed the result back in as the next
+ * `angle` to move the body over time.
  *
- * This function determines the **true anomaly** ($ν$) based on the time elapsed since the
- * **J2000 epoch** (January 1, 2000) and the provided time step. The calculation follows these steps:
+ * **Step 1: Mean anomaly after the time step ($M$)**
  *
- * **Step 1: Compute Time Elapsed Since J2000**
- *
- * The number of days since **J2000** is computed as:
- * $$ \Delta T = \frac{t - t_{J2000}}{\text{MILLISECONDS\_PER\_DAY}} $$
+ * Done by {@link advanceMeanAnomaly}:
+ * $$ M = M_0 + n \cdot \Delta t, \qquad n = \frac{2\pi}{|P|} $$
  * where:
- * - $t$ is the current date in milliseconds.
- * - $t_{J2000}$ is **J2000** (2000-01-01T00:00:00Z).
- * - **MILLISECONDS_PER_DAY** is the number of milliseconds in one day.
- *
- * **Step 2: Compute the Mean Anomaly ($M$)**
+ * - $M_0$ is the mean anomaly at the start, converted from the true anomaly
+ *   `angle` (equal to it when $e = 0$).
+ * - $n$ is the mean motion, from the period $P$ converted to days. The sign of
+ *   $P$ is ignored, so the angle always increases.
+ * - $\Delta t$ is `timeStep.value` in days, clamped to at least $10^{-5}$ days
+ *   and at most $10\,|P|$ ($|P|$ when $e \geq 0.9$).
  *
- * The **mean anomaly** is calculated as:
- * $$ M = M_0 + n \cdot (\Delta T + \text{timeStep}) $$
- * where:
- * - $M_0$ is the initial mean anomaly at **J2000**.
- * - $n$ is the mean motion (orbital angular velocity).
- * - $\text{timeStep}$ is the time step in days.
+ * The result is wrapped to $[0, 2\pi)$.
  *
- * **Step 3: Solve Kepler’s Equation for Eccentric Anomaly ($E$)**
+ * **Step 2: Solve Kepler's equation for the eccentric anomaly ($E$)**
  *
- * Kepler’s equation:
  * $$ M = E - e \sin(E) $$
- * is solved numerically to obtain the **eccentric anomaly** ($E$).
- *
- * **Step 4: Convert Eccentric Anomaly ($E$) to True Anomaly ($ν$)**
+ * is solved numerically by {@link solveKepler}.
  *
- * Using the relation:
- * $$ ν = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
+ * **Step 3: Convert the eccentric anomaly to the true anomaly ($\nu$)**
  *
- * **Step 5: Adjust the True Anomaly for Certain Bodies**
- * - For **comets**, the final **true anomaly** is adjusted by subtracting the **eccentricity** ($e$).
- * - For **stars and other celestial objects**, the computed **true anomaly** is used directly.
+ * Done by {@link eccentricToTrueAnomaly}:
+ * $$ \nu = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the true anomaly.
- * @param {TimeStepInterface} timeStep - The time step in days since the last frame.
- * @returns {Radians} The computed true anomaly in radians.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TimeStepInterface} timeStep - How far to advance, in days.
+ * @returns {Radians} The true anomaly after the time step, in radians.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
  *
  * @example
  * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
  * import { computeAngle } from '@interstellar-tools/equations';
  *
- * const earth: CelestialBodyType = {
- *   name: 'Earth',
- *   e: 0.0167, // Orbital eccentricity
- *   angle: 0,  // Initial angle
- * };
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ * const timeStep: TimeStepInterface = { value: 1, unit: 'd' };
  *
- * const timeStep: TimeStepInterface = { days: 1 }; // Time step of 1 day
- * const angle: Radians = computeAngle(earth, timeStep);
- * console.log(angle); // Output: Computed true anomaly in radians
+ * // Earth's true anomaly one day after its dataset position
+ * const angle: Radians = computeAngle(earth.orbit!, timeStep);
  * ```
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
@@ -67,13 +56,13 @@ import { solveKepler } from '../kepler/solve-kepler.js';
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Angles
  */
-export const computeAngle = (body, timeStep) => {
+export const computeAngle = (orbit, timeStep) => {
     // Compute mean anomaly (M) with time step
-    const M = meanToEccentricAnomaly(body, timeStep);
+    const M = advanceMeanAnomaly(orbit, timeStep);
     // Solve Kepler's equation for eccentric anomaly (E)
-    const E = solveKepler(M, body.e);
+    const E = solveKepler(M, orbit.e);
     // Convert eccentric anomaly (E) to true anomaly (ν)
-    const V = eccentricToTrueAnomaly(E, body.e);
+    const V = eccentricToTrueAnomaly(E, orbit.e);
     return V;
 };
 //# sourceMappingURL=compute-angle.js.map

package/dist/categories/angle/compute-angle.js.map

@@ -1 +1 @@
-{"version":3,"file":"compute-angle.js","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,WAAW,EAAE,MAAM,2BAA2B,CAAC;AAExD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiEG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAC1B,IAAuB,EACvB,QAA2B,EAClB,EAAE;IACX,0CAA0C;IAC1C,MAAM,CAAC,GAAY,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;IAC1D,oDAAoD;IACpD,MAAM,CAAC,GAAY,WAAW,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,oDAAoD;IACpD,MAAM,CAAC,GAAY,sBAAsB,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAErD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}
+{"version":3,"file":"compute-angle.js","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,kBAAkB,EAAE,MAAM,sCAAsC,CAAC;AAC1E,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,WAAW,EAAE,MAAM,2BAA2B,CAAC;AAExD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsDG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAC1B,KAAqD,EACrD,QAA2B,EAClB,EAAE;IACX,0CAA0C;IAC1C,MAAM,CAAC,GAAY,kBAAkB,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;IACvD,oDAAoD;IACpD,MAAM,CAAC,GAAY,WAAW,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAC3C,oDAAoD;IACpD,MAAM,CAAC,GAAY,sBAAsB,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAEtD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}

package/dist/categories/anomalies/advance-mean-anomaly.d.ts

@@ -0,0 +1,50 @@
+import type { OrbitInterface, Radians, TemporalInterface } from '@interstellar-tools/types';
+/**
+ * Advances an orbit by a time step and returns its new **mean anomaly** ($M$).
+ *
+ * The orbit's `angle` (true anomaly) is the starting point. To get the eccentric
+ * anomaly from the result, pass it to {@link solveKepler}.
+ *
+ * **Mathematical Explanation:**
+ *
+ * The **mean anomaly** is a measure of the position of an orbiting body as if it moved
+ * with uniform angular motion along a circular orbit with the same period as the actual
+ * elliptical orbit. It is computed as:
+ *
+ * $$ M = M_0 + n \cdot \Delta T $$
+ *
+ * where:
+ * - $M_0$ is the **initial mean anomaly** (converted from true anomaly if necessary).
+ * - $n$ is the **mean motion**, given by:
+ *   $$ n = \frac{2\pi}{P} $$
+ *   where $P$ is the orbital period in days.
+ * - $\Delta T$ is the time step, converted to days from its unit.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TemporalInterface} timeStep - The time step over which to compute the change, in any temporal unit.
+ * @returns {Radians} The computed mean anomaly in radians.
+ *
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside the range $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { advanceMeanAnomaly } from '@interstellar-tools/equations';
+ *
+ * const marsOrbit = {
+ *   e: 0.0934, // Eccentricity of Mars' orbit
+ *   angle: 1.047 as Radians, // Initial true anomaly (in radians)
+ *   period: { value: 687, unit: 'd' } // Orbital period in days
+ * };
+ *
+ * const timeStep: TemporalInterface = { value: 1, unit: 'd' }; // 1-day step
+ * const meanAnomaly = advanceMeanAnomaly(marsOrbit, timeStep);
+ * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
+ * ```
+ *
+ * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.
+ * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
+ * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
+ * @group Anomalies
+ */
+export declare const advanceMeanAnomaly: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, timeStep: TemporalInterface) => Radians;
+//# sourceMappingURL=advance-mean-anomaly.d.ts.map

package/dist/categories/anomalies/advance-mean-anomaly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"advance-mean-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/advance-mean-anomaly.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EACV,cAAc,EACd,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAQnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AAEH,eAAO,MAAM,kBAAkB,GAC7B,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,UAAU,iBAAiB,KAC1B,OAqCF,CAAC"}

package/dist/categories/anomalies/{mean-to-eccentric-anomaly.js → advance-mean-anomaly.js}

@@ -5,7 +5,10 @@ import { trueToMeanAnomaly } from './true-to-mean-anomaly.js';
 import { areEqual } from './utils/are-equal.js';
 /* node:coverage disable */
 /**
- * Computes the **mean anomaly** ($M$) of a celestial body for a given time step.
+ * Advances an orbit by a time step and returns its new **mean anomaly** ($M$).
+ *
+ * The orbit's `angle` (true anomaly) is the starting point. To get the eccentric
+ * anomaly from the result, pass it to {@link solveKepler}.
  *
  * **Mathematical Explanation:**
  *
@@ -20,52 +23,53 @@ import { areEqual } from './utils/are-equal.js';
  * - $n$ is the **mean motion**, given by:
  *   $$ n = \frac{2\pi}{P} $$
  *   where $P$ is the orbital period in days.
- * - $\Delta T$ is the time step.
+ * - $\Delta T$ is the time step, converted to days from its unit.
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the mean anomaly.
- * @param {TemporalInterface} timeStep - The time step over which to compute the change.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TemporalInterface} timeStep - The time step over which to compute the change, in any temporal unit.
  * @returns {Radians} The computed mean anomaly in radians.
  *
- * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the body's eccentricity is outside the range $0 \leq e < 1$.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside the range $0 \leq e < 1$.
  *
  * @example
  * ```ts
- * import { meanToEccentricAnomaly } from '@interstellar-tools/equations';
+ * import { advanceMeanAnomaly } from '@interstellar-tools/equations';
  *
- * const mars: CelestialBodyType = {
- *   name: 'Mars',
+ * const marsOrbit = {
  *   e: 0.0934, // Eccentricity of Mars' orbit
- *   angle: 1.047, // Initial true anomaly (in radians)
- *   period: { value: 687, unit: 'day' }, // Orbital period in days
+ *   angle: 1.047 as Radians, // Initial true anomaly (in radians)
+ *   period: { value: 687, unit: 'd' } // Orbital period in days
  * };
  *
- * const timeStep: TemporalInterface = { value: 1, unit: 'day' }; // 1-day step
- * const meanAnomaly = meanToEccentricAnomaly(mars, timeStep);
+ * const timeStep: TemporalInterface = { value: 1, unit: 'd' }; // 1-day step
+ * const meanAnomaly = advanceMeanAnomaly(marsOrbit, timeStep);
  * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
  * ```
  *
+ * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.
  * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Anomalies
  */
 /* node:coverage enable */
-export const meanToEccentricAnomaly = (body, timeStep) => {
-    const { e, angle, period } = body;
+export const advanceMeanAnomaly = (orbit, timeStep) => {
+    const { e, angle, period } = orbit;
     if (e < 0 || e >= 1) {
         throw new RangeError(`Invalid eccentricity: ${e}. Eccentricity must be in the range [0, 1].`);
     }
-    // Convert period to days and compute mean motion
+    // Convert period and time step to days and compute mean motion
     const periodInDays = convertTemporalUnit(period, 'd').value;
+    const timeStepInDays = convertTemporalUnit(timeStep, 'd').value;
     const meanMotion = TWO_PI / Math.abs(periodInDays);
     // Convert true anomaly (V) to mean anomaly (M0) if necessary
     const M0 = e === 0 ? angle : trueToMeanAnomaly(angle, e);
-    if (timeStep.value === 0) {
+    if (timeStepInDays === 0) {
         return M0;
     }
     // Clamping small time steps to prevent unnecessary updates
     const minTimeStep = 1e-5; // Minimum allowed time step
     const maxTimeStepFactor = e >= 0.9 ? 1 : 10; // Stricter limit for highly eccentric orbits
-    const clampedDeltaT = Math.max(minTimeStep, Math.min(timeStep.value, maxTimeStepFactor * Math.abs(periodInDays)));
+    const clampedDeltaT = Math.max(minTimeStep, Math.min(timeStepInDays, maxTimeStepFactor * Math.abs(periodInDays)));
     // Compute new mean anomaly with hysteresis filter and sign preservation
     const prevM = M0;
     let newM = wrapAngle(M0 + meanMotion * clampedDeltaT);
@@ -74,4 +78,4 @@ export const meanToEccentricAnomaly = (body, timeStep) => {
     }
     return newM;
 };
-//# sourceMappingURL=mean-to-eccentric-anomaly.js.map
+//# sourceMappingURL=advance-mean-anomaly.js.map

package/dist/categories/anomalies/advance-mean-anomaly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"advance-mean-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/advance-mean-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AACvD,OAAO,EAAE,mBAAmB,EAAE,MAAM,8BAA8B,CAAC;AAOnE,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAC9D,OAAO,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAEhD,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAChC,KAAqD,EACrD,QAA2B,EAClB,EAAE;IACX,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC;IAEnC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,yBAAyB,CAAC,6CAA6C,CACxE,CAAC;IACJ,CAAC;IAED,+DAA+D;IAC/D,MAAM,YAAY,GAAG,mBAAmB,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAC5D,MAAM,cAAc,GAAG,mBAAmB,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAChE,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IACnD,6DAA6D;IAC7D,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAEzD,IAAI,cAAc,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,EAAa,CAAC;IACvB,CAAC;IAED,2DAA2D;IAC3D,MAAM,WAAW,GAAG,IAAI,CAAC,CAAC,4BAA4B;IACtD,MAAM,iBAAiB,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,6CAA6C;IAC1F,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,WAAW,EACX,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,iBAAiB,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CACrE,CAAC;IACF,wEAAwE;IACxE,MAAM,KAAK,GAAG,EAAE,CAAC;IAEjB,IAAI,IAAI,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC,CAAC;IAEtD,IAAI,QAAQ,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAC1B,IAAI,GAAG,KAAK,CAAC;IACf,CAAC;IAED,OAAO,IAAe,CAAC;AACzB,CAAC,CAAC"}

package/dist/categories/kepler/solve-kepler.d.ts

@@ -6,6 +6,9 @@ import type { Radians } from '@interstellar-tools/types';
  * - **High-eccentricity solver** for extreme orbits ($e > 0.9$).
  * - **Bisection fallback** if the selected solver doesn't return a root.
  *
+ * This is the conversion from **mean anomaly** to **eccentric anomaly**. Use
+ * {@link eccentricToTrueAnomaly} on the result to get the true anomaly.
+ *
  * ---
  *
  * **Mathematical Explanation:**

package/dist/categories/kepler/solve-kepler.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"solve-kepler.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAQzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoFG;AACH,eAAO,MAAM,WAAW,GACtB,GAAG,OAAO,EACV,GAAG,MAAM,EACT,gBAAY,EACZ,kBAAgB,KACf,OAuBF,CAAC"}
+{"version":3,"file":"solve-kepler.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAQzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuFG;AACH,eAAO,MAAM,WAAW,GACtB,GAAG,OAAO,EACV,GAAG,MAAM,EACT,gBAAY,EACZ,kBAAgB,KACf,OAuBF,CAAC"}

package/dist/categories/kepler/solve-kepler.js

@@ -10,6 +10,9 @@ import { solveKeplerNewtonRaphson } from './solve-kepler-newton-raphson.js';
  * - **High-eccentricity solver** for extreme orbits ($e > 0.9$).
  * - **Bisection fallback** if the selected solver doesn't return a root.
  *
+ * This is the conversion from **mean anomaly** to **eccentric anomaly**. Use
+ * {@link eccentricToTrueAnomaly} on the result to get the true anomaly.
+ *
  * ---
  *
  * **Mathematical Explanation:**

package/dist/categories/kepler/solve-kepler.js.map

@@ -1 +1 @@
-{"version":3,"file":"solve-kepler.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAGvD,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE7C,OAAO,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AACnE,OAAO,EAAE,2BAA2B,EAAE,MAAM,qCAAqC,CAAC;AAClF,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAE5E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoFG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,CAAU,EACV,CAAS,EACT,OAAO,GAAG,EAAE,EACZ,SAAS,GAAG,IAAI,EACP,EAAE;IACX,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IAC9E,CAAC;IAED,8DAA8D;IAC9D,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;IAEjC,wDAAwD;IACxD,IAAI,CAAC,GACH,CAAC,GAAG,GAAG;QACL,CAAC,CAAC,2BAA2B,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC;QACxD,CAAC,CAAC,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAE1D,6EAA6E;IAC7E,sEAAsE;IACtE,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAEnD,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,MAAM,GAAG,QAAQ,CAAC,IAAI,SAAS,CAAC,EAAE,CAAC;QAC1D,CAAC,GAAG,oBAAoB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IACtD,CAAC;IAED,OAAO,OAAO,CAAC,CAAC,CAAY,CAAC;AAC/B,CAAC,CAAC"}
+{"version":3,"file":"solve-kepler.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAGvD,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE7C,OAAO,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AACnE,OAAO,EAAE,2BAA2B,EAAE,MAAM,qCAAqC,CAAC;AAClF,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAE5E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuFG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,CAAU,EACV,CAAS,EACT,OAAO,GAAG,EAAE,EACZ,SAAS,GAAG,IAAI,EACP,EAAE;IACX,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IAC9E,CAAC;IAED,8DAA8D;IAC9D,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;IAEjC,wDAAwD;IACxD,IAAI,CAAC,GACH,CAAC,GAAG,GAAG;QACL,CAAC,CAAC,2BAA2B,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC;QACxD,CAAC,CAAC,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAE1D,6EAA6E;IAC7E,sEAAsE;IACtE,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAEnD,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,MAAM,GAAG,QAAQ,CAAC,IAAI,SAAS,CAAC,EAAE,CAAC;QAC1D,CAAC,GAAG,oBAAoB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IACtD,CAAC;IAED,OAAO,OAAO,CAAC,CAAC,CAAY,CAAC;AAC/B,CAAC,CAAC"}

package/dist/categories/orbits/circular-speed.d.ts

@@ -21,6 +21,20 @@
  *
  * :::
  *
+ * ::: practice The International Space Station
+ *
+ * **7.66 km/s**
+ *
+ * _Circular speed at 420 km altitude_
+ *
+ * That's about 27 600 km/h, fast enough to circle Earth every 93 minutes.
+ *
+ * ```ts
+ * circularSpeed(6.371e6 + 420e3, 3.986004418e14); // ≈ 7661 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Orbital radius ($ r $) (m), must be finite and **> 0**.
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
  * @returns {number} Circular speed ($ v_c $) (m/s).

package/dist/categories/orbits/circular-speed.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"circular-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,eAAO,MAAM,aAAa,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAUrD,CAAC"}
+{"version":3,"file":"circular-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAqDG;AACH,eAAO,MAAM,aAAa,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAUrD,CAAC"}

package/dist/categories/orbits/circular-speed.js

@@ -21,6 +21,20 @@
  *
  * :::
  *
+ * ::: practice The International Space Station
+ *
+ * **7.66 km/s**
+ *
+ * _Circular speed at 420 km altitude_
+ *
+ * That's about 27 600 km/h, fast enough to circle Earth every 93 minutes.
+ *
+ * ```ts
+ * circularSpeed(6.371e6 + 420e3, 3.986004418e14); // ≈ 7661 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Orbital radius ($ r $) (m), must be finite and **> 0**.
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
  * @returns {number} Circular speed ($ v_c $) (m/s).

package/dist/categories/orbits/circular-speed.js.map

@@ -1 +1 @@
-{"version":3,"file":"circular-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC7D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;AAC3B,CAAC,CAAC"}
+{"version":3,"file":"circular-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAqDG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC7D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;AAC3B,CAAC,CAAC"}

package/dist/categories/orbits/escape-speed.d.ts

@@ -22,6 +22,22 @@
  *
  * :::
  *
+ * ::: practice Leaving Earth
+ *
+ * **11.19 km/s**
+ *
+ * _Escape speed at Earth's surface_
+ *
+ * That's √2 times the 7.91 km/s it takes to orbit just above the surface, which
+ * is why a rocket for the Moon or Mars is so much bigger than one for low Earth
+ * orbit.
+ *
+ * ```ts
+ * escapeSpeed(6.371e6, 3.986004418e14); // ≈ 11 186 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Radius ($ r $) (m), distance from the central body’s center (must be finite and **> 0**).
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
  * @returns {number} Escape speed ($ v_{\text{esc}} $) (m/s).

package/dist/categories/orbits/escape-speed.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"escape-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/escape-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+CG;AACH,eAAO,MAAM,WAAW,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAUnD,CAAC"}
+{"version":3,"file":"escape-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/escape-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+DG;AACH,eAAO,MAAM,WAAW,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAUnD,CAAC"}

package/dist/categories/orbits/escape-speed.js

@@ -22,6 +22,22 @@
  *
  * :::
  *
+ * ::: practice Leaving Earth
+ *
+ * **11.19 km/s**
+ *
+ * _Escape speed at Earth's surface_
+ *
+ * That's √2 times the 7.91 km/s it takes to orbit just above the surface, which
+ * is why a rocket for the Moon or Mars is so much bigger than one for low Earth
+ * orbit.
+ *
+ * ```ts
+ * escapeSpeed(6.371e6, 3.986004418e14); // ≈ 11 186 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Radius ($ r $) (m), distance from the central body’s center (must be finite and **> 0**).
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
  * @returns {number} Escape speed ($ v_{\text{esc}} $) (m/s).

package/dist/categories/orbits/escape-speed.js.map

@@ -1 +1 @@
-{"version":3,"file":"escape-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/escape-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+CG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC3D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACjC,CAAC,CAAC"}
+{"version":3,"file":"escape-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/escape-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+DG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC3D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACjC,CAAC,CAAC"}

package/dist/categories/orbits/kepler-period.d.ts

@@ -18,6 +18,21 @@
  *
  * :::
  *
+ * ::: practice Geostationary orbit
+ *
+ * **23.93 h**
+ *
+ * _Orbital period 42 164 km from Earth's centre_
+ *
+ * That's one sidereal day, the time Earth takes to turn once relative to the
+ * stars, so a satellite there stays over the same spot on the equator.
+ *
+ * ```ts
+ * keplerPeriod(42164e3, 3.986004418e14); // ≈ 86 164 s
+ * ```
+ *
+ * :::
+ *
  * @param {number} a  Semi-major axis (m), must be finite and **> 0**.
  * @param {number} mu Gravitational parameter (m³/s²), must be finite and **> 0**.
  * @returns {number} Orbital period ($ T $) (s).

package/dist/categories/orbits/kepler-period.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"kepler-period.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/kepler-period.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0CG;AACH,eAAO,MAAM,YAAY,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAYpD,CAAC"}
+{"version":3,"file":"kepler-period.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/kepler-period.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAyDG;AACH,eAAO,MAAM,YAAY,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAYpD,CAAC"}

package/dist/categories/orbits/kepler-period.js

@@ -18,6 +18,21 @@
  *
  * :::
  *
+ * ::: practice Geostationary orbit
+ *
+ * **23.93 h**
+ *
+ * _Orbital period 42 164 km from Earth's centre_
+ *
+ * That's one sidereal day, the time Earth takes to turn once relative to the
+ * stars, so a satellite there stays over the same spot on the equator.
+ *
+ * ```ts
+ * keplerPeriod(42164e3, 3.986004418e14); // ≈ 86 164 s
+ * ```
+ *
+ * :::
+ *
  * @param {number} a  Semi-major axis (m), must be finite and **> 0**.
  * @param {number} mu Gravitational parameter (m³/s²), must be finite and **> 0**.
  * @returns {number} Orbital period ($ T $) (s).

package/dist/categories/orbits/kepler-period.js.map

@@ -1 +1 @@
-{"version":3,"file":"kepler-period.js","sourceRoot":"","sources":["../../../src/categories/orbits/kepler-period.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0CG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;QACpC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM;IAE5B,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;AAC1C,CAAC,CAAC"}
+{"version":3,"file":"kepler-period.js","sourceRoot":"","sources":["../../../src/categories/orbits/kepler-period.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAyDG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;QACpC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM;IAE5B,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;AAC1C,CAAC,CAAC"}

package/dist/categories/orbits/mean-motion.d.ts

@@ -18,6 +18,21 @@
  * Common usage: reference rate for relative motion models (e.g., Clohessy–Wiltshire / Hill),
  * orbital period estimates, and circular-orbit dynamics.
  *
+ * ::: practice The ISS, orbit by orbit
+ *
+ * **15.51 orbits/day**
+ *
+ * _Mean motion at 420 km altitude_
+ *
+ * The function returns 0.001128 rad/s. Multiply by 86 400 s and divide by 2π for
+ * orbits per day, the unit used in a satellite's two-line elements (TLEs).
+ *
+ * ```ts
+ * meanMotion(3.986004418e14, 6.371e6 + 420e3); // ≈ 0.001128 rad/s
+ * ```
+ *
+ * :::
+ *
  * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @param a - Semi-major axis $ a $ (finite, > 0).
  * @returns Mean motion $ n $ in **rad/s**.

package/dist/categories/orbits/mean-motion.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"mean-motion.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,eAAO,MAAM,UAAU,GAAI,IAAI,MAAM,EAAE,GAAG,MAAM,KAAG,MAYlD,CAAC"}
+{"version":3,"file":"mean-motion.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,eAAO,MAAM,UAAU,GAAI,IAAI,MAAM,EAAE,GAAG,MAAM,KAAG,MAYlD,CAAC"}

package/dist/categories/orbits/mean-motion.js

@@ -18,6 +18,21 @@
  * Common usage: reference rate for relative motion models (e.g., Clohessy–Wiltshire / Hill),
  * orbital period estimates, and circular-orbit dynamics.
  *
+ * ::: practice The ISS, orbit by orbit
+ *
+ * **15.51 orbits/day**
+ *
+ * _Mean motion at 420 km altitude_
+ *
+ * The function returns 0.001128 rad/s. Multiply by 86 400 s and divide by 2π for
+ * orbits per day, the unit used in a satellite's two-line elements (TLEs).
+ *
+ * ```ts
+ * meanMotion(3.986004418e14, 6.371e6 + 420e3); // ≈ 0.001128 rad/s
+ * ```
+ *
+ * :::
+ *
  * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @param a - Semi-major axis $ a $ (finite, > 0).
  * @returns Mean motion $ n $ in **rad/s**.

package/dist/categories/orbits/mean-motion.js.map

@@ -1 +1 @@
-{"version":3,"file":"mean-motion.js","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,EAAU,EAAE,CAAS,EAAU,EAAE;IAC1D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,gCAAgC,EAAE,EAAE,CAAC,CAAC;IAE5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,+BAA+B,CAAC,EAAE,CAAC,CAAC;IAE1D,IAAI,EAAE,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,EAAE,EAAE,CAAC,CAAC;IAErE,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC,CAAC"}
+{"version":3,"file":"mean-motion.js","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,EAAU,EAAE,CAAS,EAAU,EAAE;IAC1D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,gCAAgC,EAAE,EAAE,CAAC,CAAC;IAE5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,+BAA+B,CAAC,EAAE,CAAC,CAAC;IAE1D,IAAI,EAAE,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,EAAE,EAAE,CAAC,CAAC;IAErE,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC,CAAC"}

package/dist/categories/orbits/peri-apoapsis-radii.d.ts

@@ -17,6 +17,23 @@ import { PeriApoRadiiType } from '@interstellar-tools/types';
  * - Inputs: ($ a $) in **meters (m)**, ($ e $) dimensionless.
  * - Output: ($ r_p, r_a $) in **meters (m)**.
  *
+ * ::: practice Molniya orbit
+ *
+ * **545 km / 39 913 km**
+ *
+ * _Altitude at perigee and apogee_
+ *
+ * With e = 0.74, a Molniya satellite races through perigee and spends most of its
+ * 12 hour orbit slowly crossing apogee, high above the northern latitudes it
+ * serves.
+ *
+ * ```ts
+ * periApoapsisRadii(26600e3, 0.74);
+ * // rp = 6 916 km, ra = 46 284 km from Earth's centre
+ * ```
+ *
+ * :::
+ *
  * @param {number} a Semi-major axis ($ a $) (m). Must be finite and **≠ 0**.
  * @param {number} e Eccentricity ($ e $) (dimensionless). Must be finite and **≥ 0**.
  * @returns {{ rp: number; ra: number | null }} Object with periapsis `rp` and apoapsis `ra` (or `null` if undefined for hyperbolic).

package/dist/categories/orbits/peri-apoapsis-radii.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"peri-apoapsis-radii.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/peri-apoapsis-radii.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAG7D;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0CG;AAEH,eAAO,MAAM,iBAAiB,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,gBAsDxD,CAAC"}
+{"version":3,"file":"peri-apoapsis-radii.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/peri-apoapsis-radii.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAG7D;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2DG;AAEH,eAAO,MAAM,iBAAiB,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,gBAsDxD,CAAC"}

package/dist/categories/orbits/peri-apoapsis-radii.js

@@ -17,6 +17,23 @@
  * - Inputs: ($ a $) in **meters (m)**, ($ e $) dimensionless.
  * - Output: ($ r_p, r_a $) in **meters (m)**.
  *
+ * ::: practice Molniya orbit
+ *
+ * **545 km / 39 913 km**
+ *
+ * _Altitude at perigee and apogee_
+ *
+ * With e = 0.74, a Molniya satellite races through perigee and spends most of its
+ * 12 hour orbit slowly crossing apogee, high above the northern latitudes it
+ * serves.
+ *
+ * ```ts
+ * periApoapsisRadii(26600e3, 0.74);
+ * // rp = 6 916 km, ra = 46 284 km from Earth's centre
+ * ```
+ *
+ * :::
+ *
  * @param {number} a Semi-major axis ($ a $) (m). Must be finite and **≠ 0**.
  * @param {number} e Eccentricity ($ e $) (dimensionless). Must be finite and **≥ 0**.
  * @returns {{ rp: number; ra: number | null }} Object with periapsis `rp` and apoapsis `ra` (or `null` if undefined for hyperbolic).

package/dist/categories/orbits/peri-apoapsis-radii.js.map

@@ -1 +1 @@
-{"version":3,"file":"peri-apoapsis-radii.js","sourceRoot":"","sources":["../../../src/categories/orbits/peri-apoapsis-radii.ts"],"names":[],"mappings":"AAEA,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0CG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,CAAS,EAAE,CAAS,EAAoB,EAAE;IAC1E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,KAAK,CAAC,2BAA2B,CAAC,CAAC;IAC/C,CAAC;IAED,sBAAsB;IACtB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,MAAM,IAAI,KAAK,CACb,0GAA0G,CAC3G,CAAC;YACJ,CAAC;YAED,MAAM,IAAI,KAAK,CACb,mEAAmE,CACpE,CAAC;QACJ,CAAC;QAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvB,4CAA4C;QAC5C,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;YACxB,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,CAAC;QAED,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;IACpB,CAAC;IAED,aAAa;IACb,6CAA6C;IAC7C,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,iBAAiB;IAEzC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;QACd,4BAA4B;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAE1C,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;YACd,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC;QAC7B,CAAC;QAED,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IACtE,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC;AAC1B,CAAC,CAAC"}
+{"version":3,"file":"peri-apoapsis-radii.js","sourceRoot":"","sources":["../../../src/categories/orbits/peri-apoapsis-radii.ts"],"names":[],"mappings":"AAEA,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2DG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,CAAS,EAAE,CAAS,EAAoB,EAAE;IAC1E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,KAAK,CAAC,2BAA2B,CAAC,CAAC;IAC/C,CAAC;IAED,sBAAsB;IACtB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,MAAM,IAAI,KAAK,CACb,0GAA0G,CAC3G,CAAC;YACJ,CAAC;YAED,MAAM,IAAI,KAAK,CACb,mEAAmE,CACpE,CAAC;QACJ,CAAC;QAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvB,4CAA4C;QAC5C,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;YACxB,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,CAAC;QAED,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;IACpB,CAAC;IAED,aAAa;IACb,6CAA6C;IAC7C,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,iBAAiB;IAEzC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;QACd,4BAA4B;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAE1C,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;YACd,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC;QAC7B,CAAC;QAED,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IACtE,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC;AAC1B,CAAC,CAAC"}

package/dist/categories/orbits/specific-mechanical-energy.d.ts

@@ -25,6 +25,26 @@
  *
  * :::
  *
+ * ::: practice Energy of the ISS orbit
+ *
+ * **-29.35 MJ/kg**
+ *
+ * _Specific orbital energy at 420 km altitude_
+ *
+ * Negative means the orbit is bound. Escaping Earth from there takes 29.35 MJ
+ * more for every kilogram, which brings the energy up to zero.
+ *
+ * ```ts
+ * const mu = 3.986004418e14;
+ * const r = 6.371e6 + 420e3;
+ *
+ * const v = circularSpeed(r, mu);
+ *
+ * specificMechanicalEnergy(v, r, mu); // ≈ -2.935e7 J/kg
+ * ```
+ *
+ * :::
+ *
  * @param {number} v  Speed magnitude ($ v $) (m/s).
  * @param {number} r  Radius ($ r=\|\mathbf r\| $) from the central body (m).
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²).

package/dist/categories/orbits/specific-mechanical-energy.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"specific-mechanical-energy.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/specific-mechanical-energy.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,eAAO,MAAM,wBAAwB,GACnC,GAAG,MAAM,EACT,GAAG,MAAM,EACT,IAAI,MAAM,KACT,MAWF,CAAC"}
+{"version":3,"file":"specific-mechanical-energy.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/specific-mechanical-energy.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkEG;AACH,eAAO,MAAM,wBAAwB,GACnC,GAAG,MAAM,EACT,GAAG,MAAM,EACT,IAAI,MAAM,KACT,MAWF,CAAC"}

package/dist/categories/orbits/specific-mechanical-energy.js

@@ -25,6 +25,26 @@
  *
  * :::
  *
+ * ::: practice Energy of the ISS orbit
+ *
+ * **-29.35 MJ/kg**
+ *
+ * _Specific orbital energy at 420 km altitude_
+ *
+ * Negative means the orbit is bound. Escaping Earth from there takes 29.35 MJ
+ * more for every kilogram, which brings the energy up to zero.
+ *
+ * ```ts
+ * const mu = 3.986004418e14;
+ * const r = 6.371e6 + 420e3;
+ *
+ * const v = circularSpeed(r, mu);
+ *
+ * specificMechanicalEnergy(v, r, mu); // ≈ -2.935e7 J/kg
+ * ```
+ *
+ * :::
+ *
  * @param {number} v  Speed magnitude ($ v $) (m/s).
  * @param {number} r  Radius ($ r=\|\mathbf r\| $) from the central body (m).
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²).

package/dist/categories/orbits/specific-mechanical-energy.js.map

@@ -1 +1 @@
-{"version":3,"file":"specific-mechanical-energy.js","sourceRoot":"","sources":["../../../src/categories/orbits/specific-mechanical-energy.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,MAAM,CAAC,MAAM,wBAAwB,GAAG,CACtC,CAAS,EACT,CAAS,EACT,EAAU,EACF,EAAE;IACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;QAC9B,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;IAEpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;QAC/B,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAE9D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAChC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IAEvE,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AAC9B,CAAC,CAAC"}
+{"version":3,"file":"specific-mechanical-energy.js","sourceRoot":"","sources":["../../../src/categories/orbits/specific-mechanical-energy.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkEG;AACH,MAAM,CAAC,MAAM,wBAAwB,GAAG,CACtC,CAAS,EACT,CAAS,EACT,EAAU,EACF,EAAE;IACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;QAC9B,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;IAEpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;QAC/B,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAE9D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAChC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IAEvE,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AAC9B,CAAC,CAAC"}

package/dist/categories/orbits/sphere-of-influence-radius.d.ts

@@ -14,6 +14,23 @@
  * Units:
  * - Output `rSOI` has the same distance unit as `a` (since the mass ratio is dimensionless).
  *
+ * ::: practice Earth's sphere of influence
+ *
+ * **925 000 km**
+ *
+ * _Where Earth's gravity outweighs the Sun's_
+ *
+ * The Moon, at 384 400 km, is well inside it. Beyond this distance, a
+ * spacecraft's path is best described as an orbit around the Sun.
+ *
+ * ```ts
+ * const au = 1.495978707e11; // m
+ *
+ * sphereOfInfluenceRadius(au, 5.972e24, 1.989e30); // ≈ 9.245e8 m
+ * ```
+ *
+ * :::
+ *
  * @param {number} a Semi-major axis `a` of the secondary around the primary (finite, > 0).
  * @param {number} m Mass of the secondary body `m` (finite, > 0).
  * @param {number} M Mass of the primary body `M` (finite, > 0).

package/dist/categories/orbits/sphere-of-influence-radius.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"sphere-of-influence-radius.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/sphere-of-influence-radius.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiCG;AACH,eAAO,MAAM,uBAAuB,GAClC,GAAG,MAAM,EACT,GAAG,MAAM,EACT,GAAG,MAAM,KACR,MA0BF,CAAC"}
+{"version":3,"file":"sphere-of-influence-radius.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/sphere-of-influence-radius.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkDG;AACH,eAAO,MAAM,uBAAuB,GAClC,GAAG,MAAM,EACT,GAAG,MAAM,EACT,GAAG,MAAM,KACR,MA0BF,CAAC"}

package/dist/categories/orbits/sphere-of-influence-radius.js

@@ -14,6 +14,23 @@
  * Units:
  * - Output `rSOI` has the same distance unit as `a` (since the mass ratio is dimensionless).
  *
+ * ::: practice Earth's sphere of influence
+ *
+ * **925 000 km**
+ *
+ * _Where Earth's gravity outweighs the Sun's_
+ *
+ * The Moon, at 384 400 km, is well inside it. Beyond this distance, a
+ * spacecraft's path is best described as an orbit around the Sun.
+ *
+ * ```ts
+ * const au = 1.495978707e11; // m
+ *
+ * sphereOfInfluenceRadius(au, 5.972e24, 1.989e30); // ≈ 9.245e8 m
+ * ```
+ *
+ * :::
+ *
  * @param {number} a Semi-major axis `a` of the secondary around the primary (finite, > 0).
  * @param {number} m Mass of the secondary body `m` (finite, > 0).
  * @param {number} M Mass of the primary body `M` (finite, > 0).

package/dist/categories/orbits/sphere-of-influence-radius.js.map

@@ -1 +1 @@
-{"version":3,"file":"sphere-of-influence-radius.js","sourceRoot":"","sources":["../../../src/categories/orbits/sphere-of-influence-radius.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiCG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CACrC,CAAS,EACT,CAAS,EACT,CAAS,EACD,EAAE;IACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACpC,CAAC,CAAC"}
+{"version":3,"file":"sphere-of-influence-radius.js","sourceRoot":"","sources":["../../../src/categories/orbits/sphere-of-influence-radius.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkDG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CACrC,CAAS,EACT,CAAS,EACT,CAAS,EACD,EAAE;IACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACpC,CAAC,CAAC"}

package/dist/categories/orbits/vis-viva-speed.d.ts

@@ -22,6 +22,25 @@
  *
  * :::
  *
+ * ::: practice Departing for Mars
+ *
+ * **32.73 km/s**
+ *
+ * _Speed at the start of a Hohmann transfer to Mars_
+ *
+ * Earth moves around the Sun at 29.78 km/s, so a spacecraft leaving on the
+ * transfer orbit needs 2.95 km/s on top of Earth's own speed.
+ *
+ * ```ts
+ * const au = 1.495978707e11; // m
+ * const muSun = 1.32712440018e20; // m³/s²
+ * const a = ((1 + 1.524) / 2) * au; // transfer orbit
+ *
+ * visVivaSpeed(au, a, muSun); // ≈ 32 731 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Radius ($ r $) (m), distance to the attracting body.
  * @param {number} a  Semi-major axis ($ a $) (m). Use `Infinity` for parabolic, negative for hyperbolic.
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²).

package/dist/categories/orbits/vis-viva-speed.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"vis-viva-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/vis-viva-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,eAAO,MAAM,YAAY,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAsC/D,CAAC"}
+{"version":3,"file":"vis-viva-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/vis-viva-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0DG;AACH,eAAO,MAAM,YAAY,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAsC/D,CAAC"}

package/dist/categories/orbits/vis-viva-speed.js

@@ -22,6 +22,25 @@
  *
  * :::
  *
+ * ::: practice Departing for Mars
+ *
+ * **32.73 km/s**
+ *
+ * _Speed at the start of a Hohmann transfer to Mars_
+ *
+ * Earth moves around the Sun at 29.78 km/s, so a spacecraft leaving on the
+ * transfer orbit needs 2.95 km/s on top of Earth's own speed.
+ *
+ * ```ts
+ * const au = 1.495978707e11; // m
+ * const muSun = 1.32712440018e20; // m³/s²
+ * const a = ((1 + 1.524) / 2) * au; // transfer orbit
+ *
+ * visVivaSpeed(au, a, muSun); // ≈ 32 731 m/s
+ * ```
+ *
+ * :::
+ *
  * @param {number} r  Radius ($ r $) (m), distance to the attracting body.
  * @param {number} a  Semi-major axis ($ a $) (m). Use `Infinity` for parabolic, negative for hyperbolic.
  * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²).

package/dist/categories/orbits/vis-viva-speed.js.map

@@ -1 +1 @@
-{"version":3,"file":"vis-viva-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/vis-viva-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAU,EAAE;IACvE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,iEAAiE;QACjE,IAAI,CAAC,KAAK,QAAQ,EAAE,CAAC;YACnB,kBAAkB;YAClB,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACjC,CAAC;QAED,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IAED,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;IACpE,CAAC;IAED,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,gDAAgD;IACpE,MAAM,QAAQ,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC;IAErC,oEAAoE;IACpE,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvC,IAAI,QAAQ,GAAG,CAAC,GAAG;YAAE,OAAO,CAAC,CAAC,CAAC,yBAAyB;QAExD,MAAM,IAAI,KAAK,CACb,uCAAuC,EAAE,OAAO,CAAC,OAAO,CAAC,IAAI,CAC9D,CAAC;IACJ,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;AAC7B,CAAC,CAAC"}
+{"version":3,"file":"vis-viva-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/vis-viva-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0DG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAU,EAAE;IACvE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,iEAAiE;QACjE,IAAI,CAAC,KAAK,QAAQ,EAAE,CAAC;YACnB,kBAAkB;YAClB,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACjC,CAAC;QAED,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IAED,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;IACpE,CAAC;IAED,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,gDAAgD;IACpE,MAAM,QAAQ,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC;IAErC,oEAAoE;IACpE,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvC,IAAI,QAAQ,GAAG,CAAC,GAAG;YAAE,OAAO,CAAC,CAAC,CAAC,yBAAyB;QAExD,MAAM,IAAI,KAAK,CACb,uCAAuC,EAAE,OAAO,CAAC,OAAO,CAAC,IAAI,CAC9D,CAAC;IACJ,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;AAC7B,CAAC,CAAC"}

package/dist/index.d.ts

@@ -19,8 +19,8 @@ export * from './categories/helpers/radians.js';
 export * from './categories/helpers/rot-1.js';
 export * from './categories/helpers/rot-3.js';
 export * from './categories/helpers/transpose-3.js';
+export * from './categories/anomalies/advance-mean-anomaly.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
-export * from './categories/anomalies/mean-to-eccentric-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 export * from './categories/kepler/solve-kepler-bisection.js';
 export * from './categories/kepler/solve-kepler-high-eccentricity.js';

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAIjD,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAI3E,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAIpD,cAAc,qDAAqD,CAAC;AAEpE,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAI/D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAIpD,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAIhE,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAItD,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAIjD,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAI3E,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAIpD,cAAc,gDAAgD,CAAC;AAE/D,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAI/D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAIpD,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAIhE,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAItD,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}

package/dist/index.js

@@ -23,8 +23,8 @@ export * from './categories/helpers/rot-1.js';
 export * from './categories/helpers/rot-3.js';
 export * from './categories/helpers/transpose-3.js';
 // categories/anomalies
+export * from './categories/anomalies/advance-mean-anomaly.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
-export * from './categories/anomalies/mean-to-eccentric-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 // categories/kepler (solvers)
 export * from './categories/kepler/solve-kepler-bisection.js';

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,mBAAmB;AACnB,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAEjD,yBAAyB;AAEzB,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAE3E,qBAAqB;AAErB,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAEpD,uBAAuB;AAEvB,cAAc,qDAAqD,CAAC;AAEpE,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAE/D,8BAA8B;AAE9B,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAEpD,qBAAqB;AAErB,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAEhE,oBAAoB;AAEpB,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAEtD,wBAAwB;AAExB,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,mBAAmB;AACnB,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAEjD,yBAAyB;AAEzB,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAE3E,qBAAqB;AAErB,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAEpD,uBAAuB;AAEvB,cAAc,gDAAgD,CAAC;AAE/D,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAE/D,8BAA8B;AAE9B,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAEpD,qBAAqB;AAErB,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAEhE,oBAAoB;AAEpB,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAEtD,wBAAwB;AAExB,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@interstellar-tools/equations",
-  "version": "0.60.2",
+  "version": "0.61.1",
   "description": "Curated, well-tested orbital mechanics and dynamics formulas (e.g., Newtonian gravity, vis-viva, Kepler's 3rd law, escape/circular speeds, specific energy, specific angular momentum, peri/apoapsis radii).",
   "keywords": [
     "interstellar",
@@ -110,12 +110,13 @@
     "test:ci": "glob-bin -c \"node --import tsx --test --no-warnings --experimental-test-coverage --test-coverage-exclude='**/__tests__/helpers/**' --test-reporter=cobertura --test-reporter-destination=coverage/cobertura-coverage.xml --test-reporter=spec --test-reporter-destination=stdout\" \"**/__tests__/**/*.spec.ts\""
   },
   "dependencies": {
-    "@interstellar-tools/constants": "^0.60.0",
+    "@interstellar-tools/constants": "^0.61.0",
     "@interstellar-tools/temporal": "^0.60.0",
-    "@interstellar-tools/types": "^0.60.0"
+    "@interstellar-tools/types": "^0.61.0"
   },
   "devDependencies": {
     "@types/node": "^24.19.2",
+    "fast-check": "4.10.2",
     "prettier": "^3.9.10",
     "typescript": "^6.0.3"
   },

package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts

@@ -1,47 +0,0 @@
-import type { CelestialBodyType, Radians, TemporalInterface } from '@interstellar-tools/types';
-/**
- * Computes the **mean anomaly** ($M$) of a celestial body for a given time step.
- *
- * **Mathematical Explanation:**
- *
- * The **mean anomaly** is a measure of the position of an orbiting body as if it moved
- * with uniform angular motion along a circular orbit with the same period as the actual
- * elliptical orbit. It is computed as:
- *
- * $$ M = M_0 + n \cdot \Delta T $$
- *
- * where:
- * - $M_0$ is the **initial mean anomaly** (converted from true anomaly if necessary).
- * - $n$ is the **mean motion**, given by:
- *   $$ n = \frac{2\pi}{P} $$
- *   where $P$ is the orbital period in days.
- * - $\Delta T$ is the time step.
- *
- * @param {CelestialBodyType} body - The celestial body for which to compute the mean anomaly.
- * @param {TemporalInterface} timeStep - The time step over which to compute the change.
- * @returns {Radians} The computed mean anomaly in radians.
- *
- * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the body's eccentricity is outside the range $0 \leq e < 1$.
- *
- * @example
- * ```ts
- * import { meanToEccentricAnomaly } from '@interstellar-tools/equations';
- *
- * const mars: CelestialBodyType = {
- *   name: 'Mars',
- *   e: 0.0934, // Eccentricity of Mars' orbit
- *   angle: 1.047, // Initial true anomaly (in radians)
- *   period: { value: 687, unit: 'day' }, // Orbital period in days
- * };
- *
- * const timeStep: TemporalInterface = { value: 1, unit: 'day' }; // 1-day step
- * const meanAnomaly = meanToEccentricAnomaly(mars, timeStep);
- * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
- * ```
- *
- * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
- * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
- * @group Anomalies
- */
-export declare const meanToEccentricAnomaly: (body: CelestialBodyType, timeStep: TemporalInterface) => Radians;
-//# sourceMappingURL=mean-to-eccentric-anomaly.d.ts.map

package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"mean-to-eccentric-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/mean-to-eccentric-anomaly.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EACV,iBAAiB,EACjB,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAQnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2CG;AAEH,eAAO,MAAM,sBAAsB,GACjC,MAAM,iBAAiB,EACvB,UAAU,iBAAiB,KAC1B,OAoCF,CAAC"}

package/dist/categories/anomalies/mean-to-eccentric-anomaly.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"mean-to-eccentric-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/mean-to-eccentric-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AACvD,OAAO,EAAE,mBAAmB,EAAE,MAAM,8BAA8B,CAAC;AAOnE,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAC9D,OAAO,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAEhD,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2CG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,sBAAsB,GAAG,CACpC,IAAuB,EACvB,QAA2B,EAClB,EAAE;IACX,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;IAElC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,yBAAyB,CAAC,6CAA6C,CACxE,CAAC;IACJ,CAAC;IAED,iDAAiD;IACjD,MAAM,YAAY,GAAG,mBAAmB,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAC5D,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IACnD,6DAA6D;IAC7D,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAEzD,IAAI,QAAQ,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,EAAa,CAAC;IACvB,CAAC;IAED,2DAA2D;IAC3D,MAAM,WAAW,GAAG,IAAI,CAAC,CAAC,4BAA4B;IACtD,MAAM,iBAAiB,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,6CAA6C;IAC1F,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,WAAW,EACX,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE,iBAAiB,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CACrE,CAAC;IACF,wEAAwE;IACxE,MAAM,KAAK,GAAG,EAAE,CAAC;IAEjB,IAAI,IAAI,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC,CAAC;IAEtD,IAAI,QAAQ,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAC1B,IAAI,GAAG,KAAK,CAAC;IACf,CAAC;IAED,OAAO,IAAe,CAAC;AACzB,CAAC,CAAC"}