使用梯形法求解f(x)在[a,b]上的定积分,f(x)=1/(1+x)。为了求得图形面积,先将区间[a,b]分成n等分,每个区间的宽度为h=(b-a)/n,对应将图形分成n等分,再分别求小曲边梯形的面积并累加得到积分面积。
代码如下:
real a,b,h,x,s,f0,f1,si
integer n,i
read *,a,b,n
h=(b-a)/n
x=a;s=0;f0=1.0/(1+x)
do i=1,n
x=x+h
f1=1.0/(1+x)
si=(f0+f1)*h/2
s=s+si
f0=f1
end do
print *,s
end