用四阶龙格库塔法求解
Runge_Kutta_4.m—采用龙格-库塔四阶法计算ODE
用Fehlberg公式计算六阶RK ODE。
q1a2.m- %的作者
这个程序使用线性插值来计算在一定的t值下x的近似值。
l_intp (t,输入)——线性插值