Question

如何利用fortran程序解这个积分

Answer 1

主程序这么写：
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