由度规算曲率标量

上周被某人说成是“不务正业”,这次就讲讲本行吧(忽悠新手,外行可以无视,同行尽管鄙视)。

当代宇宙学的理论基石是在Robertson-Walker度规下由Einstein场方程推导出的Fridemman方程,现在的绝大部分工作都是以此为起点,而且也得到了可靠的观测支持。从给定度规得出场方程的具体形式是广义相对论的基础内容,但是由于方法繁复,教科书中都不会给出具体的计算过程,而结果又不是一望便知的,学到这里谁都少不了课下的一番推导验证。我当年偷懒跳过,现在却发现自己怎么都算不对了……

怅惘之际在论文库里发现2000年《上海天文台年刊》第21期中有一篇《利用Mathematica软件表示真空Einstein场方程》的文章,但又不想为此学门新语言,便用Matlab仿写了一个,代码如下:

虽然编写的时候问题多多,但总算是得出正确结果了,看来Matlab的符号计算功能也不差么,为什么就没有人用呢?兴奋之余又试了一下传说中可以改善显示效果的pretty函数,看到的便是左边的结果……

我半分钟之后才反应过来这居然就是3R”/R(显然是BBS泡的太少,对ASCII码不够敏感),那一刻我终于明白为什么没人用Matlab做符号计算了……

当下决定学Mathematica,周建峰等人的代码也贴上来备查吧,shift+回车直接执行。

附:

订阅评论
提醒

4 评论
最旧
最新 最多投票
内联反馈
查看所有评论
请问下,我按周建锋的给的程序包算得不对啊
2009 年 11 月 27 日 10:12

我按周建锋的给的程序包算得不对啊,请指教下,不知什么问提哦

李豪豪
回复给  请问下,我按周建锋的给的程序包算得不对啊
2017 年 8 月 27 日 19:18

这个不是程序包,就是mma的一个程序。这个程序中没有乘号。mma可以用空格代替乘号,这个程序中该有空格的地方没有空格。

2009 年 11 月 29 日 21:52

楼主太帅啦!