plu404.m

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

 2.80629,-1.12755,-0.00402,0.00053,0.00024,-0.00043,0.00403,-0.00210,0.00603,0.00411,-0.00260,0.00416,2.29235,3.05992,2.36465,-0.58750,0.14030,0.13523,0.89998,0.70156,-0.02188,0.02003,-0.00533,0.00447,...
 2.96411,1.30183,0.01422,0.00624,-0.10737,-0.38316,-0.05968,0.04379,0.01171,0.01180,-0.00989,-0.01375,-0.00845,0.03782,0.09484,0.09909,0.07640,-0.00898,-0.01076,0.02760,0.01630,0.02198,0.05985,0.04130];
 
 tabr=[ 17990649.12487,24806479.30874,12690953.00645,2892671.69562,249947.71316,-5138.71425,1142.68629,6075.25751,-34.76785,-19.72399,-15.81516,30.47718,-11.73638,21.87955,9.42107,-10.40957,-5.59670,6.85778,-167.06735,-2.31999,-32.42575,-13.72714,130.16635,117.97555,31.33915,39.64331,-1378.54934,-395.83244,-562.79856,-167.74359,-45.12476,-17.08986,-4.20576,-16.56724,1762.12089,-1148.86987,736.55320,-423.09108,...
 56.13621,-26.26674,9.77810,-38.05151,4702224.98754,27254904.94363,5306232.25993,39518429.29982,1725110.05669,21833263.27069,46010.62605,5425411.66252,17238.09865,536771.62156,-61263.36051,66270.70142,2084.66296,-1936.71208,35898.49503,34885.28549,1.93276,10.66292,-665.11445,3.70467,-265.68478,16.16272,-19.45954,2.32738,3.04237,3.97339,-2.64312,0.66308,...
 -3207.68754,3418.03720,-2342.62310,1729.15030,-450.84643,179.00943,-13.20367,-1.86087,-4.95659,7.22347,...
 -5.08890,-1.28891,-6.21713,5.10551,13.97276,0.44529,3.25177,25.02775,-45.56672,11.58470,124443355.55450,-100018293.41775,190506421.77863,-118262753.40162,108199328.45091,-45247957.63323,27272084.41143,-4125106.01144,2583469.66051,1024678.12935,-22702.55109,199269.51481,-15783.14789,5564.52481,-427.22231,-6330.86079,-97.50757,-204.32241,-9060.54822,156661.77631,-47791.83678,59725.58975,-8807.74881,-92.38886,-28886.11572,-244419.59744,-53336.36915,-92232.16479,...
 -8724.89354,-2446.76739,889.71335,936.51108,494.80305,2252.83602,-18326.60823,-25443.13554,-3130.86382,-5426.29135,23494.08846,91.28882,4664.14726,1552.06143,-8090.43357,2843.48366,-1445.73506,1023.11482,11664.20863,-7020.08612,3100.21504,-64.16577,-9724.97938,-12261.47155,-3008.08276,-1523.06301,6788.74046,10708.27853,343.09434,1701.52760,...
 14743.99857,-4781.96586,-15922236.41469,1825172.51825,-14006084.36972,10363332.64447,-979550.91360,6542446.18797,1160614.26915,570804.88172,89912.68112,-171247.08757,-13899.52899,-6182.25841,-240.64725,412.42581,-66.24510,71.30726,-15.81125,-15.76899,-21.85515,-102.12717,-10.18287,-19.38527,...
 1.43749,-3.87533,1.97109,0.20138,0.32012,1.02928,-40.22077,20.80684,-15.69766,9.63663,-1.26010,0.56197,0.08592,0.18540,-0.07303,0.03897,0.01438,-0.08809,0.15479,0.10354,0.19052,2.08790,...
 405480.24475,-607986.83623,582811.58843,-915111.10396,258696.21023,-493391.09443,23403.62628,-119503.67282,-4036.86957,-9766.17805,-663.93268,2544.07799,40.36638,76.26390,246.67716,-13.93440,...
 0.12403,0.25378,0.14004,-0.08501,0.07904,0.12731,1.02117,-1.34663,0.25142,-0.26903,0.18135,-0.57683,-0.30092,-0.36121,...
 -0.09623,0.05873,-0.05803,0.02869,-0.01194,0.04983,0.04250,0.04894,1.34245,0.70137,0.24217,0.25912,-0.32759,-0.03575,0.06780,-0.41277,0.43865,0.17857,-763933.02226,465658.17048,-1082753.91241,593319.68634,...
 -553911.89340,274748.95145,-122250.71547,56608.95768,-9914.17300,2988.43709,707.94605,-765.01470,52.73260,-34.22263,-43.58300,-38.43647,-4.95939,-1.97173,-1.04406,-0.13072,-0.34281,4.75202,-0.35513,0.93597,...
 -0.54380,0.70536,84.83116,102.93003,26.34884,48.57746,0.02853,2.91676,-8.07116,1.66613,-2.07908,11.62592,6.64704,0.98291,-1.19192,0.93791,0.18822,0.00900,-0.03181,-0.02000,...
 0.02755,-0.01398,-0.03971,-0.03756,0.13454,-0.04193,-18672.98484,28230.75834,-28371.58823,26448.45214,-13352.09393,7461.71279,-2609.33578,726.50321,-309.72942,-86.71982,12.48589,-9.69726,...
 1.82185,14.92220,-0.04748,0.42510,-0.20047,0.00154,0.00176,-0.26262,0.78218,-0.73243,0.23694,-0.03132,-0.00290,-0.03678,14.03094,4.25948,0.79368,-0.78489,...
 -2.30962,2.31946,0.00158,-0.04125,-0.01387,0.28503,0.00892,0.05154,0.00184,-0.01727,-0.00889,0.03526,-0.00521,-0.02093,0.00200,0.04872,-0.02163,0.00578,20699.27413,-2175.57827,31177.33085,4572.02063,...
 15486.28190,8747.74091,2455.51737,3839.83609,51.31433,507.91086,15.90082,44.75942,-0.98374,-2.64477,2.52336,-3.09203,-0.08897,-0.00083,-15.91892,0.72597,14.04523,-3.16525,4.33379,-30.82980,0.40462,-0.75845,13.14831,-0.02721,-0.01779,0.00481,0.42365,-0.09048,0.08653,0.04391,...
 0.00846,0.01082,-0.04736,0.02308,6282.21778,-4952.70286,7886.57505,-5328.36122,3113.76826,-1696.84590,330.70011,-155.51989,-18.31559,-3.90798,-3.11242,1.87818,-1.05578,0.11198,0.05077,-0.01571,2.41291,2.40568,-0.01136,-0.00076,-0.00392,-0.02774,...
 634.85065,-352.21937,674.31665,-260.73473,199.16422,-28.44198,6.54187,6.44960,-1.55155,0.29755,0.16977,0.17540,-0.02652,0.03726,-0.00623,0.11777,-0.00933,0.02602,-0.13943,-0.24818,0.02876,-0.01463,-0.07166,0.06747,-0.01578,0.01628,0.00233,-0.00686,0.00431,-0.00276,0.21774,0.09735,0.07894,0.07279,...
 -0.01300,-0.00268,0.10824,0.09435,0.00720,0.02111,...
 -0.01960,0.06154,0.56867,-0.07544,0.18210,0.06343,-0.00906,0.01942,-0.00850,-0.00351,-0.06988,0.01713,-0.01110,-0.00663,0.00196,-0.02064,-0.00008,0.00043,...
 0.00375,0.00084,-0.00279,0.00100,0.00271,-0.02017,-0.00074,-0.00357,0.03793,-0.10108,-0.01083,-0.03952,...
 0.00030,0.00012,0.01576,0.01142,0.00351,0.00277,0.01409,-0.00774,-0.00065,0.01895,0.07350,-0.02519,0.01528,-0.01057,...
 -0.00099,-0.00295,0.21347,-0.17458,0.04940,-0.02757,-0.06243,0.05203,0.01055,-0.00109,0.00003,-0.04201,-0.00263,0.02387,0.00886,-0.01168,0.00479,0.00204,-0.00239,0.00022,-0.00223,-0.02029,-0.14130,-0.15237,-0.01827,-0.04877,...
 0.12104,0.06796,0.16379,0.31892,-0.15605,0.07048,-0.00700,0.07481,-0.00370,-0.00142,-0.00446,0.00329,-0.00018,0.00117,-0.00910,0.00510,-0.00055,-0.00114,0.04131,-0.04013,-0.13238,0.02680,-0.10369,1.38709,0.35515,0.41437,-0.01327,-0.02692,38.02603,13.38166,15.33389,-7.40145,-8.55293,-0.13185,-0.03316,0.13016,...
 0.04428,-1.60953,-12.87829,-76.97922,-23.96039,-22.45636,14.83309,14.09854,0.24252,0.13850,-4.16582,4.08846,0.00751,-0.00051,0.03456,0.02900,0.01625,-0.04660,0.01390,-0.00530,0.01665,-0.04571,40.90768,-14.11641,7.46071,-58.07356,-0.27859,-1.33816,23.76074,-0.03124,-0.27860,0.13654,-0.04800,0.05375,...
 4.38091,4.39337,0.02233,0.00514,-0.25616,-0.54439,-0.05155,0.11553,0.02944,-0.00818,0.00570,0.00119,-0.00733,-0.02700,-0.23759,-0.08712,-0.12433,0.07397,0.20629,0.60251,0.56512,0.14790,0.07778,0.11614];

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

%Total terms = 173, small = 156 
plu.maxargs=9;
plu.max_harmonic=[   0,  0,  0,  0,  2,  2,  9, 13, 13];
plu.max_power_of_t=7;
plu.arg_tbl=args;
plu.lon_tbl=tabl;
plu.lat_tbl=tabb;
plu.rad_tbl=tabr;
plu.distance=   3.9539999999999999e+01;
plu.timescale=3652500.0;
plu.trunclvl=1.0;

an=plu;