有谁用过Maxima （类似也可，如Matlab)

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

有谁帮忙在你的类似软件上run一下。我可以贴我的代码和结果。你 run 一下就行。就几行。

R

Octave

如果你方便的话，不知能不能帮我 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)

【 在 minquan (三民主义) 的大作中提到: 】
: R
: Octave

你到统计版问问

【 在 oldwillow (老柳树) 的大作中提到: 】
: 如果你方便的话，不知能不能帮我 run 一下 (不才对这两个软件知之甚少) ?
: 这是metric tensor,
: ( [1 - t2*C2/(M2*r2), C/M2/r*t*m , 0],
: [C/M2/r*t*m, -(m/M)2, 0],
: [0, 0, -1]
: );
: 我的 Maxima 代码是这个, 您参考：
: load(ctensor);
:
: (ratwtlvl:false,ratfac:true);
: ...................

那的人主要用 sas 和 R 吧。据愚所知这两个都没有算tensor 的 package （或叫
module）

【 在 minquan (三民主义) 的大作中提到: 】
: 你到统计版问问