Matlab 范围内的双重积分

Question

我正在尝试使用以下代码在 Matlab 中进行双重集成。

%Parameter Values
theta1=0.5;

    t1=0;
    t2=1;
    delta=0.1;
    
    f2=@(x,y)((exp(-(x./theta1)+(y./theta1))*exp(-delta*(x/theta1).*(y/theta1)).*((1+delta.*(x/theta1))*(1+delta.*(y/theta1))-delta))/theta1^2);
    
    Pd=integral2(f2,t1,inf,t2,inf);

但我不断收到错误。 我的 Matlab 错误显示“integral2 错误（第 106 行）Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct)；”

请告诉我是否犯了任何错误。 我也用 python 试过，但它说它超出了范围。

这也是我的 Python 代码：

import math
positive_infnity = float('inf') 

from scipy import integrate
theta1=0.5;
delta=0.1;

f=lambda y,x:((math.exp(-(x/0.5)+(y/0.5))*math.exp(-delta*(x/0.5)*(y/0.5))*((1+(delta*x)/0.5)*(1+(delta*y)/0.5)-delta))/0.25);

print(integrate.dblquad(f, 0,positive_infnity, lambda x: 0, lambda x: positive_infnity))

运行后我得到一个数学范围错误