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

#if 0

    在好几个网友的要求下，先给出一个简单的逐次分半法梯形求积公式的C语言直接实现
。不过还是提醒一下如果学习算法的话，最好还是自己实际进行编码的好。

#endif

#ifdef CODE1//g++ -DCODE1 thisfile.cpp
#include <stdio.h>
#include <math.h>
// 梯形公式的逐次分半法求积函数
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<