【模板】自适应辛普森积分

最新推荐文章于 2025-06-06 11:44:45 发布

原创最新推荐文章于 2025-06-06 11:44:45 发布 · 265 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#c++ #c语言

收录于

模板同时被 2 个专栏收录

35 篇文章

订阅专栏

数学/数论

6 篇文章

订阅专栏

本文介绍了一种名为自适应辛普森积分的数值积分方法，通过C++代码实现原函数、三点辛普森公式和自适应策略，适用于求解连续函数的定积分。通过实例展示了如何在给定区间[L, R]上使用此方法，精度可达1e-6。

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
double a,b,c,d,L,R;
double f(double x)//原函数
{
	return (c*x+d)/(a*x+b);
}
double simpson(double l,double r)//三点辛普森
{
	double mid = (l+r)/2;
	return (r-l)*(f(l)+f(r)+4*f(mid))/6;
}
double integral(double l,double r,double eps)//自适应辛普森
{
	double mid = (l+r)/2;
	double SL = simpson(l,mid),SR = simpson(mid,r),ST = simpson(l,r);
	if(fabs(ST-SL-SR)<=(eps*15))
	    return SL+SR+(SL+SR-ST)/15;
	else return integral(l,mid,eps/2)+integral(mid,r,eps/2);
}
int main()
{
	cin>>a>>b>>c>>d>>L>>R;
	printf("%.6lf",integral(L,R,1e-6));
	return 0;
}