数值积分（复合辛普森公式）

 

#include <stdio.h>

#include <math.h>

#define eps 1e-4//迭代精度

double f(double x)//被积函数

{

    return pow((1+pow(cos(x),2)),0.5);

}

double Simpson(double a,double b,int n,double result[2])//复合辛普森公式

{

    double h=(b-a)/n,s1=f(a+h/2),s2=0;

    for(int i=1;i<n;i+