【求助】用龙格库塔法数值计算（毕业论文相关，急。导师不管事只能网上求助了）

问题如下： x_1=x,x_2=x’, 相应初始条件为：x_1 (0)=1,x_2 (0)=0 求： 突加载荷时，也就是P（t）=p时 （1） t-x曲线； （2） 不同P/T/U下的x-x’曲线； （3） 不同U/T下P-Xmax曲线 程序如下： 主程序： function [x,y]=runge_kutta1(ufunc,y0,h2,a,b) %参数表顺序依次是微分方程组的函数名称，初始值向量

>>阅读原文<<