("下面为四阶龙格-库塔法结果:\n"); 24 double k1 = 0.0, k2 = 0.0, k3 = 0.0, k4 = 0.0; 25 for (int i = 1; i <= n; i++) 26 { 27 x[i] = 0.0, y[i] = 0.0; 28 } 29 x[0] = a; 30 y[0] = m; 31 for (int i = 1; i <= n; i++) 32 { 33 k1 = h...
那么利用标准4阶龙格-库塔法: xn+1=xn+1/6(k1+2k2+2k3+k4) yn+1=yn+1/6(l1+2l2+2l3+l4) zn+1=zn+1/6(m1+2m2+2m3+m4) 则其中y的误差函数为:△y=yn+1-yn=1/6(l1+2l2+2l3+l4) 其中:k1=△t*f(tn,xn,yn,zn) k2=△t*f(tn+△t/2,xn+k1/2,yn+l1/2,zn+m1/2) ...
c求解四阶龙格库塔法某一地区的病菌传染,三种人员的人数的状态方程,即可能受传染的人数x1,已被传染的病的人数x2,及已治愈的人数x3,并假设他们是因接触而传染。 设:α是x1中单位时间内的传染系数,β是x2中单位时间内被治愈的比例系数,则得以下状态方程: f`(x1)=-α*x1*x2; f`(x2)=α*x1*x2-β*x...
四阶龙格库塔法原理C代码差分方法数值分析简明教程2editon高等教育出版社page105算法流程图代码维护 /** ***四阶Runge-Kutta法*** 经典格式: y(n+1) = y(n) + h/6 ( K1 + 2*K2 + 2*K3 + K4 ) K1 = f( x(n) , y(n) ) K2 = f( x(n+1/2) , y(n) + h/2*K1 ) K3 = f(...
内容提示: 6 第 l期 焦作大学学报 JOURNAL OF JIAOZUO UNIVERSITY № 1 M ar 2001 四阶龙格一库塔算法的 C语言实现 毋玉芝 焦作财会 学校 ) 摘要 本文叙 述 了四阶龙 格一 库塔 算法的 C语 言 实现 过程 、数据存 储度 其蛄 果 的曲 线显 示.并以具体 实例 说明了这一过程。 关键词 龙格 库 ...
/*四阶Runge-Kutta法* 经典格式:y(n+1) = y(n) + h/6 ( K1 + 2*K2 + 2*K3 + K4 )K1 = f( x(n) , y(n) )K2 = f( x(n+1/2) , y(n) + h/2*K1 )K3 = f( x(n+1/2) , y(n) + h/2*K2 )K4 = f( x(n+1) , y(n) + h*K3 ) Runge-Kutta法是基于泰勒...
四阶龙格库塔法代码dragonlordkuttarunge /***四阶Runge-Kutta法***经典格式:y(n+1)=y(n)+h/6(K1+2*K2+2*K3+K4)K1=f(x(n),y(n))K2=f(x(n+1/2),y(n)+h/2*K1)K3=f(x(n+1/2),y(n)+h/2*K2)K4=f(x(n+1),y(n)+h*K3)Runge-Kutta法是基于泰勒展开方法,因而需要所求解具...
printf("用三阶龙格-库塔方法\n");RungeKutta(1,0,1,5,x,y,3,function);for(i=0;i<6;i++)printf("x[%d]=%f,y[%d]=%f\n",i,x[i],i,y[i]);printf("用四阶龙格-库塔方法\n");RungeKutta(1,0,1,5,x,y,4,function);for(i=0;i<6;i++)printf("x[%d]=%f,y[%d]=...
这就是四阶龙格库塔方程,据说这里的精度很好,可以精确到小数点后6位。 方程组中的h为差分的步长。 下面给出c语言实现的代码: #include<stdlib.h> #include<stdio.h> /*n表示几等分,n+1表示他输出的个数*/ int RungeKutta(double y0,double a,double b,int n,double *x,double *y,int style,double ...