moonlr.m

球面天文学（岁差、章动）和天体力学（行星星历表）的原始计算公式、算法和程序。matlab编译通过。主要包括：公历/儒略历转换为儒略日

function an=moonlr
% geocentric moon
%   polar coordinates re mean equinox and ecliptic of date

%   For latitude coefficients, see mlat404.c.


%   S. L. Moshier
%   December, 1996


%              Residuals against JPL ephemeris DE404 (arc seconds)

%First date in file = 1221000.5
%Number of samples = 1053099
%Sampling interval = 1.515625 days

%Peak excursions from these mostly different test points
%were consolidated with the above.  They added .01" to a few
%of the peak readings.
%First date in file = 1221000.50
%Number of samples = 524290.0
%Sampling interval = 3.0 days


% Julian Years             Longitude          Latitude           Distance
%                                                               1 = 1.9 km
%                       Peak  RMS   Ave    Peak  RMS   Ave    Peak  RMS   Ave
% -1369.0 to -1000.0:   0.43  0.07  0.00   0.33  0.05 -0.00   0.18  0.03  0.00 
% -1000.0 to  -500.0:   0.49  0.06 -0.00   0.33  0.04 -0.00   0.18  0.03  0.00 
%  -500.0 to     0.0:   0.48  0.06  0.00   0.32  0.04  0.00   0.15  0.03  0.00 
%     0.0 to   500.0:   0.45  0.05  0.00   0.30  0.04 -0.00   0.17  0.03 -0.00 
%   500.0 to  1000.0:   0.48  0.05 -0.00   0.29  0.04  0.00   0.17  0.03 -0.00 
%  1000.0 to  1500.0:   0.42  0.05 -0.00   0.28  0.04 -0.00   0.16  0.03  0.00 
%  1500.0 to  2000.0:   0.35  0.05 -0.00   0.26  0.04  0.00   0.15  0.03  0.00 
%  2000.0 to  2500.0:   0.39  0.06  0.00   0.25  0.04 -0.00   0.15  0.03 -0.00 
%  2500.0 to  3000.0:   0.44  0.07 -0.00   0.30  0.05 -0.00   0.19  0.03 -0.00 
%  3000.0 to  3000.8:   0.23  0.08 -0.04   0.11  0.04 -0.00   0.08  0.03 -0.00 


tabl=[ 175667,66453,5249,-42,20057,403,-2360,6148,-7644,24646,-1273,9127,-1395,1958,232,-289,-97,553,69,130,-80,6,129,-868,26,-89,1042,1172,194,-112,-47433,-241666,224626,-103752,63419,127606,2294,-691,-1827,-1254,-1,-119,1057,324,505,-195,254,-641,-36,1008,-1082,-3,...
 -87,122,161,11,2,-106,29,-123,-32,41,-524,-35,133,-595,225,837,-108,-191,-2294,841,-340,-394,-351,-1039,238,-108,-66,21,1405,869,520,2776,-174,71,425,652,-1260,-80,249,77, -192,-17,-97,134,-7,-54,-802,-7436,-2824,70869,-35,2481,1865,1749,-2166,2415,33,-183,-835,283,27,-45,56,235,2,718,-1206,275,-87,-158,-7,-2534,0,10774,1,-324,-208,821,281,1340,-797,440,...
 224,72,-65,-5,-7,-44,-48,66,-151,-40,-41,-45,76,-108,-18,1202,0,-2501,1438,-595,900,3040, -3435,-5,-100,-26,0,-13714,-183,68,453,-83,-228,325,97,13,2,105,-61,257,0,57,88,-11,-1,-8220,0,275,-43,-10,-199,105,1,-5849,2,24887,-128,48,712,970,-1407,845,-266,378,311,1526,-1751,27,...
 0,-185858,133,6383,-108,25,-7,1944,5,390,-11,31,277,-384,158,72,-81,-41,-13,-111,-2332,-65804,-698,505812,34,1676716,-72,-6664384,154,-57,52,95,...
 -4,-5,-7,37,-63,-32,4,3349,1,-14370, 16,-83,0,-401,13,3013,48,-20,0,250,51,-79,-7,-146,148,9,0,-64,-17,-59,-67,-492,-2,2116601,-12,-1848,8,-436,-6,324,0,-1363,-163,9,0,-74,63,8167,-29,37587,-22,-74501,-71,497,-1,551747,-87,-22,0,-51,-1,-463,0,-444,3,89,15,-84,-36,-6829,-5,-21663,0,86058,0,-298,-2,751,-2,-1015,0,69,1,-4989,0,21458,0,-330,0,-7,0,-226,0,-1407,0,2942,...
 0,66,0,667,0,-155,0,105,0,-107,0,-74,0,-52,0,91,0,59,0,235,-1,-1819,0,2470,71,13,0,1026,14,-54, 0,-174,-121,-19,0,-200,0,3008,-16,-8043,-10,-37136,-3,73724,-157,-5,0,-854,8,147,-13,-893,0,11869,-23,-172,89,14,-1,872,0,-3744,11,1606,0,-559,-1,-2530,0,454,0,-193,-60,-10,-82,-13,-75,6,36,81,354,-162836,148,-516569,...
 4,2054441,4,-94,39,38,61,-30,2,121,-11,590,62,2108,0,-12242,-476,-42,-84,113,-394,236,0,276,-49,31,0,86,1,-1313,1,69,-60,88,-46,18,0,-63818,14,-93,113,547,-618,17,-7,12290,-1,-25679,...
 0,92,-115,50,-48,233,4,1311,1,-5567,3,1251,29,548,-244,257,-2,1825,42,637,-46,68,-62,8,3,110,445,-100,-316,-202,2925,-621,763,1495,-169,-184,20,-76,-475,-138,8,-141,-197,1351,-1284,422,-129,1879,-102,8382,-9,45864958,-215,1350,-1285,422,-481,-136,8,-140,40,-53,2622,-543,700,1406,402,-95,-318,-194,...
 122,13,-30,147,-121,-902,61,-23,-63,7,69,479,-224,228,-7,500,0,-429,-42,193,-92,37,67,5,-350,-31,0,67,-55,-5,0,47,-36,53,5,561,0,-126,0,871,-52,4,-201,116922,-22,371352,-12,-1473285,0,87,-164,84,-3,422,30,1434,-26,38,2,-1249943,-404,-34,-57,79,5,509,1,131,-344,168,112,22540,30,71218,...
 18,-283983,0,-851,0,-1538,0,1360,-12,51,-48,68,88,-20,1,63,0,-568,303,25,0,-122,87,586,-606,-14,0,-100,-85,8,-165,54,-45,140,0,-54,4,-831,1,3495,...
 31,116,-46,-11,-371,190,-507,399,-2,57,-60,36,-198,-1174,-613,4988,-87,-4,141,560,-276,187,1876,1379,778,4386,24,-15,167,-774,-71,-9,-62,90,98,580,-663,-7,34,-112,57,15,-355,-214,-3240,-13605,12229,-5723,3496,7063,33,-51,1908,1160,-226,715,964,1170,-1264,623,14071,5280,5614,3026,...
 488,1576,-2,226395859,824,1106,-1287,617,1917,1156,-214,718,90,-97,12078,-2366,3282,6668,-219,9179,593,2015,-282,-186,57,25,31,-102,-77,-4,-268,-341,-7,-45,...
 -3,74,15,-615,-88,-7,234,-353,1,-119,-163,-1159,-601,4969,22,-58,-17,-11434,17,54,348,348,-460,434,-371,175,-11,-204,4,-6440,-5,-53,-4,-14388,-37,-45231,-7,179562,-44,136,-160,49,-101,81,-1,-188,0,2,-4,12124,-11,-25217,71,543,-557,-14,...
 -75,526,0,395274,-233,-16,93,-20,-43,61,0,-1275,0,-824,1,-415,0,1762,-261,131,-45,64,-297,-25,0,-17533,-6,-56,21,1100,1,327,1,66,23,-217,-83,-7,83,86847,49,275754,...
 -4,-1093857,-46,2,0,-24,0,-419,0,-5833,1,506,0,-827,-1,-377,-11,-78,0,131945,-2,-334,1,-75,0,-72,0,-213,-6,5564,-2,-11618,0,1790,0,-131,0,6,0,-76,0,-130,0,-1115,0,4783,...
 0,-195,0,-627,0,-55,0,-83,0,163,0,-54,0,82,0,149,0,-754,0,1578,0,138,0,68,2,-2506,0,3399,0,-125,86,16,0,-6350,0,27316,18,-63,0,-169,-1,46,-136,-21,0,-239,-30,-8788,-15,-40549,-4,80514,-46,-8,-168,-6,-1,536,0,-2314,...
 9,148,-13,-842,-1,307713,-23,-175,95,15,0,-297,11,1341,0,-106,0,5,-4,68,-114,10,32,75,159,-130487,98,-413967,2,1647339,-4,-85,100,-46,2,95,-11,461,51,1647,0,-32090,-375,-33,-65,86,-300,180,0,836,0,-3576,...
 0,-222,0,-993,-41,60,0,-4537,-431,-34,2,927,0,-1931,-79,33,-31,144,-1,284,0,-1207, 0,88,-11,315,-178,177,-1,144, -58,986,11,86,-228,-110,2636,-494,718,1474,28,-35,-24,782,-797,277,2142,-1231,856,1853,74,10797,0,23699298,-21,786,-796,277,27,-34, 2615,-494,712,1461,-226,-109,-11,663,0,-123,-169,157,-54,266,0,-76,1,-634,0,2738,-25,106,-63,24,0,-372,-221,-24,0,-5356,0,-219,0,91,-28,7684,-6,24391,...
 -1,-96795,-77,43,2,95,-47,-3,0,-84530,2,310,1,88,111,19331,32,61306,4,-243595,0,770,0,-103,0,160,0,356,0,236,-41,354,39,303,12,-56,873,-143,238,482, -28,35,-93,31,-3,7690211,-91,33,-34,43,824,-130,226,450,-39,341,-1,-687,0,-303,11,-2935,1,12618,121,924,9,-1836,-268,-1144,-678,3685,...
 -69,-261,0,-4115951,-69,-261,5,-151,0,-88,0,91,0,187,0,-1281,1,77,1,6059,3,19238, 0,-76305,0,-90,0,-238,0,-962,0,4133,0,96,0,9483,0,85,0,-688,0,-5607,0,55,0,-752,0,71,0,303,0,-288,0,57, 0,45,0,189,0,401,0,-1474,0,3087,...
 0,-71,0,2925,0,-75,0,359,0,55,1,-10155,0,43735,0,-572,0,-49,0,-660,0,-3591,0,7516,0,668,-1,-53,-2,384259,0,-163,0,-93,1,112,-95,-11528,-22,-36505, -1,145308,5,145,0,4047,1,1483,0,-6352,0,991,0,-4262,0,-93,0,-334,0,-160,0,-153,-10,127,51,185,-77,18,56,1217,6,1919574,-74,17,...
 50,180,-5,93,0,-104,0,-58,-3,-353,-1,1499,0,-229,-15,86,0,-93657,0,1561,0,-6693,0,-5839,1,6791,0,-29143,1,-701,0,3015, 0,2543,0,693,-1,361233,0,-50,0,946,-1,-140,-70,407,0,-450995,0,-368,0,54,0,-802,0,-96,0,1274,0,-5459,0,-614,0,2633,0,685,0,-915,0,-85,0,88,0,106,0,928,0,-726,0,1523,0,5715,0,-4338,0,18706,...
 0,-135,0,-132,0,-158,0,-98,0,680,-1,138968,0,-192,0,-1698,0,-2734,0,11769,0,4,0,673,0,-2891,0,889,0,-3821, 0,121,-1,143783,0,231,-9,51,0,-57413,0,-483,0,-407,0,676,0,-2902,0,531,0,445,0,672,0,19336,0,70,0,-39976,0,-68,0,4203,0,-406,0,446,0,-108,0,79,0,84,0,734,0,255,0,3944,0,-655,0,2825,...
 0,-109,0,-234,0,57,0,19773,0,-2013,0,958,0,-521,0,-757,0,10594,0,-9901,0,199,0,-275,0,64,0,54,0,165,0,1110,0,-3286,0,909,0,54,0,87,0,258,0,1261,0,-51,0,336,0,-114,0,2185,0,-850,0,75,0,-69,0,-103,0,776,0,-1238,0,137,0,67,0,-260,0,130,0,49,0,228,0,215,0,-178,0,57,0,-133];

tabb=[-1];

tabr=[ -5422,-2120,1077,772,39,75,3,10,-468,-326,-113,-78,-4,-2,1,3,29,24,4,2,1,0,-9,7,-2,0,-32,-13,-3,-3,233,126,89,77,-33,16,3,-3,0,-1,...
 2,0,0,1,4,9,1,1,16,-1,0,18, 3,2,0,0,0,0,0,0,0,0,0,-1,-22,-5,10,3,1,1,-15,7,-2,1,-8,-11,-1,-2,-1,1,46,-58,126,-23,4,8,35,8,10,-17,0,0,0,0,-10,-7,0,0,-23,3,151,10,-327,0,4,-5,6,5,...
 1,0,-1,-3,0,0,0,1,-185,0,-3,-24,-5,-2,-1062,3,4560,0,-3,0,4,1,8,-1,2,4,0,1,0,-1,0,0,-1,0,0,1,0,0,-1,-1,277,3,-583,1,-1,4,-32,7,...
 0,-34,1,-1,-23685,0,-1,-2,-1,-7,-5,-4,0,2,-2,0,-5,-1,35,0,0,2,202,0,180,0,0,-1,-3,-6,-193,0,770,-1,-2,-4,-32,23,-28,-46,-13,-9,-54,10,-1,-61,...
 -44895,0,-230,5,-1,-4,-71,0,-15,0,1,0,15,11,-3,6,2,-3,4,-1,2576,-138,-19881,-47,-65906,-1,261925,-4,-2,-7,4,-2,0,0,-1,0,1,-3,...