自适应simpson 积分

最新推荐文章于 2020-12-02 15:18:26 发布

斌同学斌同學

最新推荐文章于 2020-12-02 15:18:26 发布

阅读量607

点赞数 1

分类专栏：数论数学文章标签：数学数论

本文链接：https://blog.csdn.net/GreatJames/article/details/78152449

版权

该博客介绍了如何应用自适应Simpson积分法进行数值积分。通过示例代码展示了如何在C++中实现这个方法，从定义被积函数开始，到递归地细分区间直至满足给定的精度要求。最后，给出一个具体的积分实例，并输出了积分结果。

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

#include<bits/stdc++.h>
#include <ctime>
using namespace std;
typedef long long ll;

double f(double x)           ///被积函数
{
    return 0.5 * x * x;
}
double simpson(double a, double b)
{
    double c = a + (b - a) / 2;
    return (f(a) + 4 * f(c) + f(b)) * (b - a) / 6;
}
double asr(double a, double b, double epss, double A)
{
    double c = a + (b - a) / 2;
    double L = simpson(a, c) , R = simpson(c, b);
    if (fabs(L + R - A) <= 15 * epss)
    {
        return L + R + (L