mar404.m

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

function an=mar404

%First date in file = 1228000.50
%Number of records = 397276.0
%Days per record = 4.0
%      Julian Years      Lon    Lat    Rad
% -1349.9 to  -1000.0:   0.42   0.18   0.25 
% -1000.0 to   -500.0:   0.45   0.14   0.21 
%  -500.0 to      0.0:   0.37   0.10   0.20 
%     0.0 to    500.0:   0.33   0.09   0.22 
%   500.0 to   1000.0:   0.48   0.07   0.22 
%  1000.0 to   1500.0:   0.40   0.07   0.19 
%  1500.0 to   2000.0:   0.36   0.11   0.19 
%  2000.0 to   2500.0:   0.38   0.14   0.20 
%  2500.0 to   3000.0:   0.45   0.15   0.24 
%  3000.0 to   3000.8:  0.182  0.125  0.087 


tabl=[43471.66140,21291.11063,2033.37848,6890507597.78366,1279543.73631,317.74183,730.69258,-15.26502,277.56960,-62.96711,20.96285,1.01857,-2.19395,3.75708,3.65854,0.01049,1.09183,-0.00605,-0.04769,0.41839,0.10091,0.03887,0.11666,-0.03301,0.02664,0.38777,-0.56974,0.02974,-0.15041,0.02179,-0.00808,0.08594,0.09773,...
 -0.00902,-0.04597,0.00762,-0.03858,-0.00139,0.01562,0.02019,0.01878,-0.01244,0.00795,0.00815,0.03501,-0.00335,-0.02970,-0.00518,-0.01763,0.17257,0.14698,-0.14417,0.26028,...
  0.00062,-0.00180,13.35262,39.38771,-15.49558,22.00150,-7.71321,-4.20035,0.62074,-1.42376,0.07043,-0.06670,0.16960,-0.06859,0.07787,0.01845,-0.01608,-0.00914,5.60438,-3.44436,5.88876,6.77238,-5.29704,3.48944,0.01291,0.01280,-0.53532,0.86584,0.79604,0.31635,-3.92977,-0.94829,-0.74254,-1.37947,0.17871,-0.12477,...
  0.00171,0.11537,0.02281,-0.03922,-0.00165,0.02965,1.59773,1.24565,-0.35802,1.37272,-0.44811,-0.08611,3.04184,-3.39729,8.86270,6.65967,-9.10580,10.66103,0.02015,-0.00902,-0.01166,-0.23957,-0.12128,-0.04640,-0.07114,0.14053,-0.04966,-0.01665,0.28411,-0.37754,-1.26265,1.01377,...
 3.70433,-0.21025,-0.00972,0.00350,0.00997,0.00450,-2.15305,3.18147,-1.81957,-0.02321,-0.02560,-0.35188,0.00003,-0.01110,0.00244,-0.05083,-0.00216,-0.02026,0.05179,0.04188,5.92031,-1.61316,3.72001,6.98783,...
 -4.17690,2.61250,0.04157,2.76453,-1.34043,0.74586,-0.20258,-0.30467,0.00733,0.00376,1.72800,0.76593,1.26577,-2.02682,-1.14637,-0.91894,-0.00002,0.00036,2.54213,0.89533,-0.04166,2.36838,-0.97069,0.05486,0.46927,0.04500,0.23388,0.35005,1.61402,2.30209,-0.99859,1.63349,...
 -0.51490,-0.26112,0.27848,-0.26100,-0.07645,-0.22001,0.92901,1.12627,-0.39829,0.77120,-0.23716,-0.11245,-0.02387,0.03960,-0.00802,0.02179,2.86448,1.00246,-0.14647,2.80278,-1.14143,0.05177,1.68671,-1.23451,3.16285,0.70070,0.25817,3.17416,0.07447,-0.08116,-0.03029,-0.02795,0.00816,0.01023,0.00685,-0.01075,-0.34268,0.03680,-0.05488,-0.07430,...
 -0.00041,-0.02968,3.13228,-0.83209,1.95765,3.78394,-2.26196,1.38520,-0.00401,-0.01397,1.01604,-0.99485,0.62465,0.22431, -0.05076,0.12025,4.35229,-5.04483,14.87533,9.00826,-10.37595,19.26596,0.40352,0.19895,0.09463,-0.10774,-0.17809,-0.08979,-0.00796,-0.04313,0.01520,-0.03538,1.53301,-1.75553,4.87236,3.23662,...
 -3.62305,6.42351,-0.00439,-0.01305,0.17194,-0.64003,0.26609,0.06600,0.01767,-0.00251,-0.08871,-0.15523,0.01201,-0.03408,-0.29126,-0.07093,-0.00998,-0.07876,...
  1.05932,-25.38650,-0.29354,0.04179,-0.01726,0.07473,-0.07607,-0.08859,0.00842,-0.02359,0.47858,-0.39809,1.25061,0.87017,-0.82453,1.56864,-0.00463,0.02385,-0.29070,8.56535,-0.12495,0.06580,-0.03395,-0.02465,-1.06759,0.47004,-0.40281,-0.23957,0.03572,-0.07012,0.00571,-0.00731,0.18601,-1.34068,0.03798,-0.00532,0.00448,-0.01147,1.41208,-0.00668,0.25883,1.23788,...
 -0.57774,0.09166,-2.49664,-0.25235,-0.53582,-0.80126,  0.10827,-0.08861,-0.03577,0.06825,-0.00143,0.04633,0.01586,-0.01056,-0.02106,0.03804,-0.00088,-0.03458,-0.00033,-0.01079,0.05821,-0.02445,0.00602,0.00721,-0.00315,-0.01021,-0.65454,1.08478,-0.44593,-0.21492,-1.35004,4.47299,-4.19170,3.51236,...
  1946.04629,13960.88247,576.24572,8023.81797,2402.48512,-753.87007,-6376.99217,-10278.88014,-25743.89874,15506.87748,15609.59853,35173.63133,-3.70370,6.29538,-4.84183,-0.76942,...
 -0.02465,-0.03840,0.00565,-0.06071,0.01174,0.00253,-0.00230,0.05252,-0.02813,0.01359,0.23208,0.03393,0.01734,0.04838,-0.46340,-0.18941,0.25428,-0.56925,0.05213,0.24704,0.12922,-0.01531,0.06885,-0.08510,0.01853,-0.00390,0.01196,-0.30530,0.13117,-0.03533,1.79597,-0.42743,0.98545,2.13503,...
 -1.32942,0.68005,-0.01226,0.00571,0.31081,0.34932,0.34531,-0.32947,-0.00548,0.00186,-0.00157,-0.00065,0.30877,-0.03864,0.04921,0.06693,0.01761,-0.04119,1.28318,0.38546,0.06462,1.18337,...
 -0.48698,0.07086,0.26031,-0.22813,0.10272,0.04737,-0.04506,-0.38581,-0.16624,-0.04588,0.00992,0.00722,-0.21041,0.20560,-0.09267,-0.03438, 0.32264,-0.07383,0.09553,-0.38730,0.17109,-0.01342,-0.02336,-0.01286,0.00230,0.04626,0.01176,0.01868,-0.15411,-0.32799,0.22083,-0.14077,...
  1.98392,1.68058,-0.02526,-0.13164,-0.04447,-0.00153, 0.01277,0.00553,-0.26035,-0.11362,0.14672,-0.32242,0.16686,-0.69957,0.40091,-0.06721, 0.00837,0.09635,-0.08545,0.25178,-0.22486,16.03256,0.34130,-0.06313,0.01469,-0.09012,...
 -0.00744,-0.02510,-0.08492,-0.13733,-0.07620,-0.15329,0.13716,-0.03769,2.01176,-1.35991,-1.04319,-2.97226,-0.01433,0.61219,-0.55522,0.38579,0.31831,0.81843,...
 -0.04583,-0.14585,-0.10218,0.16039,-0.06552,-0.01802,0.06480,-0.06641,0.01672,-0.00287,0.00308,0.09982,-0.05679,-0.00249,-0.36034,0.52385,-0.29759,0.59539,-3.59641,-1.02499,-547.53774,734.11470,441.86760,-626.68255,-2255.81376,-1309.01028,-2025.69590,2774.69901,...
  1711.21478,1509.99797,-0.99274,0.61858,-0.47634,-0.33034,  0.00261,0.01183,-0.00038,0.11687,0.00994,-0.01122,0.03482,-0.01942,-0.11557,0.38237,-0.17826,0.00830,  0.01193,-0.05469,0.01557,0.01747,0.02730,-0.01182,-0.11284,0.12939,-0.05621,-0.01615,...
  0.04258,0.01058,-0.01723,0.00963,0.20666,0.11742,0.07830,-0.02922,-0.10659,-0.05407,0.07254,-0.13005, -0.02365,0.24583,0.31915,1.27060,0.00009,-0.21541,-0.55324,-0.45999,-1.45885,0.86530, 0.85932,1.92999,-0.00755,-0.00715,-0.02004,-0.00788,0.01539,0.00837,0.27652,-0.50297,-0.26703,-0.28159,...
  0.03950,0.07182,-0.07177,0.14140,0.07693,0.07564, -0.01316,-0.01259,0.01529,0.07773,-90.74225,-378.15784,-510.30190,-52.35396,-89.15267,415.56828,181.52119,54.01570, -0.01093,-0.05931,-0.01344,-0.02390,0.01432,-0.02470,-0.01509,-0.01346,0.03352,0.02248,0.02588,-0.00948,0.03610,0.17238,0.02909,-0.04065,0.00155,-0.07025,-0.09508,0.14487,0.12441,0.16451,...
  0.00001,-0.00005,-0.00982,-0.01895,-0.16968,0.36565,0.20234,0.17789, -0.04519,-0.00588,0.01268,0.00107,-56.32137,-58.22145,-80.55270,28.14532, 11.43301,52.05752,17.79480,-2.61997,-0.00005,-0.02629,0.01080,-0.00390,0.00744,0.03132,0.01156,-0.01621,0.02162,0.02552,0.00075,-0.02497,0.02495,0.00830,0.03230,0.00103,-14.84965,-4.50200,-9.73043,9.40426,4.08054,5.38571,1.53731,-1.01288,0.21076,1.74227,0.79760,0.39583,0.09879,-0.16736,-0.00723,-0.01536];


tabb=[ -364.49380,-47.17612,-554.97858,-430.63121,596.44312,-3.94434,-7.43169,-0.06665,-2.23987,0.10366,-0.05567,-0.01463,0.01908,-0.02611,-0.00350,-0.01057,-0.00610,...
 -0.00015,0.00002,0.00010,0.00033,0.00007,-0.00000,-0.00010,-0.00004,0.00012,0.00002,-0.00014,-0.00048,-0.00003,-0.00007,...
 0.00008,-0.00005,-0.00043,-0.00003,-0.00010,-0.00004,0.00001,0.00001,-0.00003,-0.00003,0.00004,0.00007,-0.00041,0.00031,0.00076,0.00062,0.00001,-0.00002,0.00035,0.00053,0.00026,0.00019,...
 0.00020,0.00010,0.02936,0.09624,-0.01153,0.01386,0.00551,-0.00690,0.00196,0.00148,-0.00408,-0.00673,-0.00067,-0.00152,...
 -0.00014,-0.00005,0.00000,0.00005,-0.00116,0.00276,-0.00391,0.00983, -0.01327,-0.01986,-0.00003,0.00001,0.01104,0.00631,-0.01364,0.01152,...
 -0.00439,0.01103,-0.00546,0.00181,-0.00039,-0.00083,0.00007,0.00002,-0.00010,-0.00008,0.00005,0.00002,-0.00584,0.00512,-0.00722,-0.00174,...
 0.00101,-0.00316,-0.02229,-0.02797,-0.10718,0.05741,0.11403,0.10033,0.00036,-0.00022,0.00787,0.01191,0.01756,-0.02121,-0.00169,-0.00364,0.00070,-0.00051,0.01850,-0.06836,0.21471,0.00162,...
 -0.29165,0.16799,-0.00002,0.00011,-0.00075,-0.00077,-0.00675,-0.00814,0.00029,-0.00599,0.00107,0.00013,0.00010,-0.00002,0.00005,0.00020,0.00355,0.00306,-0.00013,-0.00061,-0.02950,-0.00847,0.01037,-0.04783,...
 0.04237,0.11662,-0.00331,0.00207,-0.00107,-0.00264,0.00072,-0.00023,-0.00151,0.00146,-0.12847,0.02294,0.03611,0.19705,...
 0.16855,-0.28279,-0.00000,-0.00002,-0.00525,-0.03619,0.05048,-0.00481,-0.00745,0.04618,0.00286,0.00443,0.00521,-0.00351,...
 0.00200,0.00474,-0.00149,0.00031,-0.00003,0.00029,0.00686,0.02467,0.04275,-0.02223,0.02282,-0.04228,0.03312,0.01847,-0.01253,0.01601,0.00076,0.00091,0.00045,0.00035,0.00658,0.01586,-0.00310,0.00628,...
 -0.00045,0.00316,-0.01602,-0.00340,-0.01744,0.04907,0.06426,0.02275,-0.00217,-0.00377,-0.00091,0.00037,0.00040,-0.00003,-0.00017,-0.00027,0.00366,0.02693,-0.00934,0.00386,0.00616,-0.00037,0.02028,0.02120,-0.01768,0.02421,0.00102,0.00877,0.00012,0.00030,-0.00019,-0.02165,0.01245,-0.00742,...
 0.00172,0.00320,-0.17117,-0.12908,-0.43134,0.15617,0.21216,0.56432,0.01139,-0.00937,-0.00058,-0.00337,-0.00999,0.01862,-0.00621,-0.00080,-0.00025,-0.00140,0.09250,0.01173,-0.03549,0.14651,-0.01784,0.00945,0.00000,-0.00006,-0.00500,0.00086,0.01079,-0.00002,...
 -0.00012,-0.00029,-0.02661,0.00140,-0.00524,-0.00460,-0.00352,-0.00563,-0.00277,-0.00052,-0.10171,-0.02001,0.00045,0.00265,-0.00082,0.00160,-0.00302,-0.00434,-0.00022,-0.00134,0.03285,0.02964,-0.05612,-0.00668,...
 -0.01821,0.06590,0.00039,0.00061,-0.13531,-0.03831,0.02553,0.02130,-0.00336,0.00468,-0.04522,-0.05540,0.00129,-0.01767,0.00181,0.00031,-0.00011,-0.00034,-0.00146,0.01101,-0.00030,0.00240,-0.00039,0.00072,-0.01954,-0.03822,0.09682,-0.04541,-0.01567,0.09617,-0.03371,0.33028,-0.12102,0.05874,...