计算方法实验七

#include <bits/stdc++.h>
using namespace std;
double f(double x)
{
    return (4/(1+x*x));
}
int main()
{
    double a = 0,b = 1,h=(b-a),sum,x;
    double T1 = h/2*(f(a)+f(b)),T2;        //（1）梯形公式计算积分近似值
    for(int i=0;i<9;i++)                //（2）按变步长梯形公式计算积分近似值（区间分半）
    {
        printf("h=%f,T=%f\n",h,T1);
        sum = 0;
        x = a+h/2;                       //区间分半
        while(x<b)
        {
            sum+=f(x);                    //求累加和
            x+=h;
        }
        T2 = T1/2+h/2*sum;                //（2）
        h/=2;
        T1 = T2;
    }
    printf("h=%f,T=%f\n",h,T1);
    return 0;
}