rungekutta专题

龙格库塔法c语言,龙格——库塔(Rungekutta)法求解常微分方程

1.龙格库塔法的基本原理 该算法是构建在数学支持的基础之上的。对于一阶精度的拉格朗日中值定理有: 对于微分方程 y'=f(x,y) y(n+1)=y(n)+h*K1 K1=f(xn,yn) 当用点xn处的斜率近似值K1与右端点xn+1处的斜率K2的算术平均值作为平均斜率K*的近似值,那么就会得到二阶精度的改进拉格朗日中值定理: y(n+1)=y(n)+[h*( K1+ K2)/2] K1=f(x