当前在线人数15777
首页 - 分类讨论区 - 学术学科 - 数学版 - 同主题阅读文章

此篇文章共收到打赏
0

  • 10
  • 20
  • 50
  • 100
您目前伪币余额:0
未名交友
[更多]
[更多]
有谁用过Maxima (类似也可,如Matlab)
[版面:数学][首篇作者:oldwillow] , 2021年09月12日20:05:59 ,315次阅读,4次回复
来APP回复,赚取更多伪币 关注本站公众号:
[分页:1 ]
oldwillow
进入未名形象秀
我的博客
[回复] [回信给作者] [本篇全文] [本讨论区] [修改] [删除] [转寄] [转贴] [收藏] [举报] [ 1 ]

发信人: oldwillow (老柳树), 信区: Mathematics
标  题: 有谁用过Maxima (类似也可,如Matlab)
发信站: BBS 未名空间站 (Sun Sep 12 20:05:59 2021, 美东)

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

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

--
※ 修改:·oldwillow 於 Sep 13 09:14:43 2021 修改本文·[FROM: 198.]
※ 来源:·WWW 未名空间站 网址:mitbbs.com 移动:在应用商店搜索未名空间·[FROM: 2a02:c7f:6072:d]

 
minquan
进入未名形象秀
我的博客
[回复] [回信给作者] [本篇全文] [本讨论区] [修改] [删除] [转寄] [转贴] [收藏] [举报] [ 2 ]

发信人: minquan (三民主义), 信区: Mathematics
标  题: Re: 有谁用过Maxima (类似也可,如Matlab)
发信站: BBS 未名空间站 (Mon Sep 13 04:59:38 2021, 美东)

R

Octave
--
https://www.youtube.com/watch?v=9HznGTt495M
※ 来源:·WWW 未名空间站 网址:mitbbs.com 移动:在应用商店搜索未名空间·[FROM: 95.]

 
oldwillow
进入未名形象秀
我的博客
[回复] [回信给作者] [本篇全文] [本讨论区] [修改] [删除] [转寄] [转贴] [收藏] [举报] [ 3 ]

发信人: oldwillow (老柳树), 信区: Mathematics
标  题: Re: 有谁用过Maxima (类似也可,如Matlab)
发信站: BBS 未名空间站 (Mon Sep 13 07:30:39 2021, 美东)

如果你方便的话,不知能不能帮我 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 於 Sep 15 19:13:34 2021 修改本文·[FROM: 2a02:c7f:6072:d8]
※ 来源:·WWW 未名空间站 网址:mitbbs.com 移动:在应用商店搜索未名空间·[FROM: 198.]

 
minquan
进入未名形象秀
我的博客
[回复] [回信给作者] [本篇全文] [本讨论区] [修改] [删除] [转寄] [转贴] [收藏] [举报] [ 4 ]

发信人: minquan (三民主义), 信区: Mathematics
标  题: Re: 有谁用过Maxima (类似也可,如Matlab)
发信站: BBS 未名空间站 (Tue Sep 14 04:22:22 2021, 美东)

你到统计版问问

【 在 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);
: ...................



--
https://www.youtube.com/watch?v=9HznGTt495M
※ 来源:·WWW 未名空间站 网址:mitbbs.com 移动:在应用商店搜索未名空间·[FROM: 95.]

 
oldwillow
进入未名形象秀
我的博客
[回复] [回信给作者] [本篇全文] [本讨论区] [修改] [删除] [转寄] [转贴] [收藏] [举报] [ 5 ]

发信人: oldwillow (老柳树), 信区: Mathematics
标  题: Re: 有谁用过Maxima (类似也可,如Matlab)
发信站: BBS 未名空间站 (Tue Sep 14 05:29:14 2021, 美东)

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

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




--
※ 修改:·oldwillow 於 Sep 14 19:16:26 2021 修改本文·[FROM: 2a02:c7f:6072:d8]
※ 来源:·WWW 未名空间站 网址:mitbbs.com 移动:在应用商店搜索未名空间·[FROM: 2a02:c7f:6072:d]

[分页:1 ]
[快速返回] [ 进入数学讨论区] [返回顶部]
回复文章
标题:
内 容:

未名交友
将您的链接放在这儿

友情链接


 

Site Map - Contact Us - Terms and Conditions - Privacy Policy

版权所有,未名空间(mitbbs.com),since 1996