Equation of Time

Compute the equation of time with accuracy on the order of seconds. This routine can handle multiple millenniums.

Function Call With Time String:
>> EQtime = EquationOfTime('2000/01/01 00:00:00');

Function Call With Julian Date:
>> EQtime = EquationOfTime(2451544.5);

Input List:
UTCjd (Coordinated Universal Time YYYY/MM/DD hh:mm:ss or Julian Date (Fractional Days)) [N x M] (double or char)

Output List:
EQtime (Equation of Time Solution for a specified date in fractional minutes) [N x 1] (double)