ura404.m

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

 0.00083,0.00020,0.00140,-0.00176,0.00381,-0.00731,-0.01618,0.01570,-0.10201,0.05809,-0.03359,0.01024,-0.00535,0.00018,0.00024,0.00509,-0.00158,...
    -0.00466,0.00009,-0.00083,-0.00700,-0.00090,-0.00011,-0.00079,0.00133,-0.00126,0.01416,0.05553,0.04283,-0.06719,0.00119,0.00291,-0.00263,0.01282,-0.00040,0.00188,-0.00237,0.00973,-0.39533,0.18773,-0.79821,-0.40168,0.00151,-0.00161,0.00123,-0.00516,-0.01432,-0.00293,-0.05477,0.04130,-0.48837,0.18944,-0.12552,9.37098,1.02045,5.11382,0.72098,-3.70049,-5.80982,3.30105,-0.09682,0.09696,-0.00876,0.00504,0.00318,0.00245,0.00563,-0.00665,0.00108,-0.00233,-0.00117,0.00177,-0.00343,0.00503,0.01044,-0.00651,0.00296,-0.00162,0.00037,0.00028,-0.00020,-0.00786,0.00029,0.00836,0.00004,0.00033,-0.00309,-0.00086,-0.00157,-0.00086,...
 -0.00058,0.00105,-0.04557,0.01794,-0.00122,-0.00086,0.00420,-0.00285,0.00118,-0.00020,0.00743,-0.01217,0.00053,-0.00084,...
 -0.00075,0.00097,-0.00107,0.00314,0.00576,-0.00505,0.03624,-0.02546,0.05379,0.30081,0.29870,-0.22106,0.00696,-0.00801,-0.03995,-0.01808,-0.00139,0.00102,-0.00059,0.00138,0.00019,-0.00037,0.00274,0.00658,0.00672,-0.01132,...
 0.00023,0.00051,0.00031,0.00090,-0.00017,-0.00001,0.00085,0.00004,0.02221,-0.01977,0.07498,0.03025,-0.00082,-0.00022,-0.00073,-0.00028,-0.00253,0.00259,-0.01329,0.01805,0.00096,0.00833,-0.11836,0.04277,-0.10820,-0.03018,...
 0.34504,0.09834,-0.00538,-0.00231,0.00036,0.00042,-0.00023,0.00260,-0.01137,0.00036,0.01081,-0.03271,-0.00029,-0.00028,0.00018,-0.00003,0.00009,0.00012,0.00127,0.00343,0.00100,-0.00064,0.00014,0.00004,0.00150,0.00069,-0.01484,0.00135,0.03930,0.01405,0.00064,0.00029,0.00009,0.00009,0.00054,-0.00048,0.00019,0.00005,-0.00009,0.00018,0.00192,-0.00333,0.01824,0.01071,0.00107,-0.00341,0.25530,-0.18414,-0.84151,-0.31475,...
 -0.00400,-0.00010,-0.00174,0.00019,0.00006,-0.00079,0.00066,-0.00070,0.00599,0.00330,-0.00160,-0.00013,-0.00067,-0.00006,-0.00176,-0.00111,0.00652,0.00368,0.00004,0.00001,-0.00081,0.00089,0.00366,0.00139,0.00002,0.00001,-0.01870,-0.00998,-0.00020,-0.00007,0.00005,0.00003];
 
 tabr=[ 0.00000,-53.23277,-44.70609,-62.54432,-19.15218,0.10867,-1.91911,1.47517,16.51994,5.00458,3.88980,1.55740,3598.17109,1831.07574,2633.34851,1775.69482,497.10486,488.77343,6.03892,31.08365,-2.06585,-1.12599,...
       230.37762,-113.95449,162.40244,-46.57185,6.70207,17.27241,-0.66092,-14.42065,-0.01044,-0.00287,-0.03894,-0.01663,0.01629,0.00496,0.08411,0.02855,0.01795,-0.00695,0.02426,-0.03921,-0.24495,-0.77369,-0.31404,0.38668,...
       -0.05682,-0.17197,0.06145,-0.00510,0.00606,-0.00886,-0.00370,-0.00588,0.02173,-0.11909,0.00302,-0.01796,-0.01067,0.00990,0.05283,0.06517,59710.89716,-491.12783,58672.38609,19564.41947,10597.99050,14313.02561,-2585.52040,766.78396,-138.39893,-802.43403,131.35006,-31.97561,...
       7.95978,8.16075,28.72669,31.72473,6.45792,16.50701,0.01066,1.29718,0.11565,-0.13240,0.05110,-0.01543,-0.09994,0.18864,-0.01330,0.04148,0.03510,-0.00366,0.00604,-0.00604,0.03752,-0.00256,-7.00488,-21.63748,1.43064,-17.10914,-0.62987,0.48719,0.00697,-1.22665,-0.14435,-0.00550,...
       0.32008,-0.19855,-13976.73731,-3559.49432,-7709.90803,-9310.80334,749.31835,-3491.50696,540.94979,-84.57550,16.96663,35.53930,37214.64771,-36361.15845,21093.74492,-31855.33076,1500.84653,-7031.97901,-453.40865,-18.36692,...
 -2.07726,-17.92336,-56348.30507,378512.71483,-111444.43340,370543.95160,-61893.70301,112131.05507,-11977.44617,9156.15245,-567.61838,-495.25760,16.96202,-44.06279,4.24760,-48.83674,-643705.49516,-131013.09649,-838580.02217,67627.11556,-288441.70339,150227.25291,-2500.57537,42676.19888,7084.60505,2043.65642,9639.56835,...
 -1502.03390,-4126.00409,-828.73564,-2801.35204,-2293.77751,-209.23365,-1045.31476,95.57334,-102.74623,7.19216,1.89593,-0.05661,0.02166,120.38332,-141.16507,98.31386,-40.23448,10.84269,17.57713,1.69239,1.45065,...
 -0.19626,2.76108,-0.00270,0.00360,-0.02333,-0.00710,-0.01035,0.02950,0.00737,-0.06311,-0.00613,0.01407,0.01377,0.00879,-0.03287,0.00012,-0.21667,0.01793,-1.54865,0.10953,0.54543,0.12102,-9.48047,0.11477,-1.34966,0.23199,-1.50834,0.26567,-0.64503,0.10742,-0.21452,0.04428,-0.01920,-0.00906,-0.09378,0.12773,-0.02787,-0.03090,...
 0.03111,0.00140,0.03771,-0.01269,-1.94794,1.22823,0.64183,-1.11467,-0.19301,-0.27357,0.05710,-0.08115,-0.07318,0.00806,0.14286,0.20297,0.14920,-0.07897,0.09682,0.02379,-0.13928,0.01679,-0.00774,0.10060,0.24433,0.16760,-2.88905,-1.61439,2.83052,-3.41031,...
 36.37048,3.37867,0.29321,0.09687,0.29324,-0.14651,8.11116,1.79211,1.36421,0.88111,1.21683,2.37950,-357.76211,-87.84636,-117.55745,-67.18338,-5.26029,-6.27559,7509.94562,3.68942,4223.62097,-1041.13557,-74.64464,-251.41613,-166.22180,-1.68190,-214.55340,...
 62.79593,-0.08250,-0.15936,-0.03830,0.10857,0.21368,0.50812,0.00869,0.09832,0.02158,0.02045,0.01407,0.03591,0.03460,0.01171,-0.16400,0.09751,0.03521,-0.12858,0.00700,-0.00524,0.01698,-0.04796,0.04006,0.00565,...
 -0.02783,-0.00205,-0.02296,0.00153,-0.16139,0.01514,-0.78136,-0.01546,0.40374,-0.06014,0.06212,-0.01828,0.00831,-0.00173,0.06857,-0.11677,0.00028,0.05765,-0.00796,0.00691,0.03764,0.14902,-0.02653,0.02122,...
 -0.05503,0.01549,1.56630,-0.35551,-1.87960,1.14303,-0.06063,-0.03425,0.03367,-0.11969,0.04485,-0.01651,0.04647,-0.02097,0.22841,0.47362,0.99226,-0.60660,...
 -0.01249,0.00134,-0.07435,0.00722,-0.31796,-0.00015,0.20533,-0.04398,0.93944,-0.26710,-5.60051,-9.32918,-5.13538,-4.05130,-0.56529,4.34112,7.18308,-2.66103,0.13241,-0.07999,0.01046,-0.00535,-0.04037,-0.00455,-0.00510,0.00731,-0.04576,0.00513,-0.15846,-0.00236,0.04628,-0.00463,-0.01585,0.00585,-0.00213,0.00283,...
 0.00778,-0.00198,-0.17803,0.18321,0.07702,-0.12325,0.01091,0.00349,0.14211,-0.21830,0.07289,-0.00994,0.07090,-0.00079,4.18441,-0.07413,-0.06247,-0.00011,-0.15453,0.14499,-0.06557,-0.00098,0.00290,0.02921,-0.01923,0.00457,-0.07538,-0.00120,0.02263,-0.00037,-0.01061,0.00591,-0.04725,0.02364,-0.07460,-0.24108,-0.28310,0.14643,...
 -0.00700,0.00427,0.22963,0.03713,-0.02062,0.00478,0.01434,0.00095,-0.01425,0.00376,0.29611,-0.08038,-0.37811,0.21703,-0.00723,-0.00924,-0.02736,0.01814,0.00934,0.00731,0.00613,0.00686,-0.91503,...
 -0.32009,-0.15505,0.79589,-0.00555,-0.01536,-0.00698,0.00480,0.00373,-0.00046,0.00715,-0.00470,-0.01970,-0.05238,0.60649,-0.32669,0.17790,0.33383,...
 -2.74922,-0.25827,-0.07862,0.00406,-0.00948,-0.02117,0.03127,-0.04199,0.89670,-0.02413,0.01954,0.03990,0.00063,-0.00071,-0.00226,0.02009,-0.04407,-0.05069,0.38230,0.16101,0.11893,-0.06125,0.02051,-0.00046,0.39211,0.03679,0.01666,-0.31336,53.28735,-0.01791,-0.39414,...
 0.04181,-0.01885,0.00165,0.31349,-0.47359,0.16133,-0.01023,0.00007,0.01758,-0.13351,0.07249,0.00977,0.05445,0.11650,-0.00191,-0.09824,0.40106,2.41155,-0.30655,...
 0.24975,-0.01248,-0.03688,0.01097,0.00038,-0.00051,-0.04736,0.02610,0.00968,0.02634,0.07918,-0.00606,0.02735,-0.00320,-0.07544,-0.00468,0.19996,-0.01964,0.00201,0.00267,0.39562,0.43289,1.24743,0.31084,-0.00666,0.00377,0.05668,0.00148,0.03220,-0.00026,0.03717,0.01509];

args=[ 0,3,2,1,7,-2,8,0,2,2,7,-4,8,0,2,3,7,-6,8,1,2,2,5,-5,6,4,2,1,6,-3,7,3,3,1,6,-1,7,-4,8,0,3,2,5,-7,6,6,7,0,3,2,6,-6,7,1,8,0,3,2,6,-7,7,3,8,0,3,2,6,-8,7,4,8,0,3,2,6,-7,7,2,8,0,2,2,6,-6,7,2,3,1,5,-4,6,4,7,0,3,1,6,-2,7,-1,8,0,3,1,6,-3,7,1,8,0,3,1,6,-4,7,3,8,1,...
 2,5,7,-9,8,0,2,4,7,-7,8,0,2,2,7,-3,8,6,2,1,7,-3,8,2,2,2,7,-5,8,0,2,3,7,-7,8,0,3,1,6,-6,7,5,8,1,3,1,6,-5,7,3,8,0,3,2,5,-8,6,8,7,0,3,1,5,-4,6,5,7,0,2,2,6,-5,7,3,3,1,6,1,7,-9,8,0,3,2,5,-4,6,-2,7,0,...
 2,1,6,-4,8,4,2,1,6,-2,7,4,2,5,7,-8,8,5,2,3,7,-4,8,0,1,1,7,5,2,2,7,-6,8,4,3,1,6,-6,7,4,8,0,2,1,6,-4,7,4,3,2,6,-5,7,1,8,0,3,2,6,-6,7,3,8,0,2,2,6,-7,7,0,3,1,5,-4,6,3,7,0,3,1,6,-1,7,-1,8,0,...
 2,1,5,-2,6,0,2,6,7,-9,8,0,2,5,7,-7,8,0,2,4,7,-5,8,0,2,3,7,-3,8,1,2,2,7,-1,8,0,2,1,7,1,8,2,1,3,8,0,2,3,6,-7,7,1,3,2,5,-3,6,-4,7,0,3,2,6,-3,7,-2,8,0,2,2,6,-4,7,1,3,2,6,-5,7,2,8,1,3,5,5,-9,6,-8,7,0,3,2,5,-4,6,-1,7,0,3,1,6,3,7,-8,8,0,3,2,6,-8,7,1,8,0,3,2,5,-7,6,4,7,0,3,4,5,-10,6,2,7,0,...
 2,1,6,-2,8,0,2,1,6,-1,7,2,2,8,7,-12,8,0,2,7,7,-10,8,0,2,6,7,-8,8,1,2,5,7,-6,8,0,2,4,7,-4,8,2,1,2,7,4,1,4,8,0,2,1,7,-6,8,0,2,2,7,-8,8,1,2,3,7,-10,8,0,2,4,7,-12,8,0,3,1,6,-6,7,2,8,0,2,1,6,-5,7,1,3,1,6,-4,7,-2,8,0,3,1,5,-4,6,2,7,1,3,1,5,-2,6,1,7,0,2,7,7,-9,8,0,2,6,7,-7,8,0,2,5,7,-5,8,0,2,4,7,-3,8,0,2,3,7,-1,8,0,2,2,7,1,8,0,2,3,6,-6,7,1,3,3,6,-7,7,2,8,0,3,2,5,-3,6,-3,7,1,3,2,6,-2,7,-2,8,0,2,2,6,-3,7,1,3,2,6,-4,7,2,8,0,...
 3,2,5,-7,6,3,7,1,3,1,6,1,7,-2,8,0,1,1,6,1,2,8,7,-10,8,0,2,7,7,-8,8,0,2,6,7,-6,8,0,2,5,7,-4,8,0,2,4,7,-2,8,0,1,3,7,3,2,2,7,2,8,0,2,1,7,4,8,0,2,1,5,-4,7,0,2,1,6,-6,7,0,2,8,7,-9,8,0,2,7,7,-7,8,0,2,6,7,-5,8,0,2,5,7,-3,8,0,2,4,7,-1,8,0,3,3,6,-4,7,-2,8,0,2,3,6,-5,7,1,3,3,6,-6,7,2,8,0,3,2,5,-3,6,-2,7,1,3,2,6,-1,7,-2,8,0,2,2,6,-2,7,0,3,2,6,-3,7,2,8,0,3,2,5,-7,6,2,7,1,...
 2,1,6,1,7,0,2,9,7,-10,8,0,2,8,7,-8,8,0,2,7,7,-6,8,0,2,6,7,-4,8,0,2,5,7,-2,8,0,1,4,7,1,2,3,7,2,8,0,2,1,5,-3,7,0,2,9,7,-9,8,0,2,8,7,-7,8,0,3,3,6,-3,7,-2,8,0,...
 2,3,6,-4,7,1,3,3,6,-5,7,2,8,0,3,2,5,-3,6,-1,7,0,3,2,5,-8,6,4,7,0,2,2,6,-2,8,0,2,2,6,-1,7,1,3,2,6,-2,7,2,8,0,3,2,5,-7,6,1,7,0,2,6,7,-2,8,0,1,5,7,0,3,3,6,-4,7,1,8,0,...
 2,1,5,-2,7,2,3,1,5,-3,7,2,8,0,3,1,5,-1,6,1,7,0,2,4,6,-6,7,0,2,3,6,-3,7,0,1,2,6,0,3,2,5,-4,6,3,7,0,3,1,5,1,6,-4,7,0,3,3,5,-5,6,-1,7,0,1,6,7,1,3,1,5,1,7,-4,8,0,...
 2,1,5,-2,8,0,2,1,5,-1,7,1,3,1,5,-2,7,2,8,0,3,1,5,-3,7,4,8,0,3,1,5,-5,6,1,7,1,3,1,5,-1,6,2,7,0,2,4,6,-5,7,0,2,3,6,-2,7,0,3,1,5,1,7,-2,8,0,1,1,5,1,2,4,6,-4,7,0,2,3,6,-1,7,0,3,3,5,-5,6,1,7,0,...
 2,5,6,-6,7,0,2,4,6,-3,7,0,2,5,6,-5,7,0,2,6,6,-6,7,0,2,2,5,-3,7,0,2,2,5,-2,7,0,2,2,5,-2,8,0,2,2,5,-1,7,1,3,2,5,-2,7,2,8,0,1,2,5,0,2,3,5,-3,7,0,2,3,5,-1,7,0,-1];

ura.maxargs=9;
ura.max_harmonic=[  0,  0,  0,  0,  5, 10,  9, 12,  0];
ura.max_power_of_t=6;
ura.arg_tbl=args;
ura.lon_tbl=tabl;
ura.lat_tbl=tabb;
ura.rad_tbl=tabr;
ura.distance=   1.9218446061800002e+01;
ura.timescale=3652500.0;
ura.trunclvl=1.0;

an=ura;

% Total terms = 177, small = 171