四阶龙格库塔法原理C代码(共2页).doc

精选优质文档-倾情为你奉上/*四阶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法是基于泰勒展开方法，因而需要所求解具有较好的光滑性。 属性：差分方法 数值分析简明教程-2 Editon -高等教育出版社- page 105 算法流程图 代码维护：2005.6.14 DragonLord*/#include#include#include/*举例方程：y= y - 2*x / y ( 0x1 ) y(0) = 1*/double f(do