Question

这个积分fortran程序哪里有问题？求指导！

function trap(f, a, b, n) !梯形法求积分程序，F—被积函数，A、B—积分上下限，N—区间分划数
h=(b-a)/n !子区间宽度
s=0
do i=1, n
s=s+f(a+(i-1)*h)+f(a+i*h) !对每个子区间的积分采用梯形公式计算，并求和
end do
trap=s*h/2
end
function f(x)
f=3*x**2+x-15
end
function g(x)
g=1/(1+x**2)
end
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!上面是函数，下面是主程序，调用上面函数
program main
parameter(a=-1, b=1)
integer :: n
external f, g !定义作为实参的函数子程序名f、g的属性
write(*,*) "please input a number:"
read(*,*) n
write(*,*) trap(f, a, b, n), trap(g, a, b, n)
end
就是用fortran求积分的问题

Answer 1

这个程序没什么问题啊！我把主程序修改了一下，n从100到一亿的积分结果都列了出来，给你贴出来第一个函数结果做个参考吧！
主程序这么写：
program main
parameter(a=-1, b=1)
integer :: n, m
external f, g !定义作为实参的函数子程序名f、g的属性
! write(*,*) "please input a number:"
! read(*,*) n
! write(*,*) trap(f, a, b, n), trap(g, a, b, n)
do n = 1, 7
m = 100*10**(n-1)
write(1,*) m, trap(f, a, b, m)
enddo
end
输出的结果是：
100 -27.99960
1000 -27.99997
10000 -28.00006
100000 -28.00539
1000000 -27.98183
10000000 -26.84355
100000000 -2.684355