自适应Simpson积分(近似求积分)
2017-08-08 02:24
246 查看
来自 http://blog.csdn.net/xiaosshhaa/article/details/70799670 //1261: 地狱飞龙 #include <bits/stdc++.h> using namespace std; const double eps=1e-7; double v1,v2,x,k; int t; //积分函数 double F(double t) { return k/((x-v2*t)*(x-v2*t)+(v1*t)*(v1*t)); } //自适应Simpson积分 double simpson(double a,double b) { double c=a+(b-a)/2; return (F(a)+4*F(c)+F(b))*(b-a)/6; } double asr(double a,double b,double eps,double A) { double c=a+(b-a)/2; double L=simpson(a,c),R=simpson(c,b); if(fabs(L+R-A)<=15*eps) return L+R+(L+R-A)/15.0; return asr(a,c,eps/2,L)+asr(c,b,eps/2,R); } double asr(double a,double b,double eps) { return asr(a,b,eps,simpson(a,b)); } int main() { scanf("%d",&t); while(t--) { scanf("%lf%lf%lf%lf",&v1,&v2,&x,&k); printf("%.2f\n", asr(0,10000000,eps)); } return 0; }
相关文章推荐
- 自适应Simpson积分
- 洛谷.4525.[模板]自适应辛普森法1(Simpson积分)
- 【bzoj1502】[NOI2005]月下柠檬树 自适应Simpson积分
- 学习笔记: 自适应Simpson积分
- 自适应simpson积分
- 自适应辛普森(近似求积分)
- 自适应辛普森(近似求积分模板)
- acdream1234(Simpson积分)
- Laplace近似积分
- BZOJ 1502([NOI2005]月下柠檬树-Simpson积分)
- hdu 1724 Ellipse(自适应Simpson积分) (模板)
- hdu 4498 自适应simpson
- bzoj 2178 圆的面积并【simpson积分】
- HDU - 4498 Function Curve 自适应simpson
- LA 3485 自适应Simpson积分 解题报告
- hdu 4498 自适应simpson
- 数值积分之Simpson公式与梯形公式
- bzoj 1502 月下柠檬树【Simpson积分】
- HDOJ 1724 Ellipse 自适应Simpson
- Simpson 积分算法