有谁用过 R (或 Maxima，Matlab) 请帮个忙

我用Maxima （5.17.1） 算出的 Einstein tensor 给出的不是symmetric. 这不对呀。不知是不是个bug. Ricci tensor 结果是对称的。

有谁能帮忙在你的 R（或类似软件)上 run一下。

这是我的 metric tensor ,

( [1 - t^2*C^2/(M^2*r^2), C/M^2/r*t*m , 0],
[C/M^2/r*t*m, -(m/M)^2, 0],
[0, 0, -1]
);

我在 Maxima 上的代码是这个, 您参考：

load(ctensor);

(ratwtlvl:false,ratfac:true);

dim:3;

ct_coords:[t,x,r];

(ct_coords) [t,x,r];
lg:matrix([1 - t^2*C^2/(M^2*r^2), C/M^2/r*t*m , 0],
[C/M^2/r*t*m, -(m/M)^2, 0],
[0, 0, -1]
);

cmetric();

ug;

("christof() computes the Christoffel-symbols of the 1st and 2nd kind")$
christof(mcs);

ricci(true);

einstein(true);

最后的结果是这，
ein1,1 = - 5C^2*t^2/(4*M^2*r^4)
ein1,2 = C*t(2*C^2t^2 + M^2*r2)/(M^2*m*r^5)
ein1,3 = C^2*t/(2*M^2*r^3)

ein2,1 = C*m*t/(M^2*r^3)
ein2,2 = 7C^2*t^2/(4*M^2*r^4)
ein2,3 = C*m/(2*M^2*r^2)

ein3,2 = C/(2*m*r^2)
ein3,3 = C2*t^2/(4*M^2*r^4)