Question

我只是想知道我有以下代码：

[r,R] = meshgrid(2:0.01:10,-4:0.01:4);
Z=(((R).^2)-1 );
   A=2+ (R.*(r+r^-1));
   B=sqrt(Z).*(r-r^-1);
   C=((r+r^-1)./(sqrt(Z))).*acoth(A./B);
    figure
contourf(r,R,C);
    h=colorbar;
    set(get(h,'ylabel'),'string','\gamma P_P L')
    xlabel('$r$','Interpreter','latex','FontSize',18)
    ylabel('$D \over 2\sqrt{M}$','Interpreter','latex','FontSize',18)

我只是得到A和B所有答案无限但我不确定为什么？

Answer 1

你忘了点。 r^-1与inv(r)相同。并且由于r的所有行都与第一行相同，因此它的缩减行梯形将几乎全部为零，并且逆将具有无限值。我相信你的意思

A=2+ (R.*(r+r.^-1));

或者，在我看来更清楚，

A=2+ (R.*(r+ 1./r));

等等r^-1。

的剩余部分

轮廓图显示无结果

1 个答案: