梯形法求数值定积分。如F(x)=∫ (4-X*X)^0.5dx。[a,b] c语言+注释

梯形法求数值积分。F(x)=∫ (4-X*X)^0.5dx。[a,b]

#include<stdio.h>
#include<math.h>
int main()
{
    
    
    double a,b,s=0,c;
    double step=1E-7;//梯形的高，dx
    
    do
    {
    
    
       scanf("%lf%lf",&a,&b);

    }while(a>=b);//定积分区间
    for(a=a+step;a<=b;a+=step)
    {
    
    
        c=a-step;//保持稳定的step距离
        s+=(sqrt(4-a*a)+sqrt(4-c*c))*0.5*step;//面积累加
    }
    printf("%lf",s);
    return 0;
}