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

相关推荐

舂锋:不能投什么岗都用一份简历,一般都是要看企业的岗位需求来写职业技能或者是项目经历,跟岗位相关的就写多一点。
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

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