使用Quadpack进行多维集成

Question

我想使用Quadpack在Fortran中集成三维函数：

function f(x,y,z)
...

我在考虑这样的事情：

result = Integrate(B(x), a_x, b_x)
B(x0) = Integrate(C(x0,y), a_y, b_y)
C(x0, y0) = Integrate(f(x0,y0, z), a_z, b_z)

Quadpack期望只有一个参数的函数。我如何在Fortran中实现它？模块可能会起作用，但有更容易的东西吗？