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