sat404.m

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

function ans=sat404

%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:   0.78   0.26   0.55 
% -2499.7 to  -1999.7:   0.66   0.19   0.57 
% -1999.7 to  -1499.7:   0.62   0.19   0.53 
% -1499.7 to   -999.8:   0.79   0.17   0.61 
%  -999.8 to   -499.8:   0.78   0.15   0.42 
%  -499.8 to      0.2:   0.75   0.19   0.52 
%     0.2 to    500.2:   0.62   0.18   0.41 
%   500.2 to   1000.1:   0.56   0.13   0.54 
%  1000.1 to   1500.1:   0.53   0.15   0.41 
%  1500.1 to   2000.1:   0.51   0.15   0.49 
%  2000.1 to   2500.0:   0.52   0.13   0.41 
%  2500.0 to   3000.0:   0.63   0.22   0.53 
%  3000.0 to   3000.4:  0.047  0.073  0.086 


tabl=[ 1788381.26240,2460423.68044,1370113.15868,415406.99187,72040.39885,12669.58806,439960754.85333,180256.80433,18.71177,-40.37092,66531.01889,-195702.70142,57188.02694,-179110.60982,-19803.06520,-58084.15705,-9055.13344,-31146.10779,11245.43286,-3247.59575,459.48670,2912.82402,...
 -4.06749,-13.53763,-30.55598,-4.51172,1.48832,0.37139,597.35433,1193.44545,-297.50957,976.38608, -263.26842,34.84354,-6.77785,-29.92106, -0.16325,-0.18346,-0.15364,-0.08227,0.20180,0.02244,0.04672,-0.29867,-0.04143,-0.00760,-0.17046,-0.00778,0.04200,0.23937,-0.00098,-0.05236,...
 -0.02749,-0.01813,0.00637,0.01256,-0.04506,0.04448,-0.00105,0.06224,0.01157,0.17057,-0.03214,0.18178, -0.22059,-0.01472,-0.24213,0.04309,0.03436,0.44873, 0.01350,-0.01931,-0.80618,-0.56864,0.29223,-0.03101, 0.04171,0.02264,-0.01264,-0.01645,0.01774,0.06374,-0.01925,-0.03552,0.10473,-0.04119,0.08045,0.04635,-3.01112,-9.26158,8.13745,1.88838,...
 -0.15184,0.16898,-0.22091,0.29070,-0.03259,0.06938, -0.08499,-0.21688,0.01848,-0.05594,0.50100,-0.00027,0.13300,0.12055,0.03039,0.03854,-1.55287,2.55618,-0.45497,-0.29895, -0.93268,0.83518,-0.32785,7.03878,-1.66649,2.75564, -0.29459,0.01050,0.08293,-0.03161,-0.12750,-0.04359,0.04217,0.07480,-114.43467,49.47867,-66.52340,-26.27841, 15.48190,-13.06589,3.28365,5.02286, -0.17155,-0.07404,0.00924,-0.07407,-0.02922,0.06184,108.04882,86.09791,-155.12793,208.10044,...
 -311.72810,-268.92703,74.57561,-420.03057, -0.07893,0.09246,-0.66033,-0.39026,-0.13816,-0.08490, -36.79241,-78.88254,71.88167,-68.05297,51.71616,65.77970,-43.59328,23.51076,-0.02029,-0.32943, -8.82754,1.48646,-3.12794,2.12866,-0.06926,0.44979, 0.00621,-0.51720,-3.82964,-1.48596,-0.11277,-3.21677, 0.81705,-0.19487,-0.06195,0.10005,-0.02208,0.00108, 0.00455,-0.03825,0.01217,-0.00599,-0.17479,-0.47290,0.85469,1.12548,...
 -0.80648,-0.44134,-0.01559,-0.07061,0.01268,-0.01773,0.01308,-0.03461,-0.71114,1.97680,-0.78306,-0.23052, 0.94475,-0.10743,0.18252,-8.03174,0.00734,0.04779,0.12334,-0.03513,0.01341,0.02461,0.02047,-0.03454,0.02169,-0.01921,-1.12789,0.09304,0.14585,0.36365, 0.03702,0.10661,-0.00464,-1.72706, -0.00769,-0.04635,-0.01157,0.00099,10.92646,1.96174,2.91142,4.74585, -0.29832,0.75543,0.05411,1.05850,0.38846,-0.16265,1.52209,0.12185,0.18650,0.35535, -278.33587,-82.58648,-160.00093,-225.55776,35.17458,-77.56672,10.61975,3.33907, 0.06090,2.17429,-4.32981,-5.84246,11.43116,20.61395,-0.65772,1.28796,1224.46687,-3113.15508,3798.33409,-137.28735,-256.89302,2227.35649,-779.78215,-260.37372,11.73617,-13.25050,-0.75248,-2.87527,-8.38102,17.21321,...
 -61784.69616,39475.02257,-54086.68308,54550.85490,-16403.69351,29602.70098,14672.06363,16234.17489,15702.37109,-22086.30300,-22889.89844,-1245.88352,1.48864,19.75000,0.78646,3.29343,-1058.13125,4095.02368,-2793.78506,1381.93282,-409.19381,-772.54270,161.67509,-34.15910, -514.27437,27.34222,-311.04046,48.01030,-43.36486,16.19535,-0.73816,-0.81422,287.32231,-110.44135,200.43610,37.98170,17.73719,34.40023,-2.46337,1.48125, 0.09042,-0.11788,0.37284,0.51725,0.00597,0.14590, -0.01536,0.00980,0.00721,0.02023,0.00027,0.02451,-0.72448,-0.71371,0.29322,0.18359,0.72719,-0.37154,0.14854,-0.02530,0.23052,0.04258,...
 4.82082,0.01885,3.11279,-0.63338, 0.10559,-0.02146,-0.01672,0.03412,0.00605,0.06415,-0.89085,1.51929,-0.36571,0.39317, 12.05250,-3.79392,3.96557,-3.51272,-0.17953,12.30669,-0.05083,-0.11442,0.02013,-0.02837,-0.02087,-0.01599,0.49190,0.30360,0.01316,0.17649,0.21193,-0.09149,-0.07173,-0.05707,4.24196,-1.25155,1.81336,0.68887, -0.01675,0.20772,-0.04117,-0.03531,-0.02690,-0.02766,37.54264,10.95327,8.05610,30.58210,-12.68257,1.72831,0.13466,-3.27007, 0.01864,...
 -0.00595,0.03676,0.14857,-0.07223,0.06179,0.44878,-1.64901,-20.06001,0.63384,-4.97849,4.78627,29.87370,7.29899, 0.00047,-0.00155,0.00314,0.01425,-0.17842,-0.08461,-1.61020,-8.47710,6.85048,-4.38196, 1.05809,2.68088,-0.01027,-0.00833,0.06834,-0.04205,0.03330,-0.01271,0.01301,-0.01358,0.03537,0.03612,0.02962,0.62471,-0.30400,-0.64857,...
 0.01773,0.01890,0.01426,-0.00226,-0.50957,-0.01955,-0.09702,1.09983,0.64387,-0.02755,0.26604,0.30684,0.06354,0.05114,-0.00058,-0.04672, -0.00828,0.00712,-0.00440,0.00029,-0.01601,0.03566,0.13398,-0.02666,-0.06752,-0.43044,0.07172,-0.01999,-0.01761,-0.05357, 0.06104,0.29742,-0.08785,0.05241,-6.57162,-4.20103,0.03199,-6.46187,1.32846,-0.51137,0.06358,0.37309,-1.46946,2.34981,-0.18712,0.11618,240.62965,-107.21962,219.81977,84.04246,...
 -62.22931,68.35902,-9.48460,-32.62906,5.57483,-1.82396,1.00095,-0.39774,7.87054,11.45449,-432.67155,55064.72398,12444.62359,54215.28871, 8486.03749,12297.48243,-333.27968,1147.93192,1403.73797,990.40885,-3.84938,-722.43963,16.83276,96.48787,7.04834,38.22208,0.63843,2.61007, 230.73221,171.64166,1.96751,287.80846,-85.21762,31.33649,-2.25739,-11.28441,0.04699,0.06555,-0.08887,1.70919,0.09477,0.26291, -0.15490,0.16009,1.93274,1.01953,0.36380,1.29582, -0.13911,0.14169,-0.00491,-0.00030,-0.08908,-0.10216,-0.03265,-0.03889,0.40413,-1.12715,-0.94687,-0.04514, 0.02487,-0.01048,0.39729,2.82305,-0.61100,1.11728, -0.13083,-0.04965,-0.00602,-0.02952,-6.13507,13.73998,-15.70559,-1.28059,2.64422,-9.33798,3.26470,1.56984,-0.00572,0.09992,...
 -8.80458,-8.23890,-11.51628,9.47904,11.31646,4.29587,-2.41367,-0.05883,-0.80022,-1.02706,0.21461,-0.06864,0.01882,0.01798,0.27614,-0.01007,0.04362,0.07560,0.05519,0.23435,-0.09389,0.01613,0.01298,0.04691,-0.02665,-0.03582,0.60080,-4.28673,1.87316,-1.05840, 0.13248,0.40887,-0.67657,0.67732,0.05522,0.07812, -0.17707,-0.07510,0.24885,10.63974,-7.40226,-2.33827, 2.75463,-32.51518,0.05140,0.01555,180.43808,263.28252,384.50646,-76.53434,-93.50706,-220.50123,-81.91610,103.92061,...
 30.90305,-2.89292,-0.06634,-0.37717,-0.01945,-0.05936,29.27877,-59.73705,35.86569,-18.36556,3.88812,4.82090,-0.70903,0.06615,0.01558,-0.01854,0.16209,0.12682,0.02508,0.02406,-0.03078,-0.01737,-0.00033,-0.00020,0.01023,0.05972,-0.03373,-0.07289,-2.08162,-0.14717,-0.64233,-0.75397, 0.11752,-0.09202,4.42981,-4.19241,5.02542,5.03467, -4.22983,2.80794,3.03016,-2.74373,-1.11490,-2.72378,-0.63131,0.74864, -0.00759,-0.00675,0.03615,-0.01806,-2.71920,-1.50954,0.54479,-1.92088, 0.66427,0.32228,-2.55188,-0.65332,-2.73798,2.10182, 1.54407,3.01357,38.76777,23.54578,27.29884,-14.93005,...
 -7.50931,-5.66773,0.30142,1.52416,0.00634,0.09697,-0.00748,0.01433,0.02936,0.53228,-0.03603,0.06345,0.30816,-1.07925,0.46709,-0.21568, 0.01663,0.10810,-0.42511,0.35872,-0.19662,-6.74031,1.05776,1.86205,1.08919,0.10483,-0.03368,-0.21535, 0.07556,-0.27104,0.05142,-0.03812,1.20189,-1.36782,1.35764,1.39387,-1.19124,0.77347,-0.54760,-0.26295,-0.07473,0.23043, 2.82621,-0.23524,0.47352,-0.81672,-0.08515,0.04700,0.55355,-0.40138,0.22255,0.12236, -0.09110,0.31982,0.39404,-0.17898,-0.00056,0.00014,-0.02012,0.03102,0.43236,-0.10037,-0.00961,0.07440, -0.07076,-1.97272,0.25555,-0.21832,-0.00837,-0.08393,...
 0.01531,0.00627,0.33193,0.70765,-0.43556,0.28542, -0.23190,-0.04293,-0.08062,0.13427,0.23763,-0.17092,0.09259,0.05155,0.08065,-0.11943,-0.02174,-0.68899,-0.01875,-0.01746,0.13604,0.29280,-0.17871,0.11799, 0.02003,0.04065,0.01343,-0.06060,-0.01290,-0.26068,-0.09033,0.02649,-0.00092,-0.03094,-0.00770,-0.10447,-0.04113,0.01259,-0.00469,-0.04346,-0.00010,0.06547];

tabb=[ -567865.62548,-796277.29029,-410804.00791,-91793.12562,-6268.13975,398.64391,-710.67442,175.29456,-0.87260,0.18444,-1314.88121,20709.97394,-1850.41481,20670.34255,-896.96283,6597.16433,-179.80702,613.45468,17.37823,-13.62177,-0.36348,12.34740,0.47532,0.48189,0.27162,-0.20655,-0.23268,0.05992,46.94511,15.78836,21.57439,23.11342,-0.25862,5.21410,-0.22612,-0.05822,-0.00439,-0.01641,-0.01108,-0.00608,0.00957,0.00272,-0.00217,0.00001,...
      -0.00534,-0.00545,0.00277,-0.00843,0.00167,-0.00794,0.00032,-0.00242, -0.00002,-0.00041,-0.00025,0.00031,0.00062,-0.00060,0.00083,0.00032,0.00527,-0.00211,0.00054,0.00004, -0.02769,-0.01777,0.00247,0.00097,0.00020,-0.00232,0.00044,-0.00035,      -0.00072,0.01341,0.00325,-0.01159, 0.00079,-0.00078,-0.00009,0.00066,0.00222,0.00002,0.00013,-0.00161,0.01374,-0.05305,0.00478,-0.00283,0.16033,0.13859,0.33288,-0.16932,-0.00316,0.00625,...
      -0.00309,0.01687,0.00001,0.00486,0.00401,-0.01805,-0.00048,-0.00407,-0.01329,0.01311,-0.00591,0.00166,0.00830,0.00665,-0.80207,0.22994,-0.34687,0.08460,-0.11499,-0.01449,-0.01574,0.78813,-0.03063,0.28872,-0.00337,0.01801,-0.01703,-0.00929,-0.00738,0.03938,0.05616,-0.00516,-3.09497,30.13091,-3.14968,17.62201,-0.73728,2.46962,-0.11233,0.03450,-0.07837,-0.01573,-0.01595,0.00394,0.00174,0.01470,6.83560,-2.37594,4.95125,3.24711,2.44781,...
      5.17159,1.99820,-2.38419,0.00840,0.03614,-0.00209,-0.30407,      -0.02681,-0.06128,1.50134,11.82856,4.39644,6.98850,-4.17679,5.73436,-9.66087,1.98221,-0.29755,0.08019,-0.24766,-8.54956,-1.74494,-3.36794,-0.32661,-0.00722,0.14141,0.01023,-1.21541,-2.58470,0.38983,-1.70307,0.31209,-0.10345,0.02593,0.02178,0.00289,0.00393,-0.00236,-0.00373,-0.00270,-0.00049,-0.06282,-0.00443,-0.02439,-0.02254,  -0.02220,0.03532,-0.00072,0.00010,...
      -0.00049,-0.00112,0.00086,0.00112,0.10135,-0.10972,0.08357,0.00155, 0.04363,-0.00201,-0.01996,-0.01341,-0.00039,-0.00042,-0.00294,0.00070,0.00005,-0.00027,0.00070,-0.00076,0.00234,-0.00239,-0.08365,-0.08531,-0.03531,0.15012,-0.01995,-0.01731,-0.00370,-0.00745,-0.00315,-0.00079,-0.00120,-0.00145,-0.99404,-1.31859,0.03584,-0.83421, 0.10720,-0.05768,0.06664,-0.09338,-0.01814,-0.00003,-0.05371,-0.06458,-0.00100,-0.01298,-7.08710,...
      -23.13374,4.18669,-19.94756,4.85584,-3.37187,0.58851,0.31363,0.01994,0.27494,-1.37112,2.61742,0.52477,      -0.46520, -0.13183,0.26777,836.90400,-484.65861,815.99098,236.54649,-32.38814,288.95705,-68.17178,-18.87875,-1.79782,-3.68662,-1.27310,-0.65697,-3.67530,2.10471,-13758.97795,4807.62301,-14582.14552,9019.73021,-3202.60105,4570.16895,2078.68911,2892.62326,-2399.35382,3253.16198,-8182.38152,-3588.77680,-0.16505,1.08603,...
      0.53388,0.87152,61.53677,538.43813,-407.32927,322.27446, -148.71585,-179.37765,54.07268,-34.12281,-14.76569,-17.95681,-10.82061,-6.39954,-2.10954,0.67063,0.22607,-0.43648,20.90476,-45.48667,30.39436,-14.20077,5.17385,5.12726,-0.66319,0.55668,0.02269,-0.00016,0.07811,0.00111,0.01603,0.01020,-0.00107,0.00494,-0.00077,-0.00084,-0.00196,0.00081,-0.03776,0.01286,-0.00652,-0.01450,0.05942,-0.08612,0.01093,-0.01644,0.02147,-0.00592,...
      0.36350,-0.00201,0.14419,-0.10070, -0.00491,-0.01771,-0.00053,-0.00033,0.00146,0.00048,0.00582,0.04423,-0.00549,0.00983,0.27355,-0.38057,0.24001,-0.05441,-0.07706,0.14269,-0.00059,-0.00154,-0.00013,-0.00088,-0.00046,0.00029,-0.00276,-0.00507,0.00075,-0.00076,0.01806,0.00862, -0.00510,-0.01364,-0.00029,-0.12664,0.03899,-0.03562, 0.00318,0.00514,0.00057,0.00201,0.00028,0.00014,-0.47022,-0.74561,0.40155,-0.16471,-0.18445,0.34425,...
      -0.07464,-0.13709,-0.01018,-0.00748,-0.01210,-0.04274,-0.00579,-0.00692,-11.09188,-1.67755,-6.62063,-13.84023,12.75563,-6.73501,8.31662,5.40196,0.00052,0.00034,0.00128,0.00085,-0.02202,-0.00599,-0.33458,-1.65852,1.47003,-1.02434, 0.87885,1.15334,-0.00241,-0.00721,0.03154,0.00612,0.00318,-0.02521,0.00042,0.00213,-0.01094,0.05417,-0.03989,-0.00567,0.00123,-0.00244,0.00108,0.00242,-0.00138,-0.00099,0.04967,0.01643,-0.00133,0.02296, ...
      0.12207,0.05584,0.00437,-0.04432,-0.00176,-0.00922, -0.00252,0.00326,-0.00020,-0.00050,-0.00263,-0.00084,-0.01971,0.00297,0.03076,0.01736,-0.01331,0.01121,-0.00675,0.00340,-0.00256,0.00327,-0.00946,0.03377,-0.00770,0.00337,0.61383,0.71128,-0.02018,0.62097, -0.07247,0.04418,-0.02886,-0.03848,-0.44062,0.03973,-0.00999, -0.04382,57.94459,117.45112,-71.22893,126.39415,-62.33152,-31.90754,12.17738,-16.46809,-1.13298,0.08962,-0.20532,...