20xx302026王勛(編輯修改稿)
2024-08-31 07:34 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 }double f2( double x) //方程組（2）式{ return g/x。}double f3( double x, double y) //方程組（3）式{ return x*y。}double f4(double x) //方程組（4）式{ return *x。}double RungeKutta(double x0, double y0, double z0,double m0, double a, double b, int n,double x[], double y[], double z[],double m[], double (*f1)(double, double, double),double (*f2)(double), double (*f3)(double, double), double (*f4)(double )) //RungeKutta法{double k1, k2, k3, k4 ,m1, m2, m3, m4, n1, n2, n3, n4, c1, c2, c3, c4 ,h,t。int count。count =0。h =(ba)/n。x[0] = x0。 //初值y[0] = y0。z[0] = z0。m[0] = m0。FILE *file= fopen(“”,”w”)。for(int i = 1。i=n。i++) {t = a+(i1)*h。k1 = f1(x[i1], y[i1], z[i1])。m1 = f2(x[i1])。n1 = f3(x[i1], y[i1])。c1 = f4(x[i1])。k2 = f1(x[i1]+k1*h/, y[i1]+m1*h/, z[i1]+n1*h/)。m2 = f2(x[i1]+k1*h/)。n2 = f3(x
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1