微积分是许多学科的基础,在编程方面也有很大的作用,程序里写积分很多时候是用梯形法,今天学了一种叫辛普森积分的方法,感觉很好用只不过辛普森积分法在区间较大时误差比较大,所以一般用二分减小区间大小以获得更加精确的结果。现在以求两个圆柱体相交部分的体积为例来学习辛普森积分:给出的是两个圆柱的半径,已知两个圆柱体垂直相交。
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#define eps 1e-9
using namespace std;
int r1,r2;
double f(double x){
return 8.0*(sqrt(r1*r1-x*x))*(sqrt(r2*r2-x*x));
}
double simpson(double x,double y){
return (y-x)*(f(x)+f(y)+4*f((x+y)/2))/6.0;
}
double jifen(double x,double y){
double ans=simpson(x,y);//二分区间积分减小误差
double mid=(x+y)/2.0;
double left=simpson(x,mid);
double right=simpson(mid,y);
if(fabs(ans-(left+right))<eps)return ans;
else{
return jifen(x,mid)+jifen(mid,y);
}
}
int main()
{
while(~scanf("%d%d",&r1,&r2)){
if(r1<r2) swap(r1,r2);
double x=1.0*r1;
double y=1.0*r2;
printf("%.4lf\n",jifen(0.0,y));
}
return 0;
}