数值计算课程设计(陕西科技大学).docVIP

龙贝格求积分算法 数值计算课程设计 PAGE \* MERGEFORMAT- 20 - PAGE \* MERGEFORMAT- 19 - 1、 经典四阶龙格库塔法解一阶微分方程 1.1、 算法说明 龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高，采取措施对误差进行抑制，所以其实现原理也较复杂。该算法是构建在数学支持的基础之上的。 4阶龙格-库塔方法(RK4)可模拟N=4的泰勒方法的精度。这种算法可以描述为，自初始点开始，利用下面的计算方法生成近似序列 (1-1) 1.2、 经典四阶龙格库塔法解一阶微分方程算法流程图 图1-1 经典四阶龙格库塔法解一阶微分方程算法流程图 图1-1 经典四阶龙格库塔法解一阶微分方程算法流程图 1.3、 经典四阶龙格库塔法解一阶微分方程程序调试 图1-2 经典四阶龙格库塔法解一阶微分方程程序调试 1.4、 经典四阶龙格库塔法解一阶微分方程代码 #include iostream #include iomanip using namespace std; //f为函数的入口地址，x0、y0为初值，xn为所求点，step为计算次数 double Runge_Kuta( double (*f)(double x, double y), double x0, double y0, double xn