辛普森积分法

/*

  • 辛普森积分公式:sum(l,r)=( f(l)+f(r)+4*f((l+r)/2) ) * (r-l) / 6
  • 解: 对 f 函数在 (l,r) 的积分可以用 以上函数拟合
  • 欲对 (l,r) 的一段函数积分并给出了精度要求
  • 可以对此段区间,递归使用上述公式(以设置的精度作为递归终点)即可以得到满足精度要求的答案。

*/

include<bits/stdc++.h>

using namespace std;

//读入优化
void in(int &x){
x=0;char c=getchar();
int y=1;
while(c<'0'||c>'9'){if(c=='-')y=-1;c=getchar();}
while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}
x*=y;
}
void o(int x){
if(x<0){x=-x;putchar('-');}
if(x>9)o(x/10);
putchar(x%10+'0');
}

//设置需要控制的精度
const double eps=1e-7;

double a;//一个变量

//需要积分的函数 例: a^(a/x-x)
double F(double x){
return pow(x,a/x-x);
}

double simpson(double l,double r){ //辛普森拟合积分计算函数
return ( F(l)+F(r)+4.0*F((l+r)/2.0) ) * (r-l) / 6;
}

double simpson(double l, double r, double exps, double val) {//递归积分函数
double mid = (l + r) / 2; // 中点
double lval = simpson(l, mid), rval = simpson(mid, r);//lval: 左半部分的辛普森积分 , rval 右半部分的辛普森积分
if (fabs(lval + rval - val) <= 15 * exps)//val:本段区间的积分 if(精度满足要求)
return lval + rval + (lval + rval - val) / 15;//返回答案
return simpson(l, mid, exps / 2, lval) + simpson(mid, r, exps / 2, rval);
}

double simpson(double l, double r, double exps) {
return simpson(l, r, exps, simpson(l, r));//调用
}

signed main() {
scanf("%lf", &a);//
printf("%.5lf", simpson(1e-11, 30.0, eps));//eps初始化为-e7;
return 0;
}

原文地址:https://www.cnblogs.com/yesuweiYYYY/p/14313777.html