plu404.m

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

function an=plu404

%First date in file = 625296.50
%Number of records = 16731.0
%Days per record = 131.0
%      Julian Years      Lon    Lat    Rad
% -3000.0 to  -2499.7:   1.17   0.90   0.83 
% -2499.7 to  -1999.7:   0.57   0.51   0.58 
% -1999.7 to  -1499.7:   0.63   0.39   0.40 
% -1499.7 to   -999.8:   0.40   0.45   0.41 
%  -999.8 to   -499.8:   0.42   0.22   0.30 
%  -499.8 to      0.2:   0.41   0.24   0.35 
%     0.2 to    500.2:   0.58   0.24   0.26 
%   500.2 to   1000.1:   0.47   0.35   0.33 
%  1000.1 to   1500.1:   0.43   0.31   0.28 
%  1500.1 to   2000.1:   0.37   0.40   0.35 
%  2000.1 to   2500.0:   0.46   0.35   0.39 
%  2500.0 to   3000.0:   1.09   0.70   0.46 
%  3000.0 to   3000.4:  0.871  0.395  0.051 

tabl=[ 74986469.33577,100898343.73690,48199471.54076,9520076.03177,690431.67340,-427355.12716,52266623.77862,860190.70714,-21.08511,-143.39295,-126.71124,48.18528,-88.74508,40.50942,51.29367,-10.24075,0.63094,32.03258,-410.00781,399.90234,-116.54319,51.50329,116.84565,687.76781,-13.38605,182.70107,-1668.12226,-5162.22648,-585.68297,-2247.56041,-20.95875,-193.13703,-57.12097,-10.35058,-1778.01442,-6463.73779,-657.86093,-2713.44436,...
 -41.32560,-211.82042,-107.16513,-36.58628,97929588.08231,-33582903.63417,143382679.31770,-47411568.92345,79428827.73581,-24658834.87499,19861233.17488,-5755585.62084,1995060.95931,-693507.08147,135176.31467,109360.38672,-8188.00598,-1680.95072,71020.19608,-70785.39049,-24.56034,-20.34919,1618.45976,-2079.48538,635.62954,-850.87068,44.95257,-64.04459,-18.61475,-1.77734,9.38079,5.19958,17422.08783,-4872.53852,10985.76629,-1137.68569,1771.28893,288.93745,40.22664,37.90027,...
 6.81361,-32.65868,16.97268,11.76152,29.33024,-12.92289,-55.13187,-14.73791,7.52474,-102.05060,182.53144,-20.18960,-490237997.49400,486646248.63360,-781277018.26430,602300460.57290,-463787999.46420,249529525.96100,-123964746.86420,31353019.97807,-13353800.92159,-3463382.63269,-35469.17654,-1035343.45385,65076.64025,-38181.61312,-16473.33813,3928.44674,188.60263,1000.42530,-208376.39376,-700566.62363,114839.84613,-342407.71113,...
 39467.04812,-14553.84849,581895.26100,1012499.16715,406317.22416,310804.78515,43625.07033,-4157.26545,-5930.13795,-2070.62413,3348.17582,-10871.23729,144609.18550,60383.63650,27912.02226,15254.61228,-98561.37758,-67442.28158,-15573.63338,-19931.99773,24323.06905,-37473.32558,2840.64042,-8911.23694,-19636.31898,71725.21946,-12280.54554,12251.00101,88626.52260,5513.68450,18506.41546,-6513.87434,-83350.14621,44300.00743,-22075.37353,3731.57531,...
 -29167.76020,-21642.67384,56602666.72177,-22225578.01823,50576897.80669,-50319847.79086,5689259.25622,-29585299.79697,-4249711.27661,-4490830.29568,-727678.08724,366050.85631,19183.62792,55647.98226,1897.78091,-1091.03988,432.38158,-138.62556,101.38743,25.67379,320.20735,362.16615,85.06067,54.02616,...
 2.39460,18.70004,-8.43353,2.72100,-3.11205,-3.06201,136.31503,-28.33620,48.68781,-18.45285,1.15302,-1.52360,-0.13706,-0.37489,0.08750,-0.14579,-0.07051,-0.06518,0.30237,-0.00448,4.83172,6.83684,...
 1752447.78043,-945086.75857,2340978.12819,-1963675.42559,1254147.25257,-1274861.91191,279459.60237,-263954.01378,11835.62290,-16344.44434,9119.98960,-2390.44218,-23.67909,86.73916,-642.78635,-1290.12208,-0.43345,-1.85348,0.03094,-0.01655,0.12380,0.31834,5.54756,-1.63109,1.10598,-0.17578,2.66994,-2.17573,0.97360,-0.92226,...
 -0.18533,-0.39747,0.45809,-0.65286,0.26129,0.91922,0.81068,0.11183,6.32182,14.16786,0.20872,3.28489,-1.47458,-2.11724,1.70020,-1.99889,3.13003,1.90638,...
 -4483669.52795,-446454.90158,-6586256.67478,-671890.16779,-3620444.55554,-499672.41074,-855998.32655,-191073.94266,-73186.69110,-22649.38582,-2414.81729,-1294.40542,436.80907,125.48109,-81.16877,458.86508,-11.57414,-26.39114,-4.00801,-5.01054,-18.17569,20.86879,-4.80965,3.10535,-4.71122,1.18169,74.75544,649.21464,-26.55060,272.35592,-8.06982,16.86110,-26.54868,26.75711,-35.82910,38.51063,22.22814,19.38336,-6.30462,0.90602,0.62856,-0.34981,...
 -0.10232,-0.00939,0.04439,-0.18875,0.16025,0.11306,-0.06803,0.06004,-91305.66728,262370.61704,-194633.44577,304838.17733,-124505.90904,94111.75602,-22317.18255,1575.23438,748.66316,-349.78711,166.64450,-89.05045,120.76207,-100.26715,3.13806,3.71747,-1.44731,-0.35235,-0.51660,-1.50621,2.81310,-3.93573,1.20292,-0.36412,...
 -0.03340,-0.00561,-5.29764,26.02941,1.91382,3.30686,-3.35265,-3.20868,0.05807,-0.11885,-0.78588,0.34807,-0.19038,0.11295,-0.03642,-0.03794,0.00251,0.03449,-0.08426,-0.00310,0.05297,-0.09278,0.10941,0.00099,...
 -228688.56632,312567.73069,-331458.31119,328200.19460,-143760.57524,104182.01134,-17313.30132,12591.15513,-440.32735,-105.67674,104.35854,-852.84590,0.95527,0.30212,-54.63983,4.06948,0.07545,-0.13429,16.21005,29.24658,9.23410,50.48867,30.55641,12.76809,0.11781,0.70929,...
 -0.04100,13.60859,0.04976,-0.02083,0.36279,0.30130,-0.02129,0.09363,-0.07812,0.01570,-0.06217,-0.37181,-29348.55031,43889.87672,-35765.41577,33855.90070,-10128.69894,4535.32148,281.75353,-218.49194,-7.55224,134.28640,2.11319,-2.13109,15.71244,11.07183,-0.05406,-0.23337,-1.28949,1.34281,0.04212,-0.02080,0.08109,0.14820,-6010.46564,3639.41780,-5973.16000,1381.66999,-1177.36865,-501.06937,166.14792,-103.36431,14.92766,4.12877,-2.20893,-6.32033,-0.29038,-0.43172,-0.59341,0.20477,-0.13143,-0.03150,...
 0.10992,0.01976,-0.00254,0.02028,-0.30044,-0.44658,-0.03409,-0.10758,0.08349,0.06153,-0.06055,0.18249,-1.15341,-8.68699,-0.11348,-3.30688,...
 1.08604,1.04018,-0.46892,-0.69765,0.21504,0.01968,-0.00455,-0.01678,3.95643,-3.17191,3.95220,-2.12670,0.99305,-0.16651,0.34839,-0.49162,0.85744,0.20173,-0.00975,0.20225,...
 -0.02627,-0.02281,-0.18002,-0.01803,-0.06144,-0.21510,0.15935,-0.01251,-0.21378,0.44806,-0.01174,0.05779,0.07646,-0.19656,-0.04044,-0.02521,0.02996,0.06169,0.16698,-0.04710,-0.06506,-0.02114,0.05500,0.00276,0.08433,0.03160,0.08193,0.35773,0.05454,0.10718,...
 -0.02823,-0.00839,0.54078,0.49347,0.09609,0.11825,-0.16092,-0.11897,0.09059,0.08254,0.16712,0.05860,-0.09547,-0.03206,0.03876,0.04719,-0.02345,0.02240,-0.00609,-0.00649,0.03859,0.00077,0.47819,0.26196,0.09780,0.08104,-0.16919,0.05042,-0.42652,0.30810,...
 -0.03409,-0.51452,-0.23120,-0.01380,-0.01157,-0.00143,-0.00512,-0.01628,-0.00189,0.00183,-0.01427,-0.02861,0.00618,-0.00015,0.13087,0.13870,0.15158,-0.21056,-3.94829,-1.06028,-1.36602,0.77954,0.08709,-0.03118,-44.74949,91.17393,8.78173,45.84010,1.97560,-15.02849,-0.10755,-0.02884,3.38670,0.30615,130.92778,-24.33209,43.01636,-40.81327,-19.43900,22.18162,-0.12691,0.33795,...
 -6.44790,-6.23145,0.00319,0.01141,-0.03252,0.03872,0.04467,0.01614,-0.00382,-0.00019,0.05955,0.01533,16.11371,41.37565,61.44963,6.90615,1.41326,-0.73920,-0.03871,24.81978,-0.10229,-0.32775,-0.05188,-0.05628,...
 -2.33618,2.39053,-0.00584,0.00436,0.20903,0.02220,-0.01738,-0.02765,-0.00217,0.00613,-0.01772,0.01903,0.07075,-0.00530,0.15234,-0.37760,-0.11641,-0.20102,-0.63675,0.20525,-0.15783,0.58945,-0.06243,0.04306];
 
tabb=[ -35042727.30412,-49049197.81293,-25374963.60995,-5761406.03035,-467370.57540,14040.11453,2329.15763,-13978.69390,45.43441,29.70305,32.33772,-38.34012,26.43575,-28.76136,-18.59040,12.64837,5.56569,-12.51581,248.37350,-64.44466,54.02618,4.39466,-269.35114,-290.63134,-48.03841,-52.83576,...
 1508.94995,1682.78967,554.02336,715.65819,34.37602,58.44397,16.63685,16.10176,-1069.51609,2300.89166,-437.16796,927.89245,-33.17679,68.74495,18.72022,32.97640,-34004958.12619,-17758805.77098,-48416073.75788,-24973405.03542,-25374996.23732,-13351084.97340,-5738294.54942,-3082092.63350,-519989.39256,-206440.89101,44186.23548,-87639.22630,2506.47602,2327.01164,-53878.47903,-19670.13471,2.66934,-3.86086,106.32427,576.47944,46.56388,218.28339,...
 4.35402,15.04642,2.68717,-2.86835,0.81728,-2.34417,-1604.85823,-1999.24986,-631.47343,-1382.19156,-15.74075,-256.97077,6.99648,-4.54257,2.63478,1.88838,0.17628,-2.11518,-2.46735,-1.48743,1.83456,4.68487,-7.10919,3.57046,-5.36342,-7.70367,28395956.20816,-37176795.74372,48969952.83034,-48145798.96248,...
 31155823.23557,-21163596.14822,9057634.38260,-3167688.51696,1167488.70078,219103.97591,-19017.97335,107849.61195,-3814.43474,4405.92120,5800.13959,12619.88708,22.18168,-89.47801,52202.81929,55119.44083,5082.58907,37955.06062,-3165.24355,3316.67588,-113906.43970,-69279.41495,-57358.07767,-10176.17329,-4179.79867,2495.99374,787.87180,-154.35591,-1148.62509,1034.58199,-22194.95235,3341.97949,-4578.53994,108.30832,...
 7444.39789,16646.40725,509.75430,3808.92686,-179.85869,7408.76716,340.65366,1504.64227,-3783.09873,-13505.60867,875.74489,-3181.27898,-16220.93983,8041.37347,-2631.07448,2899.50781,18894.92095,-20072.81471,5925.05701,-1947.91902,-6731.56601,8014.52403,-987793.49463,6491762.34471,-279205.73643,6117135.96868,-140925.91402,2259422.06929,114028.61646,605600.90358,91858.00186,56506.65187,8949.15777,-9782.67413,-394.66541,-105.19208,-76.54752,-32.59411,-19.28741,10.40013,-107.64003,-7.36229,-22.25126,4.05952,...
 -3.74402,-2.79308,1.03337,-2.13968,1.53794,-0.02617,35.70756,12.97733,14.46213,6.20518,1.79381,1.65422,-0.31216,0.29053,-0.03538,-0.01584,-0.08934,0.00079,0.05539,-0.21591,2.86929,-2.24724,320797.07455,93342.16556,-20903.39115,79523.22083,-226588.37473,-121017.23944,-48472.25935,-74195.36778,...
 -7962.48081,-4607.76339,-4597.33274,-7983.12541,-20.34500,56.82999,-1038.19507,619.69624,1.08907,-0.91278,-0.13391,0.34956,-0.19982,-0.18296,-0.97688,2.36806,-0.30127,0.50980,0.96103,1.96432,0.43338,0.87317,0.36997,-0.01583,-0.44692,-0.25159,-0.53525,0.01154,-0.13231,0.35562,3.88928,-4.02882,1.06967,-0.56305,-0.45204,0.77213,-0.82873,-0.25854,0.21136,-1.06696,...
 458529.05491,616790.47568,698431.01349,1124501.41713,300226.10339,766533.33698,26896.22954,207880.75720,1116.29607,21793.26153,-850.64044,3528.95568,29.61278,-120.13367,376.95131,66.45758,-3.64868,2.76062,-0.85352,0.95115,5.35056,2.52803,0.90026,0.76403,0.43191,0.83605,125.81792,-39.65364,50.14425,-5.75891,...
 2.78555,2.05055,-4.27266,-4.92428,6.78868,5.73537,3.35229,-3.70143,0.08488,1.07465,0.10227,0.06074,...
 0.00291,0.01522,-0.02274,0.00297,0.01095,-0.01856,-0.02862,0.00178,143640.07486,707.21331,177163.08586,53386.52697,56856.89297,48268.74645,1764.52814,7711.76224,352.34159,-968.03169,-45.16568,-81.60481,...
 -76.35993,-98.06932,-1.42185,1.81425,-0.23427,0.59023,0.57127,-0.36335,1.89975,0.66890,0.28797,0.43592,-0.03769,0.03273,-6.06571,-2.68515,-0.55315,0.86977,1.53840,-0.59422,-0.05453,0.02447,-0.12658,0.22814,-0.01715,0.08497,...
 -0.01288,-0.00606,0.01547,-0.00692,0.01157,0.02407,-0.03883,0.00835,-0.01542,-0.04761,174386.39024,158048.26273,159192.81681,220154.55148,33716.11953,87537.86597,-116.90381,7535.83928,-962.06994,-132.28837,-644.90482,-110.52332,3.42499,3.74660,-0.94008,41.55548,...
 -0.03824,-0.05607,28.74787,-37.31399,30.87853,-26.11940,10.79742,-5.97905,1.01237,-0.04429,0.54402,0.41905,-0.02440,-0.03991,-0.00347,-0.04362,-0.00347,-0.00469,-0.02707,0.02761,-0.17773,-0.11789,26475.02580,35363.04345,19877.11475,41430.35940,2948.09998,12983.41406,281.93744,570.70054,147.83157,16.00090,-1.62814,-8.30846,...
 9.29131,-10.16496,-0.15799,0.03843,1.44716,0.46953,-0.02150,-0.02502,0.08861,-0.06690,2237.41551,3739.08722,753.74867,3460.41553,-298.69226,520.47031,-33.62615,-138.12767,3.61843,-8.29860,-4.56656,0.79553,0.20041,-0.25771,-0.35233,-0.27913,-0.02799,-0.08328,-0.06889,-0.16853,0.01701,-0.00964,-0.37737,0.18030,-0.08525,0.01906,...
 0.05236,-0.05155,0.11320,0.05991,-5.66926,-0.54402,-2.08508,-0.39407,0.82155,-0.55975,0.39168,-0.25551,0.00623,0.16162,-0.02519,0.02420,-1.23293,-3.19649,-0.60519,-2.79729,...
 0.05362,-0.61569,-0.25638,-0.27033,-0.03987,0.46623,-0.12070,0.00643,0.00849,-0.00768,-0.03687,0.10445,-0.13544,-0.00592,0.02078,0.09172,0.15824,0.15815,0.02020,0.00747,0.10919,0.09553,0.01953,-0.00135,...
 0.04266,-0.00218,0.02182,-0.13742,-0.01249,0.01724,-0.02200,0.02975,-0.01401,0.03416,-0.28873,0.04235,-0.08137,0.04223,-0.00326,0.02144,-0.40423,0.14281,-0.08256,0.02142,...
 0.08116,-0.03680,-0.02324,0.07260,-0.06746,0.11645,0.03233,-0.05997,-0.03101,0.02197,-0.00896,-0.00491,0.00574,0.00855,0.00052,0.01209,-0.31828,0.29955,-0.08133,0.04318,0.06787,-0.08865,-0.13228,-0.06507,0.34008,0.06417,-0.00177,-0.15116,...
 -0.00553,-0.01950,0.01144,-0.00309,-0.00115,-0.00153,0.02063,-0.00791,-0.00314,0.00493,-0.10614,0.08338,0.08845,0.20168,1.38955,-2.52285,-0.30475,-1.05787,0.00580,0.06623,-44.33263,-47.70073,-29.80583,-8.77838,7.02948,2.77221,0.05248,-0.13702,-0.78176,1.77489,-16.32831,46.68457,2.54516,21.78145,-5.09080,-8.42611,-0.24419,-0.03315,...