The Planrise.m computes planetary positions using JPL_Eph_DE430.m. To compute the rising and setting times of a planet we first determine its (geocentric) equatorial coordinates α and δ for an arbitrary time on the given day. From the declination and the geographical latitude, Planrise then determines whether the planet is continuously above or below the horizon.
Montenbruck O., Pfleger T.; Astronomy on the Personal Computer; Springer Verlag, Heidelberg; 4th edition (2000).
Meeus J.; Astronomical Algorithms; Willmann-Bell; Richmond, Virginia; 2nd edition (1998).
The code works with updated DE430Coeff.mat now for all years including years before 1940 and years after 2080.
DE430Coeff.mat full matrix is now available in the following link:
Use it rather than existing DE430Coeff.mat matrix.
The code not works for year before 1940 or year after 2080. The .mat file return index error for all year before 1940 or after 2080.
How to correct it?
Accuracy is improved using planetary positions with respect to true equator and equinox of date rather than mean equator and equinox of date. (TOD_PlanPos = Nutation*Precession*ICRF_PlanPos)
Condition for computation of Moon's rising/setting times is added.
Planetary positions are computed applying Jet Propulsion Laboratory Development Ephemeris (DE430). CityName.m is added to give the name of city you have entered its geographical latitude and longitude.