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

主程序这么写：

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

龙贝格积分公式Sn=(4T2n-Tn)/3,

Cn=(4^2S2n-Sn)/(4^2-1),

Rn=(4^3C2n-Cn)/(4^3-1),

其中,2n,n都是下标.

同理,依次类推.

这是在变步长求积过程中的三个加速公式,将粗糙的积分近似值迅速加工成精度较高的积分近似值的求积方法为龙为个求积算法.

这个地方输入法问题,麻烦你的眼睛辨认了！

---