自适应simpton积分模板



double f(double x){
	return b*sqrt(1-(x*x)/(a*a)); 
} 
double simpton(double a,double b){
	double c=(a)+(b-a)/2;
	return (f(a)+4*f(c)+f(b))*(b-a)/6;
}

double integral(double L,double R,double Eps){
	double mid=L+(R-L)/2;
	double ST=simpson(L,R),SL=simpson(L,mid),SR=simpson(mid,R);
	if (fabs(SL+SR-ST)<=15*Eps) return SL+SR+(SL+SR-ST)/15;
	
	return integral(L,mid,Eps/2)+integral(mid,R,Eps/2);
}
全部评论

相关推荐

这不纯纯作弊了吗😢😢😢
编程界菜鸡:信这个的这辈子有了,这智商你靠啥都没用
你找工作的时候用AI吗?
点赞 评论 收藏
分享
不愿透露姓名的神秘牛友
06-11 13:34
offe从四面八方来:我真的没时间陪你闹了
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务