diff --git a/doc/calc-help/constants.md b/doc/calc-help/constants.md
index b08ff592..1cf55d5a 100644
--- a/doc/calc-help/constants.md
+++ b/doc/calc-help/constants.md
@@ -1115,7 +1115,7 @@ the Sun. It is distinct from the rotation period `Prot☿` (the spin period).
 Computed from Meeus' Chapter 38 orbital elements and cross-checked against JPL
 Horizons. The stored value is a mean orbital period, used as an approximation of
 the anomalistic period that governs the recurrence of the perihelion passages
-exposed by `T₀☿`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+exposed by `T₀☿`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀☿ constant
 
@@ -1126,7 +1126,7 @@ reference perihelion passage (Tp) by whole orbital periods (`Porb☿`) and retur
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`, so the constant displays as a date but tracks "now". Mercury passes
-perihelion roughly four times per year. [Materials 22](#materials-22)
+perihelion roughly four times per year. [Reference 27](#reference-27)
 
 
 ## Venus constants
@@ -1259,7 +1259,7 @@ Sun. It is distinct from the rotation period `Prot♀` (the spin period). Comput
 from Meeus' Chapter 38 orbital elements and cross-checked against JPL Horizons.
 The stored value is a mean orbital period, used as an approximation of the
 anomalistic period that governs the recurrence of the perihelion passages
-exposed by `T₀♀`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+exposed by `T₀♀`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♀ constant
 
@@ -1269,7 +1269,7 @@ Computed, not stored: `T₀♀` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♀`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Venus passes perihelion about twice per year. [Materials 22](#materials-22)
+`JDN→`. Venus passes perihelion about twice per year. [Reference 27](#reference-27)
 
 
 ## Earth constants
@@ -1405,7 +1405,7 @@ Earth's sidereal orbital period — the time to complete one revolution around t
 Sun. It is distinct from the rotation period `Prot♁` (the sidereal day). Computed
 from Meeus' Chapter 38 orbital elements and cross-checked against JPL Horizons.
 The stored value is the mean (anomalistic) year used to advance `T₀♁` between
-successive perihelion passages. [Materials 22](#materials-22) [Reference 23](#reference-23)
+successive perihelion passages. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♁ constant
 
@@ -1415,7 +1415,7 @@ Computed, not stored: `T₀♁` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♁`, the
 anomalistic year) and returns the most recent perihelion passage at or before
 the current date (early January each year). The computation is carried out in
-Julian Day Number and converted to a date with `JDN→`. [Materials 22](#materials-22)
+Julian Day Number and converted to a date with `JDN→`. [Reference 27](#reference-27)
 
 
 ### a♁GPS constant
@@ -1617,7 +1617,7 @@ passages (about 27.55 days). It is distinct from the rotation period `Prot☽` a
 from the sidereal month. Its absolute uncertainty (±1.12 days) is the measured 1σ
 spread of real perigee-to-perigee intervals over 2000–2050, made large by the
 solar perturbation; this variability is why the Moon differs from the planets.
-The value is the standard anomalistic month used by `ⓁPeriSel`. [Materials 22](#materials-22)
+The value is the standard anomalistic month used by `ⓁPeriSel`. [Reference 27](#reference-27)
 
 ### T₀☽ constant
 
@@ -1636,7 +1636,7 @@ accepts either a date or a Julian Day Number as input. The perigee recurs about
 every 27.55 days (`Porb☽`). Validated against observed perigees: JDN 2457706.5 →
 2457706.974818 (2016-Nov-14, k=224) and 2464291.5 → 2464292.421948 (2034-Nov-25,
 k=463). Note: Meeus' worked example at p.357 (JDN 2447442.35) is an apogee, not a
-perigee, and must not be used to check this value. [Materials 22](#materials-22)
+perigee, and must not be used to check this value. [Reference 27](#reference-27)
 
 
 ## Mars constants
@@ -1766,7 +1766,7 @@ Sun. It is distinct from the rotation period `Prot♂` (the spin period, close t
 an Earth day). Computed from Meeus' Chapter 38 orbital elements and cross-checked
 against JPL Horizons. The stored value is a mean orbital period, used as an
 approximation of the anomalistic period that governs the recurrence of the
-perihelion passages exposed by `T₀♂`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+perihelion passages exposed by `T₀♂`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♂ constant
 
@@ -1776,7 +1776,7 @@ Computed, not stored: `T₀♂` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♂`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Mars's orbital period is approximately 1.88 years. [Materials 22](#materials-22)
+`JDN→`. Mars's orbital period is approximately 1.88 years. [Reference 27](#reference-27)
 
 
 ## Jupiter constants
@@ -1911,7 +1911,7 @@ Jupiter's orbital period — the time to complete one revolution around the Sun
 (about 11.86 years). It is distinct from the rotation period `Prot♃` (the ~10 h
 spin period). The value is the anomalistic period (perihelion to perihelion),
 computed from Meeus' Chapter 38, which governs the recurrence of the perihelion
-passages exposed by `T₀♃`. [Materials 22](#materials-22)
+passages exposed by `T₀♃`. [Reference 27](#reference-27)
 
 ### T₀♃ constant
 
@@ -1921,7 +1921,7 @@ Computed, not stored: `T₀♃` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♃`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Jupiter's orbital period is approximately 11.86 years. [Materials 22](#materials-22)
+`JDN→`. Jupiter's orbital period is approximately 11.86 years. [Reference 27](#reference-27)
 
 
 ## Saturn constants
@@ -2055,7 +2055,7 @@ Saturn's orbital period — the time to complete one revolution around the Sun
 (about 29.46 years). It is distinct from the rotation period `Prot♄` (the ~10.7 h
 spin period). The value is the anomalistic period (perihelion to perihelion),
 computed from Meeus' Chapter 38, which governs the recurrence of the perihelion
-passages exposed by `T₀♄`. [Materials 22](#materials-22)
+passages exposed by `T₀♄`. [Reference 27](#reference-27)
 
 ### T₀♄ constant
 
@@ -2065,7 +2065,7 @@ Computed, not stored: `T₀♄` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♄`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Saturn's orbital period is approximately 29.46 years. [Materials 22](#materials-22)
+`JDN→`. Saturn's orbital period is approximately 29.46 years. [Reference 27](#reference-27)
 
 
 ## Uranus constants
@@ -2200,7 +2200,7 @@ Uranus's orbital period — the time to complete one revolution around the Sun
 (about 84 years). It is distinct from the rotation period `Prot⛢` (the ~17 h spin
 period). The value is the anomalistic period (perihelion to perihelion), computed
 from Meeus' Chapter 38, which governs the recurrence of the perihelion passages
-exposed by `T₀⛢`. [Materials 22](#materials-22)
+exposed by `T₀⛢`. [Reference 27](#reference-27)
 
 ### T₀⛢ constant
 
@@ -2212,7 +2212,7 @@ the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`. Uranus's orbital period is approximately 84 years; the last perihelion
 was in 1966 and the next is around 2050, so the returned value can be decades in
-the past. [Materials 22](#materials-22)
+the past. [Reference 27](#reference-27)
 
 
 ## Neptune constants
@@ -2344,7 +2344,7 @@ Neptune's sidereal orbital period — the time to complete one revolution around
 the Sun (about 165 years). It is distinct from the rotation period `Prot♆` (the
 ~16 h spin period). Tabulated sidereal value from the NASA Planetary Fact Sheet,
 used as an approximation of the anomalistic period that governs the recurrence
-of the perihelion passages exposed by `T₀♆`. [Materials 23](#materials-23)
+of the perihelion passages exposed by `T₀♆`. [Reference 28](#reference-28)
 
 ### T₀♆ constant
 
@@ -2357,7 +2357,7 @@ Julian Day Number and converted to a date with `JDN→`. Neptune is a deliberate
 exception to the "most recent past perihelion" convention: its true last
 perihelion (~1876) is uninformative and hard to source, so Tp is set to the next
 perihelion (2042-09-04). Because that date is in the future, the floor-IFTE
-returns it unchanged, giving the useful upcoming date. [Materials 21](#materials-21)
+returns it unchanged, giving the useful upcoming date. [Reference 4](#reference-4)
 
 
 ## Pluto constants
@@ -2489,7 +2489,7 @@ Pluto's sidereal orbital period — the time to complete one revolution around t
 Sun (about 248 years). It is distinct from the rotation period `Prot♇` (the
 ~6.4 day spin period). Tabulated sidereal value from the NASA Planetary Fact
 Sheet, used as an approximation of the anomalistic period that governs the
-recurrence of the perihelion passages exposed by `T₀♇`. [Materials 23](#materials-23)
+recurrence of the perihelion passages exposed by `T₀♇`. [Reference 28](#reference-28)
 
 ### T₀♇ constant
 
@@ -2501,7 +2501,7 @@ the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`. Pluto's orbital period is approximately 248 years; the last perihelion
 was in 1989 (Sep 05) and the next is around 2237, so the returned value can be
-decades in the past. [Materials 21](#materials-21)
+decades in the past. [Reference 4](#reference-4)
 
 ## Sun constants
 
@@ -2666,7 +2666,7 @@ Milky Way dark matter fraction
 
 Measured fraction of dark matter in the Milky Way galaxy. Estimated
 from rotation curves, satellite dynamics, and gravitational lensing.
-[Reference 8](#reference-8) [Reference 9](#reference-9)
+[Reference 9](#reference-9) [Reference 32](#reference-32)
 
 ### RG constant
 
@@ -2705,21 +2705,21 @@ Andromeda galaxy total mass
 
 Measure of the total mass of the Andromeda galaxy (M31, NGC 224), including
 dark matter, estimated from satellite galaxy kinematics and gravitational
-timing arguments. [Reference 11](#reference-11) [Reference 12](#reference-12)
+timing arguments. [Reference 35](#reference-35) [Reference 36](#reference-36)
 
 ### fDMNGC224 constant
 
 Andromeda dark matter fraction
 
 Measured fraction of dark matter in the Andromeda galaxy (M31), estimated
-from dynamical modelling of its halo. [Reference 11](#reference-11) [Reference 12](#reference-12)
+from dynamical modelling of its halo. [Reference 31](#reference-31) [Reference 35](#reference-35)
 
 ### RNGC224 constant
 
 Andromeda galaxy disk radius
 
 Measurement of the mean equatorial radius of the Andromeda galaxy (M31),
-estimated from stellar population surveys. [Reference 11](#reference-11)
+estimated from stellar population surveys. [Reference 34](#reference-34)
 
 ### DNGC224 constant
 
@@ -2727,7 +2727,7 @@ Distance to Andromeda galaxy
 
 Measurement of the mean distance from the Milky Way to the Andromeda
 galaxy (M31), measured using Cepheid variable stars and tip of the red
-giant branch. [Reference 11](#reference-11)
+giant branch. [Reference 29](#reference-29)
 
 ### GM★Vega constant
 
@@ -3265,7 +3265,7 @@ Dielectric
 Fused silica density
 
 Measured. Volumetric density of amorphous fused silica (vitreous SiO2),
-the primary material for standard optical fibers. [Materials 5](#materials-5) [Materials 7](#materials-7)
+the primary material for standard optical fibers. [Materials 7](#materials-7)
 
 ### nSiO2 constant
 
@@ -3720,14 +3720,6 @@ Watkins, L.L., et al. (2019). "Evidence for an Intermediate-Mass Milky Way from
 
 Posti, L., & Helmi, A. (2019). "Mass and shape of the Milky Way's dark matter halo with globular clusters from Gaia and Hubble". Astronomy & Astrophysics, 621, A56. arXiv:1805.01408 — DOI: 10.1051/0004-6361/201833355
 
-### Reference 11
-
-Peñarrubia, J., et al. (2014). "A timing constraint on the (total) mass of the Large Magellanic Cloud". Monthly Notices of the Royal Astronomical Society, 443(3), 2204–2222. arXiv:1405.4662 — DOI: 10.1093/mnras/stu879
-
-### Reference 12
-
-Veljanoski, J., et al. (2014). "The M31 satellite plane: Evidence for a recent flyby?" Monthly Notices of the Royal Astronomical Society: Letters, 442(1), L86–L90. arXiv:1404.5988 — DOI: 10.1093/mnrasl/slu050
-
 ### Reference 13
 
 Yoon, J., et al. (2010). "A New View of Vega's Composition, Mass, and Age". The Astrophysical Journal, 708(1), 71–79. arXiv:0910.3124 — DOI: 10.1088/0004-637X/708/1/71
@@ -3736,6 +3728,30 @@ Yoon, J., et al. (2010). "A New View of Vega's Composition, Mass, and Age". The
 
 Aufdenberg, J.P., et al. (2006). "First Results from the CHARA Array. VII. Long-Baseline Interferometric Measurements of Vega Consistent with a Pole-On, Rapidly Rotating Star". The Astrophysical Journal, 645(1), 664–675. arXiv:astro-ph/0603327 — DOI: 10.1086/504149
 
+### Reference 29
+
+Li, S., Riess, A.G., & Yuan, W. (2021). The Astrophysical Journal, 920, 84. DOI: 10.3847/1538-4357/ac1597 — Cepheid distance to M31 (761 ± 11 kpc), source for the distance to the Andromeda galaxy `DNGC224`.
+
+### Reference 31
+
+Sick, J., Courteau, S., Cuillandre, J.-C., et al. (2015). IAU Symposium 311, 82 (arXiv:1410.0017) — M31 stellar mass and mass-to-light modelling, source for the Andromeda dark matter fraction `fDMNGC224`.
+
+### Reference 32
+
+Licquia, T.C., & Newman, J.A. (2015). The Astrophysical Journal, 806, 96. DOI: 10.1088/0004-637X/806/1/96 — Milky Way stellar mass, source for the Milky Way dark matter fraction `fDMG`.
+
+### Reference 34
+
+Nieten, C., Neininger, N., Guélin, M., et al. (2006). Astronomy & Astrophysics, 453, 459. DOI: 10.1051/0004-6361:20035672 — M31 disk extent from CO and HI surveys, source for the Andromeda galaxy disk radius `RNGC224`.
+
+### Reference 35
+
+Watkins, L.L., Evans, N.W., & An, J.H. (2010). Monthly Notices of the Royal Astronomical Society, 406, 264. DOI: 10.1111/j.1365-2966.2010.16708.x — Tracer mass estimators applied to the Andromeda satellite system, source for the Andromeda galaxy total mass `MNGC224` and the dynamical mass in `fDMNGC224`.
+
+### Reference 36
+
+Zhang, X., Chen, B., Chen, P., Sun, J., & Tian, Z. (2024). Monthly Notices of the Royal Astronomical Society, 528, 2653. DOI: 10.1093/mnras/stae025 — Updated dynamical mass of the Andromeda galaxy, corroborating `MNGC224`.
+
 ### Particle physics data
 
 ### Reference 15
@@ -3790,6 +3806,14 @@ Williams, J.G., et al. (2014). "Lunar interior properties from the GRAIL mission
 
 Nimmo, F., et al. (2017). "Mean radius and shape of Pluto and Charon from New Horizons images". Icarus, 287, 12–29. arXiv:1603.00821 — DOI: 10.1016/j.icarus.2016.06.027
 
+### Reference 27
+
+Meeus, J. (1998). Astronomical Algorithms, 2nd ed. Willmann-Bell Inc., Richmond, Virginia. ISBN: 978-0-943396-61-3. Chapter 38: Perihelion and Aphelion of the Planets; Chapter 50: Perigee and Apogee of the Moon, pp. 355–358. (Planetary orbital periods derived from the Chapter 38 elements; lunar perigee formula from Chapter 50.)
+
+### Reference 28
+
+Williams, D.R. NASA Planetary Fact Sheet. NASA Space Science Data Coordinated Archive (NSSDCA), Goddard Space Flight Center. [Source](https://nssdc.gsfc.nasa.gov/planetary/factsheet/) Accessed: June 2026. (Sidereal orbital periods of the giant planets, used as an approximation of the anomalistic period.)
+
 ## Materials & Fluids data
 
 ### Materials 1
@@ -3812,10 +3836,6 @@ Wagner, W., & Pruss, A. (2002). "The IAPWS Formulation 1995 for the Thermodynami
 
 Malitson, I.H. (1965). "Interspecimen Comparison of the Refractive Index of Fused Silica". J. Opt. Soc. Am., 55(10), 1205–1209. DOI: 10.1364/JOSA.55.001205
 
-### Materials 6
-
-Palik, E.D. (Ed.) (1985). Handbook of Optical Constants of Solids. Academic Press, New York. ISBN: 0-12-544420-6
-
 ### Materials 7
 
 Corning Incorporated (2023). Corning HPFS Fused Silica — Standard Grade Optical Properties. Product datasheet. [Source](https://www.corning.com)
@@ -3866,19 +3886,3 @@ Luke, K., Okawachi, Y., Lamont, M.R.E., Gaeta, A.L., & Lipson, M. (2015). "Broad
 ### Materials 19
 
 NIST Materials Measurement Laboratory (2024). Materials Properties Database. [Source](https://trc.nist.gov) — Accessed: March 2026
-
-### Materials 20
-
-Espenak, F. (2025–2026). Sky Event Almanacs — Greenwich Mean Time. AstroPixels.com. [Source](https://www.astropixels.com/almanac/almanac21) Accessed: April–May 2026. (Source for perihelion/perigee dates of inner planets and Moon.)
-
-### Materials 21
-
-Park, R.S., et al. (2021). "The JPL Planetary and Lunar Ephemerides DE440 and DE441". The Astronomical Journal, 161(3), 105. DOI: 10.3847/1538-3881/abd414 (Source for perihelion dates of outer planets: Jupiter to Pluto.)
-
-### Materials 22
-
-Meeus, J. (1998). Astronomical Algorithms, 2nd ed. Willmann-Bell Inc., Richmond, Virginia. ISBN: 978-0-943396-61-3. Chapter 38: Perihelion and Aphelion of the Planets; Chapter 50: Perigee and Apogee of the Moon, pp. 355–358. (Planetary orbital periods derived from the Chapter 38 elements; lunar perigee formula from Chapter 50.)
-
-### Materials 23
-
-Williams, D.R. NASA Planetary Fact Sheet. NASA Space Science Data Coordinated Archive (NSSDCA), Goddard Space Flight Center. [Source](https://nssdc.gsfc.nasa.gov/planetary/factsheet/) Accessed: June 2026. (Sidereal orbital periods of the giant planets, used as an approximation of the anomalistic period.)
diff --git a/help/db48x.md b/help/db48x.md
index 54de6f61..9e27e2a4 100644
--- a/help/db48x.md
+++ b/help/db48x.md
@@ -7418,7 +7418,7 @@ the Sun. It is distinct from the rotation period `Prot☿` (the spin period).
 Computed from Meeus' Chapter 38 orbital elements and cross-checked against JPL
 Horizons. The stored value is a mean orbital period, used as an approximation of
 the anomalistic period that governs the recurrence of the perihelion passages
-exposed by `T₀☿`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+exposed by `T₀☿`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀☿ constant
 
@@ -7429,7 +7429,7 @@ reference perihelion passage (Tp) by whole orbital periods (`Porb☿`) and retur
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`, so the constant displays as a date but tracks "now". Mercury passes
-perihelion roughly four times per year. [Materials 22](#materials-22)
+perihelion roughly four times per year. [Reference 27](#reference-27)
 
 
 ## Venus constants
@@ -7562,7 +7562,7 @@ Sun. It is distinct from the rotation period `Prot♀` (the spin period). Comput
 from Meeus' Chapter 38 orbital elements and cross-checked against JPL Horizons.
 The stored value is a mean orbital period, used as an approximation of the
 anomalistic period that governs the recurrence of the perihelion passages
-exposed by `T₀♀`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+exposed by `T₀♀`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♀ constant
 
@@ -7572,7 +7572,7 @@ Computed, not stored: `T₀♀` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♀`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Venus passes perihelion about twice per year. [Materials 22](#materials-22)
+`JDN→`. Venus passes perihelion about twice per year. [Reference 27](#reference-27)
 
 
 ## Earth constants
@@ -7708,7 +7708,7 @@ Earth's sidereal orbital period — the time to complete one revolution around t
 Sun. It is distinct from the rotation period `Prot♁` (the sidereal day). Computed
 from Meeus' Chapter 38 orbital elements and cross-checked against JPL Horizons.
 The stored value is the mean (anomalistic) year used to advance `T₀♁` between
-successive perihelion passages. [Materials 22](#materials-22) [Reference 23](#reference-23)
+successive perihelion passages. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♁ constant
 
@@ -7718,7 +7718,7 @@ Computed, not stored: `T₀♁` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♁`, the
 anomalistic year) and returns the most recent perihelion passage at or before
 the current date (early January each year). The computation is carried out in
-Julian Day Number and converted to a date with `JDN→`. [Materials 22](#materials-22)
+Julian Day Number and converted to a date with `JDN→`. [Reference 27](#reference-27)
 
 
 ### a♁GPS constant
@@ -7920,7 +7920,7 @@ passages (about 27.55 days). It is distinct from the rotation period `Prot☽` a
 from the sidereal month. Its absolute uncertainty (±1.12 days) is the measured 1σ
 spread of real perigee-to-perigee intervals over 2000–2050, made large by the
 solar perturbation; this variability is why the Moon differs from the planets.
-The value is the standard anomalistic month used by `ⓁPeriSel`. [Materials 22](#materials-22)
+The value is the standard anomalistic month used by `ⓁPeriSel`. [Reference 27](#reference-27)
 
 ### T₀☽ constant
 
@@ -7939,7 +7939,7 @@ accepts either a date or a Julian Day Number as input. The perigee recurs about
 every 27.55 days (`Porb☽`). Validated against observed perigees: JDN 2457706.5 →
 2457706.974818 (2016-Nov-14, k=224) and 2464291.5 → 2464292.421948 (2034-Nov-25,
 k=463). Note: Meeus' worked example at p.357 (JDN 2447442.35) is an apogee, not a
-perigee, and must not be used to check this value. [Materials 22](#materials-22)
+perigee, and must not be used to check this value. [Reference 27](#reference-27)
 
 
 ## Mars constants
@@ -8069,7 +8069,7 @@ Sun. It is distinct from the rotation period `Prot♂` (the spin period, close t
 an Earth day). Computed from Meeus' Chapter 38 orbital elements and cross-checked
 against JPL Horizons. The stored value is a mean orbital period, used as an
 approximation of the anomalistic period that governs the recurrence of the
-perihelion passages exposed by `T₀♂`. [Materials 22](#materials-22) [Reference 23](#reference-23)
+perihelion passages exposed by `T₀♂`. [Reference 27](#reference-27) [Reference 23](#reference-23)
 
 ### T₀♂ constant
 
@@ -8079,7 +8079,7 @@ Computed, not stored: `T₀♂` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♂`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Mars's orbital period is approximately 1.88 years. [Materials 22](#materials-22)
+`JDN→`. Mars's orbital period is approximately 1.88 years. [Reference 27](#reference-27)
 
 
 ## Jupiter constants
@@ -8214,7 +8214,7 @@ Jupiter's orbital period — the time to complete one revolution around the Sun
 (about 11.86 years). It is distinct from the rotation period `Prot♃` (the ~10 h
 spin period). The value is the anomalistic period (perihelion to perihelion),
 computed from Meeus' Chapter 38, which governs the recurrence of the perihelion
-passages exposed by `T₀♃`. [Materials 22](#materials-22)
+passages exposed by `T₀♃`. [Reference 27](#reference-27)
 
 ### T₀♃ constant
 
@@ -8224,7 +8224,7 @@ Computed, not stored: `T₀♃` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♃`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Jupiter's orbital period is approximately 11.86 years. [Materials 22](#materials-22)
+`JDN→`. Jupiter's orbital period is approximately 11.86 years. [Reference 27](#reference-27)
 
 
 ## Saturn constants
@@ -8358,7 +8358,7 @@ Saturn's orbital period — the time to complete one revolution around the Sun
 (about 29.46 years). It is distinct from the rotation period `Prot♄` (the ~10.7 h
 spin period). The value is the anomalistic period (perihelion to perihelion),
 computed from Meeus' Chapter 38, which governs the recurrence of the perihelion
-passages exposed by `T₀♄`. [Materials 22](#materials-22)
+passages exposed by `T₀♄`. [Reference 27](#reference-27)
 
 ### T₀♄ constant
 
@@ -8368,7 +8368,7 @@ Computed, not stored: `T₀♄` evaluates an expression (IFTE) that advances fro
 reference perihelion passage (Tp) by whole orbital periods (`Porb♄`) and returns
 the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
-`JDN→`. Saturn's orbital period is approximately 29.46 years. [Materials 22](#materials-22)
+`JDN→`. Saturn's orbital period is approximately 29.46 years. [Reference 27](#reference-27)
 
 
 ## Uranus constants
@@ -8503,7 +8503,7 @@ Uranus's orbital period — the time to complete one revolution around the Sun
 (about 84 years). It is distinct from the rotation period `Prot⛢` (the ~17 h spin
 period). The value is the anomalistic period (perihelion to perihelion), computed
 from Meeus' Chapter 38, which governs the recurrence of the perihelion passages
-exposed by `T₀⛢`. [Materials 22](#materials-22)
+exposed by `T₀⛢`. [Reference 27](#reference-27)
 
 ### T₀⛢ constant
 
@@ -8515,7 +8515,7 @@ the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`. Uranus's orbital period is approximately 84 years; the last perihelion
 was in 1966 and the next is around 2050, so the returned value can be decades in
-the past. [Materials 22](#materials-22)
+the past. [Reference 27](#reference-27)
 
 
 ## Neptune constants
@@ -8647,7 +8647,7 @@ Neptune's sidereal orbital period — the time to complete one revolution around
 the Sun (about 165 years). It is distinct from the rotation period `Prot♆` (the
 ~16 h spin period). Tabulated sidereal value from the NASA Planetary Fact Sheet,
 used as an approximation of the anomalistic period that governs the recurrence
-of the perihelion passages exposed by `T₀♆`. [Materials 23](#materials-23)
+of the perihelion passages exposed by `T₀♆`. [Reference 28](#reference-28)
 
 ### T₀♆ constant
 
@@ -8660,7 +8660,7 @@ Julian Day Number and converted to a date with `JDN→`. Neptune is a deliberate
 exception to the "most recent past perihelion" convention: its true last
 perihelion (~1876) is uninformative and hard to source, so Tp is set to the next
 perihelion (2042-09-04). Because that date is in the future, the floor-IFTE
-returns it unchanged, giving the useful upcoming date. [Materials 21](#materials-21)
+returns it unchanged, giving the useful upcoming date. [Reference 4](#reference-4)
 
 
 ## Pluto constants
@@ -8792,7 +8792,7 @@ Pluto's sidereal orbital period — the time to complete one revolution around t
 Sun (about 248 years). It is distinct from the rotation period `Prot♇` (the
 ~6.4 day spin period). Tabulated sidereal value from the NASA Planetary Fact
 Sheet, used as an approximation of the anomalistic period that governs the
-recurrence of the perihelion passages exposed by `T₀♇`. [Materials 23](#materials-23)
+recurrence of the perihelion passages exposed by `T₀♇`. [Reference 28](#reference-28)
 
 ### T₀♇ constant
 
@@ -8804,7 +8804,7 @@ the most recent perihelion passage at or before the current date. The
 computation is carried out in Julian Day Number and converted to a date with
 `JDN→`. Pluto's orbital period is approximately 248 years; the last perihelion
 was in 1989 (Sep 05) and the next is around 2237, so the returned value can be
-decades in the past. [Materials 21](#materials-21)
+decades in the past. [Reference 4](#reference-4)
 
 ## Sun constants
 
@@ -8969,7 +8969,7 @@ Milky Way dark matter fraction
 
 Measured fraction of dark matter in the Milky Way galaxy. Estimated
 from rotation curves, satellite dynamics, and gravitational lensing.
-[Reference 8](#reference-8) [Reference 9](#reference-9)
+[Reference 9](#reference-9) [Reference 32](#reference-32)
 
 ### RG constant
 
@@ -9008,21 +9008,21 @@ Andromeda galaxy total mass
 
 Measure of the total mass of the Andromeda galaxy (M31, NGC 224), including
 dark matter, estimated from satellite galaxy kinematics and gravitational
-timing arguments. [Reference 11](#reference-11) [Reference 12](#reference-12)
+timing arguments. [Reference 35](#reference-35) [Reference 36](#reference-36)
 
 ### fDMNGC224 constant
 
 Andromeda dark matter fraction
 
 Measured fraction of dark matter in the Andromeda galaxy (M31), estimated
-from dynamical modelling of its halo. [Reference 11](#reference-11) [Reference 12](#reference-12)
+from dynamical modelling of its halo. [Reference 31](#reference-31) [Reference 35](#reference-35)
 
 ### RNGC224 constant
 
 Andromeda galaxy disk radius
 
 Measurement of the mean equatorial radius of the Andromeda galaxy (M31),
-estimated from stellar population surveys. [Reference 11](#reference-11)
+estimated from stellar population surveys. [Reference 34](#reference-34)
 
 ### DNGC224 constant
 
@@ -9030,7 +9030,7 @@ Distance to Andromeda galaxy
 
 Measurement of the mean distance from the Milky Way to the Andromeda
 galaxy (M31), measured using Cepheid variable stars and tip of the red
-giant branch. [Reference 11](#reference-11)
+giant branch. [Reference 29](#reference-29)
 
 ### GM★Vega constant
 
@@ -9568,7 +9568,7 @@ Dielectric
 Fused silica density
 
 Measured. Volumetric density of amorphous fused silica (vitreous SiO2),
-the primary material for standard optical fibers. [Materials 5](#materials-5) [Materials 7](#materials-7)
+the primary material for standard optical fibers. [Materials 7](#materials-7)
 
 ### nSiO2 constant
 
@@ -10023,14 +10023,6 @@ Watkins, L.L., et al. (2019). "Evidence for an Intermediate-Mass Milky Way from
 
 Posti, L., & Helmi, A. (2019). "Mass and shape of the Milky Way's dark matter halo with globular clusters from Gaia and Hubble". Astronomy & Astrophysics, 621, A56. arXiv:1805.01408 — DOI: 10.1051/0004-6361/201833355
 
-### Reference 11
-
-Peñarrubia, J., et al. (2014). "A timing constraint on the (total) mass of the Large Magellanic Cloud". Monthly Notices of the Royal Astronomical Society, 443(3), 2204–2222. arXiv:1405.4662 — DOI: 10.1093/mnras/stu879
-
-### Reference 12
-
-Veljanoski, J., et al. (2014). "The M31 satellite plane: Evidence for a recent flyby?" Monthly Notices of the Royal Astronomical Society: Letters, 442(1), L86–L90. arXiv:1404.5988 — DOI: 10.1093/mnrasl/slu050
-
 ### Reference 13
 
 Yoon, J., et al. (2010). "A New View of Vega's Composition, Mass, and Age". The Astrophysical Journal, 708(1), 71–79. arXiv:0910.3124 — DOI: 10.1088/0004-637X/708/1/71
@@ -10039,6 +10031,30 @@ Yoon, J., et al. (2010). "A New View of Vega's Composition, Mass, and Age". The
 
 Aufdenberg, J.P., et al. (2006). "First Results from the CHARA Array. VII. Long-Baseline Interferometric Measurements of Vega Consistent with a Pole-On, Rapidly Rotating Star". The Astrophysical Journal, 645(1), 664–675. arXiv:astro-ph/0603327 — DOI: 10.1086/504149
 
+### Reference 29
+
+Li, S., Riess, A.G., & Yuan, W. (2021). The Astrophysical Journal, 920, 84. DOI: 10.3847/1538-4357/ac1597 — Cepheid distance to M31 (761 ± 11 kpc), source for the distance to the Andromeda galaxy `DNGC224`.
+
+### Reference 31
+
+Sick, J., Courteau, S., Cuillandre, J.-C., et al. (2015). IAU Symposium 311, 82 (arXiv:1410.0017) — M31 stellar mass and mass-to-light modelling, source for the Andromeda dark matter fraction `fDMNGC224`.
+
+### Reference 32
+
+Licquia, T.C., & Newman, J.A. (2015). The Astrophysical Journal, 806, 96. DOI: 10.1088/0004-637X/806/1/96 — Milky Way stellar mass, source for the Milky Way dark matter fraction `fDMG`.
+
+### Reference 34
+
+Nieten, C., Neininger, N., Guélin, M., et al. (2006). Astronomy & Astrophysics, 453, 459. DOI: 10.1051/0004-6361:20035672 — M31 disk extent from CO and HI surveys, source for the Andromeda galaxy disk radius `RNGC224`.
+
+### Reference 35
+
+Watkins, L.L., Evans, N.W., & An, J.H. (2010). Monthly Notices of the Royal Astronomical Society, 406, 264. DOI: 10.1111/j.1365-2966.2010.16708.x — Tracer mass estimators applied to the Andromeda satellite system, source for the Andromeda galaxy total mass `MNGC224` and the dynamical mass in `fDMNGC224`.
+
+### Reference 36
+
+Zhang, X., Chen, B., Chen, P., Sun, J., & Tian, Z. (2024). Monthly Notices of the Royal Astronomical Society, 528, 2653. DOI: 10.1093/mnras/stae025 — Updated dynamical mass of the Andromeda galaxy, corroborating `MNGC224`.
+
 ### Particle physics data
 
 ### Reference 15
@@ -10093,6 +10109,14 @@ Williams, J.G., et al. (2014). "Lunar interior properties from the GRAIL mission
 
 Nimmo, F., et al. (2017). "Mean radius and shape of Pluto and Charon from New Horizons images". Icarus, 287, 12–29. arXiv:1603.00821 — DOI: 10.1016/j.icarus.2016.06.027
 
+### Reference 27
+
+Meeus, J. (1998). Astronomical Algorithms, 2nd ed. Willmann-Bell Inc., Richmond, Virginia. ISBN: 978-0-943396-61-3. Chapter 38: Perihelion and Aphelion of the Planets; Chapter 50: Perigee and Apogee of the Moon, pp. 355–358. (Planetary orbital periods derived from the Chapter 38 elements; lunar perigee formula from Chapter 50.)
+
+### Reference 28
+
+Williams, D.R. NASA Planetary Fact Sheet. NASA Space Science Data Coordinated Archive (NSSDCA), Goddard Space Flight Center. [Source](https://nssdc.gsfc.nasa.gov/planetary/factsheet/) Accessed: June 2026. (Sidereal orbital periods of the giant planets, used as an approximation of the anomalistic period.)
+
 ## Materials & Fluids data
 
 ### Materials 1
@@ -10115,10 +10139,6 @@ Wagner, W., & Pruss, A. (2002). "The IAPWS Formulation 1995 for the Thermodynami
 
 Malitson, I.H. (1965). "Interspecimen Comparison of the Refractive Index of Fused Silica". J. Opt. Soc. Am., 55(10), 1205–1209. DOI: 10.1364/JOSA.55.001205
 
-### Materials 6
-
-Palik, E.D. (Ed.) (1985). Handbook of Optical Constants of Solids. Academic Press, New York. ISBN: 0-12-544420-6
-
 ### Materials 7
 
 Corning Incorporated (2023). Corning HPFS Fused Silica — Standard Grade Optical Properties. Product datasheet. [Source](https://www.corning.com)
@@ -10169,22 +10189,6 @@ Luke, K., Okawachi, Y., Lamont, M.R.E., Gaeta, A.L., & Lipson, M. (2015). "Broad
 ### Materials 19
 
 NIST Materials Measurement Laboratory (2024). Materials Properties Database. [Source](https://trc.nist.gov) — Accessed: March 2026
-
-### Materials 20
-
-Espenak, F. (2025–2026). Sky Event Almanacs — Greenwich Mean Time. AstroPixels.com. [Source](https://www.astropixels.com/almanac/almanac21) Accessed: April–May 2026. (Source for perihelion/perigee dates of inner planets and Moon.)
-
-### Materials 21
-
-Park, R.S., et al. (2021). "The JPL Planetary and Lunar Ephemerides DE440 and DE441". The Astronomical Journal, 161(3), 105. DOI: 10.3847/1538-3881/abd414 (Source for perihelion dates of outer planets: Jupiter to Pluto.)
-
-### Materials 22
-
-Meeus, J. (1998). Astronomical Algorithms, 2nd ed. Willmann-Bell Inc., Richmond, Virginia. ISBN: 978-0-943396-61-3. Chapter 38: Perihelion and Aphelion of the Planets; Chapter 50: Perigee and Apogee of the Moon, pp. 355–358. (Planetary orbital periods derived from the Chapter 38 elements; lunar perigee formula from Chapter 50.)
-
-### Materials 23
-
-Williams, D.R. NASA Planetary Fact Sheet. NASA Space Science Data Coordinated Archive (NSSDCA), Goddard Space Flight Center. [Source](https://nssdc.gsfc.nasa.gov/planetary/factsheet/) Accessed: June 2026. (Sidereal orbital periods of the giant planets, used as an approximation of the anomalistic period.)
 # Equations library
 
 The DB48X calculator features a library of equations covering mathematics,
diff --git a/src/constants.cc b/src/constants.cc
index 428f25fe..7874dfb6 100644
--- a/src/constants.cc
+++ b/src/constants.cc
@@ -901,7 +901,7 @@ static const cstring basic_constants[] =
                 "  0_S "
                 "  0 "
                 "  7.748E-5_S ]",
-    // *Fermi reduced coupling constant - Measurement [2]
+    // *Fermi reduced coupling constant - Measurement [0]
     "G0F",       "[ 1.1663787E-5_GeV^-2 "
                 "  0.0000006E-5_GeV^-2 "
                 "  'ROUND(UBASE(ABS(ⓈG0F/ⒸG0F));-2)' ]",
@@ -944,7 +944,7 @@ static const cstring basic_constants[] =
     "ΔfCs",     "[ 9192631770_Hz "
                 "  0_Hz "
                 "  0 ]",
-    // *Weak mixing angle - Calculation from measurement [2]
+    // *Weak mixing angle - Calculation from measurement [0]
     "θw",       "[ 'ROUND(CONVERT(ASIN(√(0.22305));1_°);XPON(UVAL(Ⓡθw*ASIN(√(0.22305))))-XPON(UVAL(ASIN(√(0.22305))))-2)' "
                    "'CONVERT(ROUND((ASIN(√(0.22305+0.00023))-ASIN(√(0.22305-0.00023)))/2;-2);1_°)' "
                 "  'ROUND(UBASE(ABS(Ⓢθw/CONVERT(ASIN(√(0.22305));1_°)));-2)' "
@@ -1252,7 +1252,7 @@ static const cstring basic_constants[] =
                 "  0.1_s "
                 "  'ROUND(ⓈProt☿/ⒸProt☿;-2)' "
                 "  5.067E6_s ]",
-    // *Mercury orbital period - Computed [M22]
+    // *Mercury orbital period - Computed [27]
     "Porb☿",    "[ 7600551.8_s "
                 "  0_s "
                 "  0 ]",
@@ -1290,7 +1290,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω☿/Ⓒω☿;-2)' "
                 "  29.12_° ]",
-    // *Mercury time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Mercury time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀☿",      "[ 'JDN→(IFTE((JDN(DateTime)-2461178.915934)/UVAL(CONVERT(ⒸPorb☿;1_d))≥1;2461178.915934+IP((JDN(DateTime)-2461178.915934)/UVAL(CONVERT(ⒸPorb☿;1_d)))*UVAL(CONVERT(ⒸPorb☿;1_d));2461178.915934))' "
                 "  0.0080_d 0 ]",
 
@@ -1339,7 +1339,7 @@ static const cstring basic_constants[] =
                 "  8.64_s "
                 "  'ROUND(ⓈProt♀/ⒸProt♀;-2)' "
                 "  2.100E7_s ]",
-    // *Venus orbital period - Computed [M22]
+    // *Venus orbital period - Computed [27]
     "Porb♀",    "[ 19414071.4_s "
                 "  0_s "
                 "  0 ]",
@@ -1378,7 +1378,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♀/Ⓒω♀;-2)' "
                 "  54.88_° ]",
-    // *Venus time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Venus time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀♀",      "[ 'JDN→(IFTE((JDN(DateTime)-2461175.615653)/UVAL(CONVERT(ⒸPorb♀;1_d))≥1;2461175.615653+IP((JDN(DateTime)-2461175.615653)/UVAL(CONVERT(ⒸPorb♀;1_d)))*UVAL(CONVERT(ⒸPorb♀;1_d));2461175.615653))' "
                 "  0.11_d 0 ]",
 
@@ -1424,7 +1424,7 @@ static const cstring basic_constants[] =
                 "  0.0001_s "
                 "  'ROUND(ⓈProt♁/ⒸProt♁;-2)' "
                 "  8.616E4_s ]",
-    // *Earth orbital period - Computed [M22]
+    // *Earth orbital period - Computed [27]
     "Porb♁",    "[ 31558956.5_s "
                 "  0_s "
                 "  0 ]",
@@ -1463,7 +1463,7 @@ static const cstring basic_constants[] =
                 "  0.00001_° "
                 "  'ROUND(Ⓢω♁/Ⓒω♁;-2)' "
                 "  114.2_° ]",
-    // *Earth time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Earth time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀♁",      "[ 'JDN→(IFTE((JDN(DateTime)-2461044.220333)/UVAL(CONVERT(ⒸPorb♁;1_d))≥1;2461044.220333+IP((JDN(DateTime)-2461044.220333)/UVAL(CONVERT(ⒸPorb♁;1_d)))*UVAL(CONVERT(ⒸPorb♁;1_d));2461044.220333))' "
                 "  0.83_d 0 ]",
 
@@ -1554,7 +1554,7 @@ static const cstring basic_constants[] =
                 "  0.1_s "
                 "  'ROUND(ⓈProt☽/ⒸProt☽;-2)' "
                 "  2.361E6_s ]",
-    // *Moon orbital period - anomalistic month (perigee to perigee) [M22]
+    // *Moon orbital period - anomalistic month (perigee to perigee) [27]
     "Porb☽",    "[ 2380713.11_s "
                 "  97165_s "
                 "  0 ]",
@@ -1593,7 +1593,7 @@ static const cstring basic_constants[] =
                 "  'ROUND(Ⓢω☽/Ⓒω☽;-2)' "
                 "  318.2_° ]",
     // *Moon time of perigee passage - PeriSel (full Meeus ch.50); unc = ch.50
-    //  method accuracy vs true perigee ~5 min (1σ), up to ~30 min rare [M22]
+    //  method accuracy vs true perigee ~5 min (1σ), up to ~30 min rare [27]
     "T₀☽",      "[ 'JDN→(ⓁPeriSel(→Num(JDN(DateTime))))' "
                 "  300_s "
                 "  0 ]",
@@ -1643,7 +1643,7 @@ static const cstring basic_constants[] =
                 "  0.1_s "
                 "  'ROUND(ⓈProt♂/ⒸProt♂;-2)' "
                 "  8.864E4_s ]",
-    // *Mars orbital period - Computed [M22]
+    // *Mars orbital period - Computed [27]
     "Porb♂",    "[ 59356065.6_s "
                 "  0_s "
                 "  0 ]",
@@ -1680,7 +1680,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♂/Ⓒω♂;-2)' "
                 "  286.5_° ]",
-    // *Mars time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Mars time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀♂",      "[ 'JDN→(IFTE((JDN(DateTime)-2461125.798009)/UVAL(CONVERT(ⒸPorb♂;1_d))≥1;2461125.798009+IP((JDN(DateTime)-2461125.798009)/UVAL(CONVERT(ⒸPorb♂;1_d)))*UVAL(CONVERT(ⒸPorb♂;1_d));2461125.798009))' "
                 "  0.16_d 0 ]",
 
@@ -1725,7 +1725,7 @@ static const cstring basic_constants[] =
     "Prot♃",    "[ 35730_s "
                 "  1_s "
                 "  'ROUND(ⓈProt♃/ⒸProt♃;-2)' ]",
-    // *Jupiter orbital period - Computed, anomalistic perigee interval [M22]
+    // *Jupiter orbital period - Computed, anomalistic perigee interval [27]
     "Porb♃",    "[ 374360783_s "
                 "  0_s "
                 "  0 ]",
@@ -1762,7 +1762,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♃/Ⓒω♃;-2)' "
                 "  273.9_° ]",
-    // *Jupiter time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Jupiter time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀♃",      "[ 'JDN→(IFTE((JDN(DateTime)-2459964.991260)/UVAL(CONVERT(ⒸPorb♃;1_d))≥1;2459964.991260+IP((JDN(DateTime)-2459964.991260)/UVAL(CONVERT(ⒸPorb♃;1_d)))*UVAL(CONVERT(ⒸPorb♃;1_d));2459964.991260))' "
                 "  8.9_d 0 ]",
 
@@ -1811,7 +1811,7 @@ static const cstring basic_constants[] =
                 "  50_s "
                 "  'ROUND(ⓈProt♄/ⒸProt♄;-2)' "
                 "  3.836E4_s ]",
-    // *Saturn orbital period - Computed, anomalistic perigee interval [M22]
+    // *Saturn orbital period - Computed, anomalistic perigee interval [27]
     "Porb♄",    "[ 928565359_s "
                 "  0_s "
                 "  0 ]",
@@ -1848,7 +1848,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♄/Ⓒω♄;-2)' "
                 "  339.4_° ]",
-    // *Saturn time of perihelion passage - Computed via IFTE from Tp and Porb [M22]
+    // *Saturn time of perihelion passage - Computed via IFTE from Tp and Porb [27]
     "T₀♄",      "[ 'JDN→(IFTE((JDN(DateTime)-2452847.154242)/UVAL(CONVERT(ⒸPorb♄;1_d))≥1;2452847.154242+IP((JDN(DateTime)-2452847.154242)/UVAL(CONVERT(ⒸPorb♄;1_d)))*UVAL(CONVERT(ⒸPorb♄;1_d));2452847.154242))' "
                 "  17_d 0 ]",
 
@@ -1897,7 +1897,7 @@ static const cstring basic_constants[] =
                 "  10_s "
                 "  'ROUND(ⓈProt⛢/ⒸProt⛢;-2)' "
                 "  6.206E4_s ]",
-    // *Uranus orbital period - Computed, anomalistic perigee interval [M22]
+    // *Uranus orbital period - Computed, anomalistic perigee interval [27]
     "Porb⛢",    "[ 2658520424_s "
                 "  0_s "
                 "  0 ]",
@@ -1934,7 +1934,7 @@ static const cstring basic_constants[] =
                 "  0.000001_° "
                 "  'ROUND(Ⓢω⛢/Ⓒω⛢;-2)' "
                 "  97.00_° ]",
-    // *Uranus time of perihelion passage - Computed via IFTE from Tp and Porb [M22]  (last perihelion 1966; next ~2050)
+    // *Uranus time of perihelion passage - Computed via IFTE from Tp and Porb [27]  (last perihelion 1966; next ~2050)
     "T₀⛢",      "[ 'JDN→(IFTE((JDN(DateTime)-2439264.256742)/UVAL(CONVERT(ⒸPorb⛢;1_d))≥1;2439264.256742+IP((JDN(DateTime)-2439264.256742)/UVAL(CONVERT(ⒸPorb⛢;1_d)))*UVAL(CONVERT(ⒸPorb⛢;1_d));2439264.256742))' "
                 "  40_d 0 ]",
 
@@ -1982,7 +1982,7 @@ static const cstring basic_constants[] =
     "Prot♆",    "[ 58000_s "
                 "  100_s "
                 "  'ROUND(ⓈProt♆/ⒸProt♆;-2)' ]",
-    // *Neptune orbital period - Sidereal [M23]
+    // *Neptune orbital period - Sidereal [28]
     "Porb♆",    "[ 5200416000_s "
                 "  0_s "
                 "  0 ]",
@@ -2019,7 +2019,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♆/Ⓒω♆;-2)' "
                 "  273.2_° ]",
-    // *Neptune time of perihelion passage - Computed via IFTE from Tp and Porb [M21]  (EXCEPTION: Tp = next perihelion 2042; floor-IFTE returns it unchanged)
+    // *Neptune time of perihelion passage - Computed via IFTE from Tp and Porb [4]  (EXCEPTION: Tp = next perihelion 2042; floor-IFTE returns it unchanged)
     "T₀♆",      "[ 'JDN→(IFTE((JDN(DateTime)-2467131.5)/UVAL(CONVERT(ⒸPorb♆;1_d))≥1;2467131.5+IP((JDN(DateTime)-2467131.5)/UVAL(CONVERT(ⒸPorb♆;1_d)))*UVAL(CONVERT(ⒸPorb♆;1_d));2467131.5))' "
                 "  10_d 0 ]",
 
@@ -2067,7 +2067,7 @@ static const cstring basic_constants[] =
                 "  0.1_s "
                 "  'ROUND(ⓈProt♇/ⒸProt♇;-2)' "
                 "  5.519E5_s ]",
-    // *Pluto orbital period - Sidereal [M23]
+    // *Pluto orbital period - Sidereal [28]
     "Porb♇",    "[ 7824384000_s "
                 "  0_s "
                 "  0 ]",
@@ -2105,7 +2105,7 @@ static const cstring basic_constants[] =
                 "  0.001_° "
                 "  'ROUND(Ⓢω♇/Ⓒω♇;-2)' "
                 "  113.8_° ]",
-    // *Pluto time of perihelion passage - Computed via IFTE from Tp and Porb [M21]  (last perihelion 1989)
+    // *Pluto time of perihelion passage - Computed via IFTE from Tp and Porb [4]  (last perihelion 1989)
     "T₀♇",      "[ 'JDN→(IFTE((JDN(DateTime)-2447774.5)/UVAL(CONVERT(ⒸPorb♇;1_d))≥1;2447774.5+IP((JDN(DateTime)-2447774.5)/UVAL(CONVERT(ⒸPorb♇;1_d)))*UVAL(CONVERT(ⒸPorb♇;1_d));2447774.5))' "
                 "  5_d 0 ]",
 
@@ -2202,9 +2202,9 @@ static const cstring basic_constants[] =
     "MG",       "[ 1.5E12_M☉ "
                 "  0.5E12_M☉ "
                 "  'ROUND(ⓈMG/ⒸMG;-2)' ]",
-    // *Milky Way dark matter fraction - Measurement [8] [9]
-    "fDMG",     "[ 0.85 "
-                "  0.05 "
+    // *Milky Way dark matter fraction - Measurement [9] [32]
+    "fDMG",     "[ 0.95 "
+                "  0.02 "
                 "  'ROUND(ⓈfDMG/ⒸfDMG;-2)' ]",
 
     // ------------------------------------------------------------------------
@@ -2225,26 +2225,26 @@ static const cstring basic_constants[] =
     "M●SgrA*",  "[ 4.297E6_M☉ "
                 "  0.012E6_M☉ "
                 "  'ROUND(0.012E6/4.297E6;-2)' ]",
-    // *Andromeda galaxy total mass - Measurement [11] [12]
-    "MNGC224",  "[ 1.5E12_M☉ "
-                "  0.5E12_M☉ "
+    // *Andromeda galaxy total mass - Measurement [35] [36]
+    "MNGC224",  "[ 1.4E12_M☉ "
+                "  0.4E12_M☉ "
                 "  'ROUND(ⓈMNGC224/ⒸMNGC224;-2)' ]",
 
     // ------------------------------------------------------------------------
-    // *Andromeda dark matter fraction - Measurement [11] [12]
-    "fDMNGC224", "[ 0.85 "
-                "  0.05 "
+    // *Andromeda dark matter fraction - Measurement [31] [35]
+    "fDMNGC224", "[ 0.92 "
+                "  0.03 "
                 "  'ROUND(ⓈfDMNGC224/ⒸfDMNGC224;-2)' ]",
-    // *Andromeda galaxy disk radius - Measurement [11]
+    // *Andromeda galaxy disk radius - Measurement [34]
     "RNGC224",  "[ 110000_ly "
                 "  10000_ly "
                 "  'ROUND(ⓈRNGC224/ⒸRNGC224;-2)' "
                 "  1.100E5_ly ]",
-    // *Distance to Andromeda galaxy - Measurement [11]
-    "DNGC224",  "[ 2537000_ly "
-                "  50000_ly "
+    // *Distance to Andromeda galaxy - Measurement [29]
+    "DNGC224",  "[ 2482000_ly "
+                "  36000_ly "
                 "  'ROUND(ⓈDNGC224/ⒸDNGC224;-2)' "
-                "  2.537E6_ly ]",
+                "  2.482E6_ly ]",
     // *Vega gravitational parameter - Measurement [13] [14]
     "GM★Vega",  "[ 1.89E19_m³/s² "
                 "  0.05E19_m³/s² "