利用integral2求二重积分遇到问题。

代码如下：

syms y z;

alpha=0.1;

ita=1;

PS1=1;

PS2=1;

theta1=0.01;

theta2=0.01;

d1=1;

d2=1;

a=ita*PS1^2*alpha;

b=ita*d1*PS1*theta1*alpha;

c=d1*d2*PS1*theta2*(1-alpha);

fun=@(y,z)exp(-z-c*y/(a*z-b*y))*(a*c*z^2)/(a*z-b*y)^2;

zmin=@(y)b*y/a;

C12=integral2(fun,0,inf,zmin,inf)

报的错误如下：错误使用 ^

输入必须为标量和方阵。

要按元素进行 POWER 计算，请改用 POWER (.^)。

明明没有用到矩阵和数组啊 我按照他的错误改正之后 又出现错误为：*两边的矩阵维度不一致

请帮忙大神指导一下

Answer 1

[腾龙公司开户【TL9704.com】]

fun=@(y,z)exp(-z-c*y./(a*z-b*y)).*(a*c*z.^2)./(a*z-b*y).^2;

C12=integral2(fun,0,inf,zmin,inf)