復化simpson求積算法 c++實現

考慮到浮點數，採用double類型。

double f(double x)
{
    //公式
}
double count(double n,double a,double b)
{
    double sum_in=0;
    double h=(b-a)/n;
    for(double i=0;i<n;i++)
    {
        if((int)i%2==0)
        {
            sum_in+=2*(f(i*h+a));
        }
        else
        {
            sum_in+=4*(f(i*h+a));
        }
    }
    double sum=0;
    sum=h/3*(f(a)+f(b)+sum_in);
    return sum;
}