sat404.m

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

      0.16320,-1.55110,-1.44757,-3102.08749,-7452.61957,-5009.53858,-7216.29165,-2476.87148,-1880.58197,-574.49433,227.45615,144.50228,379.15791,225.36130,-443.47371,-8.51989,-3.75208,-4.25415,-1.59741, -0.43946,-0.06595,150.42986,6.54937,87.67736,92.32332,-21.97187,29.87097,-4.21636,-5.72955,-0.03879, -0.01071,-0.45985,0.02679,-0.02448,0.02397, -0.06551,-0.01154,1.97905,-0.82292,1.10140,0.30924, 0.03389,0.14230,0.00003,0.00119,-0.01117,...
      0.00665,-0.00132,-0.00576,-0.08356,0.08556,-0.26362,-0.12450, 0.00509,0.00165,0.02591,0.16200,-0.03318,0.06463, -0.00899,-0.00462,0.00102,0.00004,-0.73102,0.08299,-0.52957,-0.35744,0.14119,-0.24903,0.20843,0.14143,0.00031,-0.00234,-0.42643,-2.02084,1.58848,-1.57963,0.68418,2.07749,-0.45888,0.19859,-0.30277,-0.22591, 0.11607,-0.09705,0.00040,0.00431,-0.02683, 0.03158,-0.01302,-0.00541, 0.01742,-0.00006,-0.02231,-0.01128,-0.00800,...
      0.02055,-0.00346,0.00151,0.56732,-0.68995,0.27701,-0.16748, 0.01002,0.00043,0.26916,-0.57751,0.15547,-0.15825,-0.02074,-0.07722,-8.23483,-4.02022,0.69327,-5.91543,1.72440,1.02090,0.00024,-0.00053,20.03959,14.79136,76.43531,-14.42019,-7.82608,-69.96121,-54.94229,23.55140, 26.60767,14.68275,0.05118,-0.10401,-0.00075,-0.01942,-3.84266,-26.23442,10.20395,-14.77139,3.40853,2.07297,-0.53348,0.40635,0.00716,-0.00189,0.12472,-0.02903,...
      0.02254,-0.00183,-0.00175,-0.01522,0.00003,-0.00339,0.00383,-0.00168,0.01327,-0.03657,-0.08458,-0.00115,-0.03991,-0.02629, 0.00243,-0.00505,0.33875,-0.16744,0.05183,0.01744,  -0.24427,0.15271,0.37550,-0.17378,0.09198,-0.27966,-0.22160,0.16426,0.00032,-0.00310,-0.00022,-0.00144,-0.06170,-0.01195,-0.00918,0.02538, 0.03602,0.03414,-0.14998,-0.44351,0.45512,-0.11766,0.35638,0.27539,5.93405,10.55777,12.42596,-1.82530,-2.36124,...
      -6.04176,-0.98609,1.67652,-0.09271,0.03448,-0.01951,0.00108,0.33862,0.21461,0.02564,0.06924,0.01126,-0.01168,-0.00829,-0.00740,0.00106,-0.00854,-0.08404,0.02508,-0.02722,-0.06537,0.01662,0.11454,0.06747,0.00742, -0.01975,-0.02597,-0.00097,-0.01154,0.00164,-0.00274,0.02954,-0.05161,-0.02162,-0.02069,-0.06369,0.03846,0.00219,-0.01634,-0.04518,0.06696,1.21537,0.99500,0.68376,-0.28709, -0.11397,-0.06468,0.00607,-0.00744,0.01531,...
      0.00975,-0.03983,0.02405,0.07563,0.00356,-0.00018,-0.00009,0.00172,-0.00331,0.01565,-0.03466,-0.00230,0.00142,-0.00788,-0.01019,0.01411,-0.01456,-0.00672,-0.00543, 0.00059,-0.00011,-0.00661,-0.00496,-0.01986,0.01271,-0.01323,-0.00764,0.00041, 0.01145,0.00378,-0.00137,0.00652,0.00412,0.01946,-0.00573,-0.00326,-0.00257,-0.00225,0.00090,-0.00292,-0.00317,-0.00719,0.00468, 0.00245,0.00189,0.00565,-0.00330,-0.00168,-0.00047,...
      -0.00256,0.00220,0.00180,-0.00162,-0.00085,-0.00003,-0.00100,0.00098,-0.00043,0.00007,-0.00003,-0.00013];
        
tabr=[ -38127.94034,-48221.08524,-20986.93487,-3422.75861,-8.97362,53.34259,-404.15708,-0.05434,0.46327,0.16968,-387.16771,-146.07622,103.77956,19.11054,-40.21762,996.16803,-702.22737,246.36496,-63.89626,-304.82756,78.23653,-2.58314, -0.11368,-0.06541,-0.34321,0.33039,0.05652,-0.16493,67.44536,-29.43578,50.85074,18.68861,0.39742,13.64587,-1.61284,0.11482,0.01668,-0.01182,-0.00386,0.01025,0.00234,-0.01530,-0.02569,-0.00799,-0.00429,-0.00217,-0.00672,0.00650,...
      0.01154,0.00120,-0.00515,0.00125,0.00236,-0.00216,-0.00098,0.00009,-0.00460,-0.00518,0.00600,0.00003,0.00834,0.00095,0.01967,0.00637, -0.00558,-0.06911,-0.01344,-0.06589,-0.05425,-0.00607, -0.00247,-0.00266,0.08790,-0.08537,-0.00647,0.04028,-0.00325,0.00488,0.00111,-0.00044,-0.00731,0.00127,-0.00417,0.00303,0.05261,0.01858,-0.00807,0.01195,1.26352,-0.38591,-0.34825,1.10733, -0.02815,-0.02148,-0.05083,-0.04377,-0.01206,-0.00586,0.03158,-0.01117,0.00643,0.00306,...
      -0.01186,-0.05161,0.01136,-0.00976, -0.00536,0.01949,-1.41680,-0.81290,-0.09254,-0.24347, -0.14831,-0.34381,-2.44464,0.41202,-0.99240,-0.33707, -0.01930,-0.08473,0.00830,0.01165,-0.01604,-0.02439,0.00227,0.04493,-42.75310,-22.65155,-9.93679,-18.36179,2.73773,3.24126,-1.20698,1.07731,0.00434,-0.10360,-0.02359,0.00054,-0.02664,-0.00122,-19.79520,33.11770,-53.56452,-35.41902,67.95039,-82.46551,117.31843,14.08609,...
 0.06447,0.03289,0.40365,-0.33397,0.07079,-0.09504,-30.36873,6.23538,-14.25988,-44.91408,38.53146,-16.31919,6.99584,22.47169, -0.13313,0.28016,6.83715,-6.01384,1.68531,-3.62443,-0.22469,-0.29718,0.25169,0.13780,-3.64824,1.22420,-2.48963,-1.12515, -0.01510,-0.56180,-0.03306,0.01848,-0.00103,-0.00077,-0.01681,-0.00227,-0.00402,-0.00287,0.04965,-0.16190,-0.40025,0.20734, 0.15819,-0.25451,0.02467,-0.00495,0.00597,0.00490,-0.01085,-0.00460,-0.71564,-0.26624,0.03797,-0.28263, ...
 0.03510,0.30014,2.79810,0.07258,-0.01618,0.00337,0.00876,0.04438,0.00742,-0.00455,-0.01163,-0.00683,0.00950,0.01275,-0.02124,-0.67527,-0.23635,0.06298,-0.03844,0.01010,0.73588,-0.00271, 0.01742,-0.00467,0.00017,-0.00505,-0.27482,5.00521,-1.92099,1.55295,-0.35919,-0.09314,-0.47002,0.06826,0.07924,0.16838,-0.04221,0.71510,-0.16482,0.08809,41.76829,-125.79427,106.65271,-71.30642,36.18112,17.36143,-1.63846,5.02215,-1.08404,0.00061,2.45567,-2.42818,-9.88756,5.36587,...
 -0.61253,-0.35003,1523.54790,602.82184,68.66902,1878.26100,-1098.78095,-120.72600,127.30918,-383.96064,-7.00838,-6.09942,-1.54187,0.34883,-9.47561,-4.35408,-21541.63676,-32542.09807,-29720.82604,-28072.21231, -15755.56255,-8084.58657,-8148.87315,7434.89857,11033.30133,7827.94658,610.18256,-11411.93624,-9.87426,0.94865,-1.63656,0.41275,1996.57150,511.48468,669.78228,1363.67610,-379.72037,198.84438,-16.63126,-79.37624,-2.30776,-246.07820,-16.85846,-148.18168,-6.89632,...
 -20.49587,0.39892,-0.34627,-57.81309,-136.96971,15.25671,-96.61153,16.09785,-8.79091,0.70515,1.16197,0.05647,0.04684,0.25032,-0.19951,0.07282,-0.00696,0.00493,0.00733,-0.01085,0.00422,-0.01309,0.00262,0.37616,-0.36203,-0.11154,0.18213,0.15691,0.29343,0.00485,0.06106,-0.01492,0.09954,0.28486,2.27190,0.33102,1.50696,...
 -0.01926,0.04901,0.01827,0.00863,-0.03315,0.00178,-0.77600,-0.48576,-0.21111,-0.19485,1.90295,6.44856,1.71638,2.12980,-7.19585,-0.08043,0.07004,-0.02764,0.01604,0.01158,0.00936,-0.01199,0.18396,-0.29234,0.10422,-0.00720,0.05196,0.10753,0.02859,-0.03602,0.63828,1.96280,-0.31919,0.85859,-0.10218,-0.00673,0.01748,-0.02190,0.01266,-0.02729,-4.80220,8.90557,-5.94059,2.28577,-0.19687,-1.28666,0.32398,0.14879,...
 -0.02619,-0.02056,-0.04872,-0.07011,-0.04082,-0.04740,0.60167,-2.20365,-0.27919,-0.45957,-1.31664,-2.22682,176.89871,13.03918,0.00568,0.00560,0.01093,0.00486,-0.00948,-0.31272,-11.87638,-3.68471,-1.74977,-9.60468, 2.94988,-0.57118,0.00307,-0.01636,0.02624,0.03032,-0.00464,-0.01338,0.00935,0.00530,-0.11822,0.03328,-0.41854,0.04331,0.41340,-0.21657,-0.00865,0.00849,-0.00374,-0.00899,0.01227,-0.23462,-0.71894,-0.04515, 0.00047,0.28112,-0.12788,0.11698,-0.02030,0.02759,0.02967,-0.00092,0.00454,0.00565,-0.00026,0.00164,-0.01405,-0.00862,0.01088,0.05589,0.18248,-0.06931,-0.00011,0.03713,...
 0.01932,-0.00982,-0.13861,0.09853,-0.03441,-0.02492,2.26163,-5.94453,4.14361,-0.94105,0.39561,0.75414,-0.17642,0.03724,-1.32978,-0.56610,-0.03259,-0.06752,39.07495,80.25429,-28.15558,82.69851,-37.53894,-17.88963,6.98299,-13.04691,-0.48675,-1.84530,-0.07985,-0.33004, -3.39292,2.73153,-17268.46134,1144.22336,-16658.48585,5252.94094,-3461.47865,2910.56452,-433.49442,-305.74268,-383.45023,545.16136,313.83376,27.00533,-31.41075,7.90570,-12.40592,3.01833,-0.83334,0.23404,59.26487,-112.74279,113.29402,-15.37579,14.03282,32.74482,-4.73299,1.30224,-0.00866,0.01232,-0.53797,0.00238,-0.07979,0.04443,...
 -0.05617,-0.05396,0.10185,-1.05476,0.43791,-0.32302,0.06465,0.03815,0.00028,-0.00446,0.09289,-0.06389,0.01701,-0.01409,0.47101,0.16158,0.01036,-0.39836,0.00477,0.01101,-2.06535,0.33197,-0.82468,-0.41414,0.03209,-0.09348,0.00843,-0.00030,-9.49517,-3.82206,0.66899,-10.28786,6.33435,1.73684,-0.98164,2.25164,-0.07577,-0.00277, 1.02122,0.75747,1.79155,-0.77789, -2.56780,-2.07807,0.19528,0.77118,-0.28083,0.32130,-0.04350,-0.07428,-0.01161,0.01387,0.02074,0.19802,-0.03600,0.04922,-0.19837,0.02572,-0.00682,-0.04277,-0.01805,0.00299,0.03283,-0.02099,3.57307,1.17468,0.65769,1.88181,...
 -0.39215,0.08415,-0.53635,-0.19087,-0.12456,0.02176, 0.01182,-0.07941,-2.43731,2.44464,1.03961,-1.81936, 30.33140,0.92645,0.00508,-0.01771,-81.06338,66.43957,33.16729,131.44697,76.63344,-34.34324,-35.33012,-28.04413,-1.47440,13.09015,0.13253,-0.01629,0.02187,-0.00963,-21.47470,-9.44332,-7.21711,-12.59472,1.76195,-1.63911,0.09060,0.28656,0.00635,0.00536,0.03470,-0.06493,0.00666,-0.01084,0.01116,-0.01612,-0.00102,0.00208,...
 -0.05568,0.00628,0.02665,-0.01032,0.21261,-1.90651,0.72728,-0.57788,0.08662,0.10918,3.39133,3.97302,-4.63381,4.26670,-2.50873,-3.76064,1.28114,1.81919,1.48064,-0.37578,-0.26209,-0.47187,0.00282,-0.00499,0.01749,0.03222,1.60521,-1.79705,1.61453,0.68886, -0.29909,0.55025,-0.07894,0.19880,-0.15635,0.46159,2.09769,1.52742,-7.60312,11.34886,4.35640,8.61048,2.15001,-2.15303,-0.61587,-0.11950,-0.03289,-0.00520,-0.00501,-0.00445,0.15294,-0.05277,0.02455,0.00408,1.19601,0.43479,0.20422,0.57125,-0.12790,0.01318,-0.15275,-0.43856,6.99144,-0.08794,-1.69865,0.82589,-0.20235,0.97040,...
 0.20903,0.00675,0.26943,0.08281,0.03686,0.05311,1.28468,1.21735,-1.38174,1.29570,-0.75899,-1.17168,0.44696,-0.32341,-0.06378,-0.27573,-0.06406,0.87186,0.21069,0.19724, 0.00119,-0.04147,0.39279,0.51437,-0.11035,0.21450,-0.04309,0.02359,0.20490,0.14210,0.00007,-0.00017,-0.03529,-0.02644,0.10710,0.44476,-0.02632,-0.01817, 2.11335,-0.04432,0.18206,0.27335,0.08867,0.00313, -0.00692,0.01595,-0.72957,0.32080,-0.29291,-0.44764,0.12767,-0.05778,0.04797,-0.00223,0.17661,0.22427,-0.04914,0.09114,0.12236,0.00708,0.74315,-0.01346,0.02245,-0.02555,-0.30446,0.13947,-0.12340,-0.18498,...
 -0.04099,0.02103,0.06337,-0.01224,0.28181,-0.01019,-0.02794,-0.09412,0.03272,-0.01095,0.11247,-0.00650,-0.01319,-0.04296,0.04653,-0.00423,0.06535,0.00014];

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

%Total terms = 215, small = 211 

sat.maxargs=9;
sat.max_harmonic=[ 0,  0,  1,  0,  8, 18,  9,  5,  0];
sat.max_power_of_t=7;
sat.arg_tbl=args;
sat.lon_tbl=tabl;
sat.lat_tbl=tabb;
sat.rad_tbl=tabr;
sat.distance=  9.5575813548599999e+00;
sat.timescale=3652500.0;
sat.trunclvl=1.0;

ans=sat;