simpson公式求定积分（模板）

 #include<cstdio>
 #include<cmath>
 #include <algorithm>
 using namespace std;

 double r,R1;//公式需要用到的变量

 // simpson公式用到的函数，就是待积分函数
 double F(double x)
 {
     //return sqrt(r*r-x*x)*sqrt(R1*R1-x*x);
     return f(x);
 }

 // 三点simpson法。这里要求F是一个全局函数
 double simpson(double a, double b)
 {
     double c = a + (b-a)/;
     return (F(a)+*F(c)+F(b))*(b-a)/;
 }

 // 自适应Simpson公式（递归过程）。已知整个区间[a,b]上的三点simpson值A
 double asr(double a, double b, double eps, double A)
 {
     double c = a + (b-a)/;
     double L = simpson(a, c), R = simpson(c, b);
     if(fabs(L+R-A) <= *eps) return L+R+(L+R-A)/15.0;
     return asr(a, c, eps/, L) + asr(c, b, eps/, R);
 }

 // 自适应Simpson公式（主过程）
 double asr(double a, double b, double eps)
 {
     return asr(a, b, eps, simpson(a, b));
 }

 // 用自适应Simpson公式计算积分数值
 double getValue()
 {
     return asr(, r, 1e-)*; // 第一第二个参数为积分区间，第三个参数为精度
 }

 int main()
 {
     while(~scanf("%lf%lf",&r,&R1))
     {
         if(r>R1)
             swap(r,R1);
         printf("%.10f\n",getValue());
     }
     return ;
 }