采用系数矩阵 通用化实现 1到4阶“龙格库塔法”.zip

轨道力学第3版《Orbital Mechanics for Engineering Students, Third Edition.pdf》改写的 1到4阶龙格库塔法： 1阶龙格库塔法： RK1 RK1 步长= 0.010 t = 0.000 x(t) = 0.00000000 Vx(t) = 0.00000000 RK1 步长= 0.010 t = 0.010 x(t) = 0.00000000 Vx(t) = 0.00000000 RK1 步长= 0.010 t = 0.020 x(t) = 0.00000000 Vx(t) = 0.00004000 RK1 步长= 0.010 t = 0.030 x(t) = 0.00000040 Vx(t) = 0.00011998 RK1 步长= 0.010 t = 0.040 x(t) = 0.00000160 Vx(t) = 0.00023990 RK1 步长= 0.010 t = 0.050 x(t) = 0.00000400 Vx(t) = 0.00039973 RK1 步长= 0.010 t = 0.060 x(t) = 0.00000800 Vx(t) = 0.00059944 RK1 步长= 0.010 t = 0.070 x(t) = 0.00001399 Vx(t) = 0.00083897 RK1 步长= 0.010 t = 0.080 x(t) = 0.00002238 Vx(t) = 0.00111829 RK1 步长= 0.010 t = 0.090 x(t) = 0.00003356 Vx(t) = 0.00143734 RK1 步长= 0.010 t = 0.100 x(t) = 0.00004794 Vx(t) = 0.00179607 RK1 步长= 0.010 t = 0.110 x(t) = 0.00006590 Vx(t) = 0.00219440 RK1 步长= 0.010 t = 0.120 x(t) = 0.00008784 Vx(t) = 0.00263229 RK1 步长= 0.010 t = 0.130 x(t) = 0.00011416 Vx(t) = 0.00310965 RK1 步长= 0.010 t = 0.140 x(t) = 0.00014526 Vx(t) = 0.00362640 RK1 步长= 0.010 t = 0.150 x(t) = 0.00018152 Vx(t) = 0.00418248 RK1 步长= 0.010 t = 0.160 x(t) = 0.00022335 Vx(t) = 0.00477780 RK1 步长= 0.010 t = 0.170 x(t) = 0.00027113 Vx(t) = 0.00541226 RK1 步长= 0.010 t = 0.180 x(t) = 0.00032525 Vx(t) = 0.00608578 RK1 步长= 0.010 t = 0.190 x(t) = 0.00038611 Vx(t) = 0.00679825 RK1 步长= 0.010 t = 0.200 x(t) = 0.00045409 Vx(t) = 0.00754958 RK1 步长= 0.010 t = 0.210 x(t) = 0.00052959 Vx(t) = 0.00833966 RK1 步长= 0.010 t = 0.220 x(t) = 0.00061298 Vx(t) = 0.00916837 RK1 步长= 0.010 t = 0.230 x(t) = 0.00070467 Vx(t) = 0.01003560 RK1 步长= 0.010 t = 0.240 x(t) = 0.00080502 Vx(t) = 0.01094124 RK1 步长= 0.010 t = 0.250 x(t) = 0.00091443 Vx(t) = 0.01188515 RK1 步长= 0.010 t = 0.260 x(t) = 0.00103329 Vx(t) = 0.01286721 RK1 步长= 0.010 t = 0.270 x(t) = 0.00116196 Vx(t) = 0.01388728 RK1 步长= 0.010 t = 0.280 x(t) = 0.00130083 Vx(t) = 0.01494523 RK1 步长= 0.010 t = 0.290 x(t) = 0.00145028 Vx(t) = 0.01604092 RK1 步长= 0.010 t = 0.300 x(t) = 0.00161069 Vx(t) = 0.01717419 RK1 步长= 0.010 t = 0.310 x(t) = 0.00178243 Vx(t) = 0.01834490 RK1 步长= 0.010 t = 0.320 x(t) = 0.00196588 Vx(t) = 0.01955289 RK1 步长= 0.010 t = 0.330 x(t) = 0.00216141 Vx(t) = 0.02079801 RK1 步长= 0.010 t = 0.340 x(t) = 0.00236939 Vx(t) = 0.02208009 RK1 步长= 0.010 t = 0.350 x(t) = 0.00259019 Vx(t) = 0.02339896 RK1 步长= 0.010 t = 0.360 x(t) = 0.00282418 Vx(t) = 0.02475445 RK1 步长= 0.010 t = 0.370 x(t) = 0.00307173 Vx(t) = 0.02614638 RK1 步长= 0.010 t = 0.380 x(t) = 0.00333319 Vx(t) = 0.02757458 RK1 步长= 0.010 t = 0.390 x(t) = 0.00360894 Vx(t) = 0.02903886 RK1 步长= 0.010 t = 0.400 x(t) = 0.00389933 Vx(t) = 0.03053902 RK1 步长= 0.010 t = 0.410 x(t) = 0.00420472 Vx(t) = 0.03207489 RK1 步长= 0.010 t = 0.420 x(t) = 0.00452546 Vx(t) = 0.03364625 RK1 步长= 0.010 t = 0.430 x(t) = 0.00486193 Vx(t) = 0.03525292 RK1 步长= 0.010 t = 0.440 x(t) = 0.00521446 Vx(t) = 0.03689468 RK1 步长= 0.010 t = 0.450 x(t) = 0.00558340 Vx(t) = 0.03857133 RK1 步长= 0.010 t = 0.460 x(t) = 0.00596912 Vx(t) = 0.04028265 RK1 步长= 0.010 t = 0.470 x(t) = 0.00637194 Vx(t) = 0.04202842 RK1 步长= 0.010 t = 0.480 x(t) = 0.00679223 Vx(t) = 0.04380843 RK1 步长= 0.010 t = 0.490 x(t) = 0.00723031 Vx(t) = 0.04562245 RK1 步长= 0.010 t = 0.500 x(t) = 0.00768654 Vx(t) = 0.04747025 RK1 步长= 0.010 t = 0.510 x(t) = 0.00816124 Vx(t) = 0.04935159 RK1 步长= 0.010 t = 0.520 x(t) = 0.00865475 Vx(t) = 0.05126625 RK1 步长= 0.010 t = 0.530 x(t) = 0.00916742 Vx(t) = 0.05321398 RK1 步长= 0.010 t = 0.540 x(t) = 0.00969956 Vx(t) = 0.05519453 RK1 步长= 0.010 t = 0.550 x(t) = 0.01025150 Vx(t) = 0.05720766 RK1 步长= 0.010 t = 0.560 x(t) = 0.01082358 Vx(t) = 0.05925312 RK1 步长= 0.010 t = 0.570 x(t) = 0.01141611 Vx(t) = 0.06133064 RK1 步长= 0.010 t = 0.580 x(t) = 0.01202942 Vx(t) = 0.06343998 RK1 步长= 0.010 t = 0.590 x(t) = 0.01266382 Vx(t) = 0.06558087 RK1 步长= 0.010 t = 0.600 x(t) = 0.01331962 Vx(t) = 0.06775303 RK1 步长= 0.010 t = 0.610 x(t) = 0.01399715 Vx(t) = 0.06995621 RK1 步长= 0.010 t = 0.620 x(t) = 0.01469672 Vx(t) = 0.07219013 RK1 步长= 0.010 t = 0.630 x(t) = 0.01541862 Vx(t) = 0.07445450 RK1 步长= 0.010 t = 0.640 x(t) = 0.01616316 Vx(t) = 0.07674906 RK1 步长= 0.010 t = 0.650 x(t) = 0.01693065 Vx(t) = 0.07907351 RK1 步长= 0.010 t = 0.660 x(t) = 0.01772139 Vx(t) = 0.08142756 RK1 步长= 0.010 t = 0.670 x(t) = 0.01853566 Vx(t) = 0.08381093 RK1 步长= 0.010 t = 0.680 x(t) = 0.01937377 Vx(t) = 0.08622332 RK1 步长= 0.010 t = 0.690 x(t) = 0.02023601 Vx(t) = 0.08866443 RK1 步长= 0.010 t = 0.700 x(t) = 0.02112265 Vx(t) = 0.09113397 RK1 步长= 0.010 t = 0.710 x(t) = 0.02203399 Vx(t) = 0.09363162 RK1 步长= 0.010 t = 0.720 x(t) = 0.02297031 Vx(t) = 0.09615707 RK1 步长= 0.010 t = 0.730 x(t) = 0.02393188 Vx(t) = 0.09871003 RK1 步长= 0.010 t = 0.740 x(t) = 0.02491898 Vx(t) = 0.10129017 RK1 步长= 0.010 t = 0.750 x(t) = 0.02593188 Vx(t) = 0.10389717 RK1 步长= 0.010 t = 0.760 x(t) = 0.02697085 Vx(t) = 0.10653071 RK1 步长= 0.010 t = 0.770 x(t) = 0.02803616 Vx(t) = 0.10919048 RK1 步长= 0.010 t = 0.780 x(t) = 0.02912806 Vx(t) = 0.11187613 RK1 步长= 0.010 t = 0.790 x(t) = 0.03024682 Vx(t) = 0.11458735 RK1 步长= 0.010 t = 0.800 x(t) = 0.03139270 Vx(t) = 0.11732381 RK1 步长= 0.010 t = 0.810 x(t) = 0.03256594 Vx(t) = 0.12008515 RK1 步长= 0.010 t = 0.820 x(t) = 0.03376679 Vx(t) = 0.12287105 RK1 步长= 0.010 t = 0.830 x(t) = 0.03499550 Vx(t) = 0.12568116 RK1 步长= 0.010 t = 0.840 x(t) = 0.03625231 Vx(t) = 0.12851514 RK1 步长= 0.010 t = 0.850 x(t) = 0.03753746 Vx(t) = 0.13137264 RK1 步长= 0.010 t = 0.860 x(t) = 0.03885119 Vx(t) = 0.13425332 RK1 步长= 0.010 t = 0.870 x(t) = 0.04019372 Vx(t) = 0.13715681 RK1 步长= 0.010 t = 0.880 x(t) = 0.04156529 Vx(t) = 0.14008276 RK1 步长= 0.010 t = 0.890 x(t) = 0.04296612 Vx(t) = 0.14303082 RK1 步长= 0.010 t = 0.900 x(t) = 0.04439642 Vx(t) = 0.14600062 RK1 步长= 0.010 t = 0.910 x(t) = 0.04585643 Vx(t) = 0.14899179 RK1 步长= 0.010 t = 0.920 x(t) = 0.04734635 Vx(t) = 0.15200398 RK1 步长= 0.010 t = 0.930 x(t) = 0.04886639 Vx(t) = 0.15503682 RK1 步长= 0.010 t = 0.940 x(t) = 0.05041676 Vx(t)