I want to use planetEphemeris to obtain daily noon (X,Y,Z) positions of former planet Pluto wrt to Mercury, starting at midnight in
about Year 0 AD and continuing to about year 4000 AD. An example of the code I have written to do this looks like:
***********************************************************************************************************************************************
for n=1721027.5:1725027 % This for loop uses Julian date info and does about 4000 loops for testing purposes
% Get Pluto's XYZ position relative to Mercury at this time
position_merc = planetEphemeris([n 0.5],'Mercury','Pluto'); % This was line 24 in my original code. "0.5" means get the noon position.
end
************************************************************************************************************************************************
However, right away I get the following error message. No idea what is going on. Am I not using the database correctly? Does the database DE405 only address a limited time period? How do I get (X,Y,Z) data over longer time periods? is there another database I should be using? :
***************************************************************************************************************************************************************
Error using ephRegister
The requested Julian date (1721028) is outside the range of the ephemerides database. For the database DE405, the Julian dates must
be between 2305424.5 and 2525008.5.
Error in planetEphemeris (line 197)
[registerNumber,aufac,t,warnFlag] = ephRegister(ephemerisTime(m,1:end),units,...
Error in test_mercury_pluto (line 24)
position_merc = planetEphemeris([n 0.5],'Mercury','Pluto');
****************************************************************************************************************************************************************
