MATLAB四阶龙格库塔算法

用MATLAB编写的四阶龙格库塔算法，可以直接调用状态微分方程，但是需要满足其格式，可以修改算法的步长

function dx= apolloeq（tx）
% Apollo卫星的轨迹运动方程，参考《高等应用数学问题的MATLAB求解（第三版）》薛定宇
% p251
mu = 1/82.45;
mu1 = 1-mu;
r1 = sqrt（（x（1）+mu）^2+x（3）^2）;
r2 = sqrt（（x（1）-mu1）^2+x（3）^2）;
dx = [x（2）;
    2*x（4）+x（1）-mu1*（x（1）+mu）/r1^3-mu*（x（1）-mu1）/r2^3;
    x（4）;
    -2*x（2）+x（3）-mu1*x（3）/r1^3-mu*x（3）/r2^3];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        338  2017-11-11 10:57  RK_4\apolloeq.m

     文件        773  2017-11-11 20:33  RK_4\example7_14.m

     文件        570  2017-11-11 11:23  RK_4\rk_4.m

     目录          0  2017-11-22 22:00  RK_4

----------- ---------  ---------- -----  ----

                 1681                    4