在Maxima中遇到计算问题

问题描述

在计算某些度量的Kretschmann标量时遇到了一些困难。这是我的代码：

kill(all);
load("ctensor");

dim: 4;
ct_coords: [t,r,theta,phi];

depends(m,r);
depends(A,[r,theta]);
depends(B,theta]);
depends(C,theta]);
depends(G,theta]);
depends(H,theta]);

lg: matrix([-(1-2*m(r)/r+a^2*A(r,theta)*m(r)),-a*m(r)/r*G(r,theta)],[0,1/(1-2*m(r)/r)+a^2*B(r,theta),0],r^2+a^2*C1(theta),[-a*m(r)/r*G(r,r^2*sin(theta)^2]);

cmetric(false);
christof(false);
riemann(false);
lriemann(false);
uriemann(false);
rinvariant();

这正常工作，但是如果我通过在最后一项上添加一个常量来稍微改变指标lg，则如下所示：

lg: matrix([-(1-2*m(r)/r+a^2*A(r,r^2*sin(theta)^2+a^2]);

（注意a^2）不再起作用，并产生以下输出

PQUOTIENT：较高次多项式的商（案例2b）- 错误。要调试，请尝试：debugmode（true）

错误是什么意思？错误在哪里？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）