#include <stdio.h>
float f(float x)
{
    return(exp(x)/(4+x*x));
}
main()
{
    float a=0,b=1;//積分上下界
    float h=b-a,T1,T2,s,x;//h是區間長度，T1是分為n等分的近似值，T2是分為2n等分的近似值
    int i;//次數
    T1=h/2*(f(a)+f(b));//公式
    for(i=0;i<4;i++)
    {
        printf("h=%f,T=%f\n",h,T1);
        s=0;
        x=a+h/2;
        while(x<b){
            s+=f(x);
            x+=h;
        }
        T2=T1/2+h/2*s;//公式
        h/=2;//將區間逐次分半
        T1=T2;//賦值後繼續迴圈
    }
    printf("h=%f,T=%f\n",h,T1);
    return;
}