View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
A MATLAB Script for Propagating Trajectories from the Earth to the Moon

5.0 | 1 rating Rate this file 15 Downloads (last 30 days) File Size: 861 KB File ID: #43067 Version: 1.1

A MATLAB Script for Propagating Trajectories from the Earth to the Moon


David Eagle (view profile)


14 Aug 2013 (Updated )

Integrate n-body trajectories from the Earth to the Moon. Includes non-spherical Earth gravity.

| Watch this File

File Information

PDF document and a MATLAB script called pprop_lunar that can be used to numerically integrate the geocentric orbital equations of motion of a trajectory from the Earth to the Moon. This scientific simulation begins at a user-defined epoch and geocentric state vector (position and velocity vectors) representing the point of trans-lunar injection (TLI) and ends at (1) closest approach to the Moon, (2) a user-defined Moon-centered (selenocentric) distance, or (3) at a user-defined final epoch.

MATLAB release MATLAB 7.14 (R2012a)
Other requirements Requires JPL binary ephemeris file. Available from Orbital Mechanics with MATLAB area of the Orbital and Celestial Mechanics Website located at
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
14 May 2014 Vincenzo Capuano

I soleved the problem; it didn't work because I didm't dwonload the JPL binary ephemeris file.

Anyway, I have another question; have you any idea about the accuracy in position and velocity of this propagator?

Thank you very much,

14 May 2014 Vincenzo Capuano


I tried to run your script, but it doesn't work:

Error using fread
Invalid file identifier. Use fopen to generate a valid file identifier.

Error in jplephem (line 71)
ttl = fread(fid, 252);

Error in pprop_eqm (line 44)
svsun = jplephem(jdate, 11, 3);

Error in odearguments (line 88)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 114)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in findca (line 33)
[t, ysol, tevent, yevent, ie] = ode45(@pprop_eqm, [0 tend], [ri vi], options);

Error in pprop_lunar (line 101)
[jdtdb_ca, re2sc, ve2sc, rm2sc, vm2sc] = findca(jdtdb_tli, ri, vi);

Can you help me please?

Comment only
28 May 2014 1.1

Includes new option to model the effect of the Moon on the trajectory using non-spherical lunar gravity.

Contact us