0

私は次のコードがあると思っていました：

[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 のすべての答えが無限になるだけですが、理由がわかりません。

4

1 に答える 1

2

あなたはドットを忘れました。 r^-1と同じinv(r)です。のすべてrの行は最初の行と同じであるため、削減された行のエシュロン形式はほぼすべてゼロになり、逆は無限の値になります。私はあなたが意味したと信じています

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

または、より明確に私の意見では、

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

の残りについても同様ですr^-1。

于 2016-08-08T22:44:03.917 に答える