欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

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

程序员文章站 2022-03-10 12:09:36
...

梯形法求数值积分。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;
}