变步长梯形求解微分方程

最新推荐文章于 2022-06-09 10:43:18 发布

DragonLord0323

最新推荐文章于 2022-06-09 10:43:18 发布

阅读量2.8k

点赞数

分类专栏：计算方法文章标签：出版算法教育

本文链接：https://blog.csdn.net/DragonLord0323/article/details/394234

版权

该博客介绍了一种使用变步长梯形公式进行微分方程数值求解的方法，通过调整步长以提高精度。代码示例展示了如何实现这一算法，并从《数值分析简明教程》中引用了相关算法流程图。程序中包含误差控制，以确保达到预期的精度。

摘要由CSDN通过智能技术生成

/**
***变步长梯形求解微分方程***

获取初值： T1 = h/2［ f(x(k)) + f(x(k+1)) ］

                                         n-1
变步长梯形公式： T2n = 1/2*Tn + h/2 * ∑ f ( x(k+1/2) )
                                         k=0
步长： h=b-a/n

属性：数值积分法

《数值分析简明教程》-2 Editon -高等教育出版社- page 67 算法流程图

代码维护：2005.6.14 DragonLord
**/
#include<iostream.h>
#include<math.h>
#include<stdio.h>

double f(double x)
{
    double f;
if(x==0)f=1;      //                            1
else f=sin(x)/x; //   举例方程 I = ∫ sin(x)/x dx
return f;             //                              0
}