fortran 编写一个利用梯形法求任意函数定积分的通用外部函数 integral(fun, lower, upper)，其中，形参fun为被积函数名,lower与upper分别为定积分的下限和上限。在主程序中调用integral获得以下定积分的值。

好的，我理解了你的问题。以下是使用 Fortran 编写的利用梯形法求任意函数定积分的通用外部函数 integral(fun, lower, upper)的代码：

function fun(x)
implicit none
real :: x
fun = !这里输入被积函数!
end function fun

function integral(fun, lower, upper)
implicit none
external fun
real :: lower, upper
real :: integral, h, x
integer :: i, n
n = 1000 !将积分区间分成1000个小梯形
h = (upper - lower) / n
integral = 0.0
do i = 1, n
    x = lower + h * (i - 0.5)
    integral = integral + h * fun(x)
end do
end function integral

!在主程序中调用integral函数
program main
implicit none
real :: result
result = integral(fun, 0.0, 1.0) !计算定积分的值
write(*,*) "The result is: ", result
end program main

在代码中，首先定义了被积函数fun(x)和梯形法求积分的函数integral(fun, lower, upper)。其中，fun(x)函数需要根据不同的被积函数自行定义。在integral函数中，使用梯形法将积分区间分成1000个小梯形，使用循环计算每个小梯形的面积，并将其加入到积分值中。最后，在主程序中调用integral函数，并输入被积函数fun(x)的名字和积分区间的下限和上限，计算定积分的值。

需要注意的是，这里的被积函数fun(x)需要根据实际情况自行定义，代码中只是一个样例。另外，在使用梯形法求积分时，积分区间的划分数n应该根据被积函数的特点进行调整，以保证精度。