如何利用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

（计算0～π亦即0°～180°范围内的积分，精确结果为2）

        s=0

        pai=3.1415926

        step=pai/180

        a=sin(0.0)

   悉稿型     do 10 i=0,179

        b=sin(step*(i+1))

        s=s+0.5*(a+b)*step

        a=b

10    敬伏  continue

        write(*,*)s

        stop

        睁猜end