分半梯形公式c语言,梯形求积公式的逐次分半法(C语言实现)

#if0

在好几个网友的要求下，先给出一个简单的逐次分半法梯形求积公式的C语言直接实现

。不过还是提醒一下如果学习算法的话，最好还是自己实际进行编码的好。

#endif

#ifdef CODE1//g++ -DCODE1 thisfile.cpp

#include

#include

// 梯形公式的逐次分半法求积函数

double integral(

double a,//积分区间起始值

double b,//积分区间终止值

double(*f)(double),//被积函数

double eps//积分允许误差

)

{

double n = 1;

double h = (b-a)/2;

double Tn = h*(f(a)+f(b));

double T2n = 0;

while(true)

{

double sumf = 0;

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

{

sumf += f(a+(2*i-1)*h);

}

T2n = Tn/2+h*sumf;

if(fabs(T2n-Tn)<3*eps)break;

Tn = T2n;

n *= 2;

h /= 2;

}

return T2n;

}

// 自定义的被积函数，只要满足