自适应辛普森公式求积分模板

最新推荐文章于 2023-04-23 09:08:31 发布

岁月匆匆·

最新推荐文章于 2023-04-23 09:08:31 发布

阅读量954

点赞数 1

分类专栏：模板匹配

本文链接：https://blog.csdn.net/vagebird/article/details/70544179

版权

模板匹配专栏收录该内容

2 篇文章 0 订阅

订阅专栏

话不多说直接上代码：

#include <bits/stdc++.h>

using namespace std;

double v1, v2, x, k;

double F(double t)//定义一个自己的公式F(x)
{
   return k/((x-v2*t)*(x-v2*t)+v1*v1*t*t);
}
double simpson(double a, double b)
{
    double c = a + (b-a)/2.0;
    return (F(a) + 4*F(c) + F(b)) * (b-a)/6.0;
}
double asr(double a, double b, double eps, double A)
{
    double c = a + (b-a)/2.0;
    double L = simpson(a, c), R = simpson(c,b);
    if(fabs(L + R - A) <= 15*eps)  return  L + R + (L + R - A)/15.0;
    return asr(a, c, eps/2.0, L) + asr(c, b, eps/2.0, R);
}
double ac(double a, double b, double eps)//a,b表示范围,eps表示精度
{
    return asr(a, b, eps, simpson(a, b));
}
int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        scanf("%lf %lf %lf %lf", &v1, &v2, &x, &k);
        printf("%.2f\n", ac(0, 10000000, 0.0000001));
    }
    return 0;
}