-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmercury-900-1c-plot.py
12 lines (10 loc) · 16.2 KB
/
mercury-900-1c-plot.py
1
2
3
4
5
6
7
8
9
10
11
12
import matplotlib.pyplot as plt
for year, r_max, u_min, x, y in [(1, 47383274.116350025, 0.00012947307541310896, [47279994.64428885, 47202795.33765246, 46971500.75688101, 46587020.6035792, 46050867.640217535, 45365151.7183375, 44532571.45629386, 43556403.59943748, 42440490.1048493, 41189223.00137663, 39807527.08488578, 38300840.5168364, 36675093.4027919, 34936684.43511126, 33092455.69223043, 31149665.693617843, 29115960.816598054, 26999345.187831312, 24808149.16801832, 22550996.55404874, 20236770.627767723, 17874579.18534265, 15473718.68466145, 13043637.65239746, 10593899.49458748, 8134144.857410505, 5674053.686316681, 3223307.1329371133, 791549.4597808219, -1611649.9073271349, -3976834.0296973246, -6294695.660951313, -8556113.860583164, -10752189.87187545, -12874282.12311519, -14914040.214595005, -16863437.75779436, -18714803.93780399, -20460853.67508916, -22094716.268224653, -23609962.405170932, -25000629.437238224, -26261244.816616952, -27386847.6055029, -28373007.972578794, -29215844.60044934, -29912039.93562573, -30458853.22143961, -30854131.262485135, -31096316.878528778, -31184455.014506105, -31118196.482556608, -30897799.321386162, -30524127.76738736, -29998648.84164798, -29323426.56606657, -28501113.831415303, -27534941.949128635, -26428707.928175956, -25186759.52700825, -23813978.139588818, -22315759.58295231, -20697992.86210015, -18967036.99609464, -17129695.9966938, -15193192.0982511, -13165137.344425768, -11053503.64385146, -8866591.41266356, -6612996.927471864, -4301578.5173543105, -1941421.7280135648, 458196.40485142067, 2887843.832384018, 5337970.576141291, 7798946.272604263, 10261098.03105425, 12714748.455993885, 15150253.684771597, 17558041.29096807, 19928647.904630706, 22252756.40170323, 24521232.51644159, 26725160.73326539, 28855879.317010976, 30905014.34406184, 32864512.60077633, 34726673.22013037, 36484177.932335965, 38130119.81068818, 39658030.399897635, 41061905.12042885, 42336226.84899522, 43475987.58296369, 44476708.10348314, 45334455.560236126, 46045858.908751205, 46608122.13997849, 47019035.25004673, 47276982.90714966, 47380950.78202721], [0, 2406960.4240300623, 4804466.736305154, 7183080.054842239, 9533435.964368118, 11846281.380279185, 14112510.97273811, 16323203.007579194, 18469654.462851606, 20543415.282838684, 22536321.634742595, 24440528.037039474, 26248538.233198084, 27953234.689108152, 29547906.598274685, 31026276.284520905, 32382523.898339514, 33611310.30977041, 34707798.10777805, 35667670.6235128, 36487148.90264915, 37163006.5602864, 37692582.460040346, 38073791.16733715, 38305131.13615685, 38385690.59724913, 38315151.12476166, 38093788.86735335, 37722473.43947326, 37202664.47728115, 36536405.87297237, 35726317.71066836, 34775585.93580036, 33687949.79885838, 32467687.12347703, 31119597.45684862, 29648983.169265576, 28061628.57713627, 26363777.172065794, 24562107.045674976, 22663704.60697784, 20676036.695731502, 18606921.201598205, 16464496.304551909, 14257188.457576018, 11993679.237363836, 9682871.193304265, 7333852.828937554, 4955862.853395409, 2558253.8432762846, 150455.45770937268, -2258062.6488064006, -4657828.070027827, -7039402.901593918, -9393420.858129876, -11710624.107771851, -13981899.679127075, -16198315.297381176, -18351154.50837803, -20431950.952455863, -22432521.652991228, -24344999.188688386, -26161862.62289775, -27875967.068336375, -29480571.770848542, -30969366.60189128, -32336496.85559492, -33576586.25304956, -34684758.06355479, -35656654.25996256, -36488452.63312509, -37176881.79833441, -37719234.03519742, -38113375.910713226, -38357756.643947825, -38451414.18006984, -38393978.94976932, -38185675.29986129, -37827320.58963626, -37320321.95675006, -36666670.765515, -35868934.759843454, -34930247.95156614, -33854298.284384705, -32645313.121826403, -31308042.616500624, -29847741.026069343, -28270146.049537912, -26581456.265074648, -24788306.75815929, -22897743.03587601, -20917193.329925653, -18854439.397220146, -16717585.93284678, -14515028.715652417, -12255421.611552514, -9947642.564361284, -7600758.707687969, -5223990.735051262, -2826676.66819166, -418235.16595541715]), (300, 83953957.865674, 9.545149352207785e-05, [83812258.1048585, 83678557.28912325, 83278180.5459894, 82612703.18876709, 81684744.27205463, 80497956.25871424, 79057010.61931075, 77367579.41912018, 75436312.96653338, 73270813.61042805, 70879605.79019774, 68272102.4559076, 65458567.99169996, 62450077.7880645, 59258474.62237524, 55896322.01960514, 52376854.777207464, 48713926.84877988, 44921956.792095505, 41015870.9964473, 37011044.912962794, 32923242.51947551, 28768554.25828428, 24563333.69170907, 20324133.124513913, 16067638.447204866, 11810603.4566395, 7569783.912747666, 3361871.5910189287, -796571.4095024857, -4889177.844130231, -8899839.729767917, -12812771.72795829, -16612573.254419988, -20284289.070691608, -23813468.119721163, -27186220.374115806, -30389271.473633807, -33410014.937402505, -36236561.74578999, -38857787.097299516, -41263374.15732032, -43443854.62682752, -45390645.971934564, -47096085.168521605, -48553458.8293596, -49757029.595797256, -50702058.6900548, -51384824.53981972, -51802637.40191602, -51953849.92684599, -51837863.62349351, -51455131.197411485, -50807154.75371196, -49896479.871044785, -48726685.5701004, -47302370.215735696, -45629133.40794757, -43713553.933195025, -41563163.86291352, -39186418.90101555, -36592665.09774748, -33792102.0611297, -30795742.811046917, -27615370.434543412, -24263491.71335575, -20753287.90655423, -17098562.88241809, -13313688.804145662, -9413549.583235493, -5413482.32342508, -1329216.9857278918, 2823185.4879101785, 7027396.348882773, 11266883.265586356, 15524975.34390126, 19784928.68881099, 24029992.249443557, 28243473.68875144, 32408805.018956352, 36509607.74477013, 40529757.25852191, 44453446.23420339, 48265246.77132555, 51950171.0443584, 55493730.21977416, 58881991.409055375, 62101632.43398215, 65139994.18897889, 67985130.39521909, 70625854.55061626, 73051783.89164048, 75253380.19405977, 77221987.25246386, 78949864.89139993, 80430219.37467334, 81657230.09289972, 82626072.42564295, 83332936.68700519, 83775043.08154033, 83950652.61076826], [-6611.882279032151, 4130432.444510821, 8251216.556393094, 12339515.821930718, 16379233.621290537, 20354464.71877934, 24249557.879891377, 28049177.486397732, 31738363.906895082, 35302592.3850398, 38727830.21372627, 42000591.970186, 45107992.59450139, 48037798.102466665, 50778473.73330928, 53319229.34297612, 55650061.863912486, 57761794.66458175, 59646113.65366134, 61295599.987053506, 62703759.24892253, 63865046.99204711, 64774890.53711261, 65429706.94546075, 65826917.09463785, 65964955.80137688, 65843277.95306443, 65462360.623045154, 64823701.16251949, 63929811.276416935, 62784207.107283704, 61391395.36658522, 59756855.56813299, 57887018.43457492, 55789240.56214178, 53471775.44371209, 50943740.96480785, 48215083.50086351, 45296538.757360175, 42199589.50731273, 38936420.392457515, 35519869.966223486, 31963380.16711487, 28280943.42140297, 24487047.582969874, 20596618.9267251, 16624963.419443188, 12587706.498898877, 8500731.597714005, 4380117.653484329, 242075.85064780738, -3897114.1572345286, -8021168.294085015, -12113862.124863595, -16159094.684021544, -20140951.81751427, -24043768.789039683, -27852191.903781362, -31551238.906983748, -35126357.91915543, -38563484.67595962, -41849097.84714655, -44970272.21657197, -47914729.51407265, -50670886.69912171, -53227901.50628562, -55575715.07356469, -57705091.48616417, -59607654.080216184, -61275918.36417214, -62703321.42897991, -63884247.731335014, -64814051.149627164, -65489073.22605843, -65906657.52349486, -66065160.041331746, -65963955.64986992, -65603440.517857365, -64985030.52478979, -64111155.66372397, -62985250.45767632, -61611740.42707066, -59996024.66185453, -58144454.567054726, -56064308.86547591, -53763764.95612285, -51251866.74097227, -48538489.04653106, -45634298.78058191, -42550712.97673023, -39299853.89137809, -35894501.33041904, -32348042.392705422, -28674418.827772267, -24888072.215228103, -21003887.181031663, -17037132.873950493, -13003402.932340851, -8918554.177355412, -4798644.273365416, -659868.6008334861]), (600, 132636221.70894425, 7.462883194947908e-05, [132453020.58875743, 132244585.39936492, 131622372.31958774, 130588836.44883879, 129148056.64565292, 127305719.40003337, 125069096.35994917, 122447015.60023502, 119449826.748104, 116089360.10251954, 112378879.9092498, 108333031.97675101, 103967785.83908358, 99300371.69404301, 94349212.36704087, 89133850.56809539, 83674871.73012067, 77993822.7331356, 72113126.83522381, 66055995.14682812, 59846334.99742672, 53508655.55702644, 47067971.08510044, 40549702.189364605, 33979575.48412194, 27383522.044535916, 20787575.05829755, 14217767.078540409, 7700027.283840084, 1260079.1511394829, -5076661.054553618, -11285184.572885042, -17340988.827457264, -23220174.10235751, -28899537.83456485, -34356666.15027994, -39570022.28459459, -44519031.53567361, -49184162.419131145, -53547003.70315927, -57590337.02079131, -61298204.77393551, -64655973.06197552, -67650389.38789226, -70269634.9141384, -72503371.06362255, -74342780.28196289, -75780600.80042511, -76811155.26338878, -77430373.10632259, -77635806.59578422, -77426640.46874432, -76803695.13192472, -75769423.40790941, -74327900.84063756, -72484809.59792013, -70247416.03377481, -67624541.99967356, -64626530.01681275, -61265202.447533086, -57553814.82741252, -53507003.54235163, -49140728.05783967, -44472207.92912291, -39519854.8414374, -34303199.949465156, -28842816.803120088, -23160240.16493567, -17277881.039502632, -11218938.25101197, -5007306.918033184, 1332515.8132460224, 7775527.405008559, 14296318.551896196, 20869173.404149435, 27468170.994537145, 34067287.46936865, 40640498.72029567, 47161883.01238051, 53605723.204388164, 59946608.15803819, 66159532.93655553, 72219997.39817315, 78104102.79544011, 83788646.00040339, 89251210.98395011, 94470257.18915103, 99425204.4497646, 104096514.12016529, 108465766.09677209, 112515731.4271172, 116230440.22097312, 119595244.59594661, 122596876.40915696, 125223499.54787439, 127464756.57269509, 129311809.53053834, 130757374.77595355, 131795751.66309263, 132422844.99639162, 132636181.15159096], [-4513.310923675655, 6361865.509513918, 12703095.454435958, 18994142.975689966, 25210172.766766995, 31326645.801193826, 37319416.198999375, 43164826.53899717, 48839801.24029557, 54321937.644615345, 59589594.43976027, 64621977.0755478, 69399219.8344447, 73902464.23331505, 78113933.4466849, 82017002.45819016, 85596263.662638, 88837587.66044201, 91728179.00397652, 94256626.67610559, 96412949.10236302, 98188633.51782066, 99576669.53526358, 100571576.78077249, 101169426.4887385, 101367856.97140233, 101166082.90157276, 100564898.37311447, 99566673.72687653, 98175346.15537435, 96396404.12440012, 94236865.67220348, 91705250.67391737, 88811547.18091401, 85567171.96848333, 81984925.44859506, 78078941.12533379, 73864629.79448801, 69358618.70758057, 64578685.93990708, 59543690.22285554, 54273496.51802113, 48788897.62525923, 43111532.13480719, 37263799.04728323, 31268769.396917105, 25150095.226803977, 18931916.27427119, 12638764.734070556, 6295468.474448808, -72946.91295548121, -6441356.8431510255, -12784636.858627934, -19077761.75009478, -25295904.28593196, -31414533.160473406, -37409509.77432332, -43257183.4642865, -48934484.806412645, -54419016.62360904, -59689142.33789681, -64724071.31813072, -69503940.88627136, -74009894.65826318, -78224156.91045669, -82130102.67829259, -85712323.31003341, -88956687.21849357, -91850395.59152542, -94382032.84121561, -96541611.59438884, -98320612.04668787, -99712015.52644612, -100710332.13639724, -101311622.36481133, -101513512.58193927, -101315204.36061358, -100717477.58493711, -99722687.33584899, -98334754.56552884, -96559150.59834366, -94402875.51955153, -91874430.53754334, -88983784.42821877, -85742334.19468139, -82162860.09746253, -78259475.2326069, -74047569.8563357, -69543750.67636181, -64765775.34908061, -59732482.44064885, -54463717.12873713, -48980252.93798114, -43303709.8170532, -37456468.88092519, -31461584.15429689, -25342691.66435101, -19123916.24122346, -12829776.393418068, -6485087.63339515, -114864.63459793178]), (900, 194217488.09202757, 6.067137464739563e-05, [193989796.81056884, 193689559.48818207, 192794717.79239967, 191308801.869625, 189237674.5135265, 186589507.9970176, 183374751.7914632, 179606091.2994766, 175298397.76315534, 170468669.5477752, 165135965.03129458, 159321327.36577722, 153047701.4059593, 146339843.1366721, 139224221.95249686, 131728916.17810602, 123883502.24066506, 115718937.93401341, 107267440.23253825, 98562358.13933325, 89638041.07029215, 80529703.2932706, 71273284.95833284, 61905310.26792182, 52462743.34595975, 42982842.376137525, 33503012.584790155, 24060658.64906133, 14693037.11286572, 5437109.393368079, -3670604.041426469, -12594167.748845354, -21298373.07281313, -29748877.030781828, -37912337.778612405, -45756546.11918689, -53250552.53638557, -60364789.25374724, -67071186.83758633, -73343284.8846832, -79156336.36031766, -84487405.17448455, -89315456.61415243, -93621440.27529037, -97388365.16881213, -100601366.70554036, -103247765.29599792, -105317116.3352238, -106801251.37493606, -107694310.32113217, -107992764.52905507, -107695430.70445253, -106803475.55435857, -105320411.16860932, -103252081.15003158, -100606637.54517792, -97394508.66812989, -93628357.94257645, -89323033.92435212, -84495511.70177451, -79164825.90610535, -73351995.59641767, -67079941.31652565, -60373394.65234576, -53258800.64649599, -45764213.45837055, -37919185.680185, -29754651.74821503, -21302805.909375843, -12596975.224977659, -3671488.1133941547, 5438461.049553374, 14696950.7382911, 24067473.945816375, 33513082.153864246, 42996531.03756338, 52480427.33014865, 61927376.26873888, 71300129.04022269, 80561729.64650494, 89675660.61015296, 98605986.94636177, 107317497.83388, 115775845.42743333, 123947680.26347873, 131800782.72674184, 139304190.0596619, 146428318.41331586, 153145079.4604182, 159427991.10956213, 165252281.88517296, 170594988.5613122, 175435046.66691288, 179753373.502249, 183532943.34179297, 186758854.52716365, 189418388.186189, 191501058.34510195, 192998653.23728725, 193905267.64713606, 194217326.1601825], [-6607.147787926816, 9079347.1408798, 18129380.564889174, 27107776.915897626, 35979102.86037062, 44708347.772493556, 53261061.896335185, 61603492.29181973, 69702716.02820404, 77526770.09885132, 85044777.54559684, 92227069.29412287, 99045301.21954162, 105472565.98098895, 111483499.1829858, 117054379.44514199, 122163221.98540305, 126789865.34737563, 130916050.92955653, 134525495.00267267, 137603952.93101895, 140139275.3454592, 142121456.0448264, 143542671.43735126, 144397311.36896953, 144682001.21436948, 144395615.14533734, 143539280.5245827, 142116373.40745398, 140132505.1709839, 137595500.32203233, 134515365.5759112, 130904250.32518353, 126776398.6560513, 122148093.10536791, 117037590.37705307, 111465049.2755355, 105452451.14044067, 99023513.09661461, 92203594.46385197, 85019596.69425674, 77499857.23418891, 69674037.72747499, 61573007.00254817, 53228719.30366689, 44674088.2453722, 35942856.98738509, 27069465.140457258, 18088912.92677496, 9036623.130342476, -51698.61989861533, -9140204.683666881, -18193046.804449946, -27174517.50877965, -36049190.951742075, -44782062.65261152, -53338687.56853771, -61685315.96098118, -69789026.51780616, -77617856.20441243, -85140926.33132412, -92328564.33924492, -99152420.82067095, -105585581.31784236, -111602672.45211911, -117179961.97006108, -122295452.30996579, -126928967.32053468, -131062231.79104783, -134678943.47905433, -137764837.35370958, -140307741.80199754, -142297626.57649645, -143726642.29695287, -144589151.35186934, -144881750.0777366, -144603282.12892342, -143754842.98717755, -142339775.59460768, -140363657.12447128, -137834276.94681036, -134761605.87255666, -131157756.79876068, -127036936.91162232, -122415391.634531, -117311340.54318897, -111744905.49987023, -105738031.2907031, -99314399.07957476, -92499333.01804546, -85319700.38131319, -77803805.62356417, -69981278.76999587, -61882958.58531118, -53540770.98059766, -44987603.13640162, -36257173.83905156, -27383900.541045547, -18402763.670205228, -9349168.721893905, -258806.6781576846])]:
plt.plot(x, y, label=f'{year}, {r_max / 47279994.64428885:.4g}, {u_min / 0.00012962300739400188:.3g}')
plt.title("""Retarded interaction - Sun-Mercury
Changes of orbit after many orbital period
u[0] = 1.296e-04 r%[0] = 4.728e+07""")
plt.legend(title='period, r_max / r[0], u_min / u[0]')
plt.xlabel('r_x%')
plt.ylabel('r_y%')
plt.show()