【求助】用龍格庫塔法數值計算（畢業論文相關，急。導師不管事只能網上求助了）

問題如下： 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) %參數表順序依次是微分方程組的函數名稱，初始值向量

>>阅读原文<<