Modified Julian Date

The most accurate calculation of modified Julian date (modified Julian date = Julian date - 2400000.5)

For more than six years I used my previous version of Mjday.m and I thought it was the most accurate because its results were the same as juliandate function of MATLAB. Today, I computed JPL ephemerides with my newer version of Mjday.m and I could get to an accuracy of one centimeter for the position of the major planets, the Sun, and the Moon. If you are doing computations in astrodynamics, aerospace, astronomy, and related fields, I strongly recommend you to use my Mjday.m rather than other versions.