C++利用复化Simpson公式构造的自适应算法计算积分

最新推荐文章于 2022-10-16 16:44:23 发布

神捉羊

最新推荐文章于 2022-10-16 16:44:23 发布

阅读量2.4k

点赞数 3

分类专栏：数值分析文章标签： Simpson公式数值分析复化梯形公式自适应

本文链接：https://blog.csdn.net/weixin_42373618/article/details/80557107

版权

本文介绍了使用C++编程实现复化Simpson公式来构造一个自适应算法，用于精确计算积分。通过引用具体的函数库和调用Simpson函数，程序能够根据需要自动调整细分区间，最终得出积分值0.4987111175。参考文献为黄云清和舒适的《数值计算方法》。

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

具体算法详见参考文献，本文仅仅给出C++利用复化Simpson公式构造的自适应算法计算积分的具体程序。

simpson.h:

#ifndef SIMPSON_H_INCLUDED
#define SIMPSON_H_INCLUDED
#include <iostream>
#include <math.h>
#include"fuction.h"
using namespace std;
double Simpson(double eps,double lower_bd,double upper_bd)
{
    double h=upper_bd-lower_bd,s,s1,s2;//h为步长
    int n=1,k ;
    s1=h/6*(fuction(lower_bd)+4*fuction((lower_bd+upper_bd)/2)+fuction(upper_bd));
    s=2*fuction(lower_bd+0.25*h)-fuction(lower_bd+0.25*h)+2*fuction(lower_bd+0.75*h);
    s2=0.5*s1+s*h/6;//初始化s、s1、s2
    do
    {
       h=0.5*h;
        n=2*n;
        s1=s2;
        s=0;//为方便求和
        for(k=0;k<n;k++)
        {
             s=s+2*fuction(lower_bd+(k+0.25)*h)-fuction(lower_bd+(k+0.5)*h)+2*fuction(lower_bd+(k+0.75)*h);
        }
        s2=0.5*s1+s*h/6;
    }w