r - R中的双重积分（在双重拉普拉斯反演中）

Question

我收到一条错误消息“积分可能发散”。我认为集成过程在 F_3 中中断，但我不知道具体原因。学习在 R 中执行这种双重积分会很棒。对于那些感兴趣的人，我正在尝试使用 Bromwich 轮廓执行双重拉普拉斯反演。

#  Bromwich inversion integral

        # Laplace transform of f(x,y) with parameters s and t respectively. 
        F <- function(s,t) {1/((s+3)*(t+3))}

        b <- 1.45

        F_1 <- function(x,t,u) {Re(F(complex(real=b, imaginary=u),t))*cos(u*x)}
        F_2 <- function(x,t) {2*exp(b*x)/pi*integrate(F_1, lower=0, upper=Inf, x=x, t=t, subdivisions=10000)$value}

        F_3 <- function(x,y,w) {Re(F_2(x, complex(real=b, imaginary=w)))*cos(w*y)}
        f2 <- function(x,y) {2*exp(b*y)/pi*integrate(F_3, lower=0, upper=Inf, x=x, y=y, subdivisions=10000)$value}

        # Original function
        f1 <- function(x,y) {exp(-3*(x+y))}