File Exchange

image thumbnail

Moon Position

version (17.4 MB) by Meysam Mahooti
Position of the Moon referred to the mean equator and equinox of J2000


Updated 14 Nov 2019

View License

Five different approaches are used in the test_MoonPos.m for the computation of lunar coordinates; NASA JPL Development Ephemerides (DE436), very accurate ELP2000-82, high-precision analytical series (Brown's theory), Simpson analytical method and low-precision analytical series.
1. Montenbruck O., Gill E.; Satellite Orbits: Models, Methods and Applications; Springer Verlag, Heidelberg; Corrected 3rd Printing (2005).
2. Montenbruck O., Pfleger T.; Astronomy on the Personal Computer; Springer Verlag, Heidelberg; 4th edition (2000).
3. Vallado D. A; Fundamentals of Astrodynamics and Applications; McGraw-Hill; New York; 3rd edition (2007).
4. van Flandern T. C., Pulkkinen K. F.; Low precision formulae for planetary positions; Astrophysical Journal Supplement Series 41, 391 (1979).

Cite As

Meysam Mahooti (2020). Moon Position (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (11)

Dear Serge,
The coordinate systems in all functions is J2000.
Best regards,


It would be useful to state coordinate systems in every function.
Most functions, including demo say simply "Moon position (km)".
MoonBrown.m says J2000, but is that true for all functions?


Would it be useful to state coordinate systems in every function?
Most functions, including the demo simply say "r_moon Moon position (km)", but does not state coordinate system.

liuang zhou

Xinyuan Mao

Meg Noah

Demi Moore


JPL Development Ephemerides (DE430) is replaced by DE436.

Ephemeris Time (ET) is introduced as the best approximation of Barycentric Dynamical Time (TDB) and Terrestrial Time (TT) for prediction purposes. Moreover, very accurate ELP2000-82 lunar coordinates is computed.


Revised on 2016-12-17.

Description is updated.

MATLAB Release Compatibility
Created with R2019a
Compatible with any release
Platform Compatibility
Windows macOS Linux