Code covered by the BSD License  

Highlights from
A MATLAB Script for Predicting Lunar Eclipses

  • aberatthis function corrects position vector for aberration of light.
  • anglesthis function converts a vector to angular quantities
  • applan1this function computes the apparent geocentric or topocentric place
  • atan3four quadrant inverse tangent
  • brentsolve for a single real root of a nonlinear equation
  • brootbracket a single root of a nonlinear equation
  • deg2dmsconvert decimal degrees to degrees,
  • ecleventpredict lunar eclipse events
  • eclfunc(x)lunar eclipse objective function
  • eclprint(iflag, jdutc)print lunar eclipse conditions
  • etilt1this function computes quantities related to the orientation
  • findleap(jdate)find number of leap seconds for utc julian date
  • funargthis subroutine computes fundamental arguments (mean elements)
  • gast2this function computes the greenwich sidereal time
  • gdateconvert Julian date to Gregorian (calendar) date
  • geocenthis function moves the origin of coordinates from the
  • getdateinteractive request and input of calendar date
  • getobsinteractive request of observer coordinates
  • hrs2hmsconvert decimal hours to hours,
  • jd2str(jdate)convert Julian date to string equivalent
  • jplephemreads the jpl planetary ephemeris and gives
  • julianJulian date
  • minimaone-dimensional minimization
  • nutate1this function nutates equatorial rectangular coordinates from
  • oevent1predict minimization/root-finding orbital events
  • precessthis function precesses equatorial rectangular coordinates from
  • readleapread leap seconds data file
  • solsyspurpose
  • sunfldthis function corrects position vector for the deflection
  • tdtimesthis function computes the terrestrial time (tt)
  • terrathis functions computes the position and velocity vectors of
  • utc2tdtconvert UTC julian date to TDT julian date
  • leclipse.m
  • View all files
Be the first to rate this file! 6 Downloads (last 30 days) File Size: 469 KB File ID: #39344 Version: 1.1

A MATLAB Script for Predicting Lunar Eclipses


David Eagle (view profile)


06 Dec 2012 (Updated )

Predicting the local circumstances of lunar eclipses.

| Watch this File

File Information

PDF document and MATLAB script named leclipse.m that can be used to predict local circumstances of lunar eclipses. This software provides the eclipse type, the universal times and topocentric coordinates of the Moon at the beginning and end of the penumbra contacts, and the time and coordinates at maximum eclipse. The source ephemeris for this routine is a JPL binary ephemeris file. This application uses several functions ported to MATLAB from the Fortran version of the NOVAS (Naval Observatory Vector Astrometry Subroutines) source code developed at the United States Naval Observatory ( JPL binary ephemeris files for Windows compatible computers can be downloaded at

MATLAB release MATLAB 7.12 (R2011a)
Other requirements Requires platform specific version of MICE mex file and de421.bsp binary ephemeris file. Both are available at
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
04 Feb 2013 1.1

Added MICE version of this script in the archive.

Contact us