自适应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);
}
全部评论

相关推荐

不愿透露姓名的神秘牛友
07-02 17:58
点赞 评论 收藏
分享
代码飞升:别用口语,后端就写后端,前端就写前端,最后别光后悔
点赞 评论 收藏
分享
练习生懒羊羊:开飞机把这个公司创飞吧
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

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