自适应辛普森模版

最新推荐文章于 2022-09-02 11:45:59 发布

黄二十九

最新推荐文章于 2022-09-02 11:45:59 发布

阅读量226

点赞数

分类专栏：算法刷题文章标签：数学

本文链接：https://blog.youkuaiyun.com/qq_30798083/article/details/121402800

版权

算法刷题专栏收录该内容

26 篇文章

订阅专栏

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
double a,b,c,d,l,r;
const double eps = 1e-10;
double f(double x)
{
    return (c*x+d)/(a*x+b);//函数
}
double simpson(double ll,double rr)
{
    return  (rr-ll)/6.0*( f(ll)+f(rr)+4.0*f( (ll+rr)/2.0) );
}
double integral(double ll,double rr)
{
    double mid=(ll+rr)/2.0;
    double res=simpson(ll,rr);
    if(fabs( res-simpson(ll,mid)-simpson(mid,rr) ) <eps) return res;
    else return integral(ll,mid)+integral(mid,rr);
}
int main()
{
    scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&l,&r);
    printf("%.6lf",integral(l,r));
    return 0;
}

洛谷P4525题解

同解决 P4526（这个要证明a<0时发散）