用梯形法求定积分程序实现例题——c——新人文

 

 

以下为梯形法原理。

 

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

#include<stdio.h>
#include<math.h>

int main()
{
	double step=1E-7,a,b,c,sum=0;
	scanf("%lf%lf",&a,&b);
	
	for( a=a+step;a<=b;a+=step ){
		c=a-step;
		sum+=(sqrt( 4-a*a )+sqrt( 4-c*c ))*0.5*step;
	}
	printf("%lf\n",sum);
	
	return 0;
}

其中step为步长，可以理解为精度控制