| 查看: 405 | 回复: 4 | ||
蝶舞天涯ljw捐助贵宾 (小有名气)
|
[求助]
求助大神，下面定义的函数已经执行通过了，但是调用的时候，有一个错误信息 已有1人参与
|
|
BeginPackage["Einstein`"] Einstein::usage = "hhhhhhh" Begin["`private`"] Einstein[g_, v_] := Block[ {invg, dg1, dg2, dg3, G, dG, Ruv1, Ruv2, Ruv3, Ruv4, Ruv, Rrr, Result} Result = {}; (*Result save the return Values*)(*Calculate the inverse metric of \ g.*) invg = Inverse[g]; (*Calculate the affine connection.*) dg1 = Outer[D, g, v]; dg2 = Transpose[dg1, {1, 3, 2}]; dg3 = Transpose[dg1, {2, 3, 1}]; G = (1/2) invg.(dg1 + dg2 - dg3); (*Calculate the Ricci tensor.*) dG = Outer[D, G, v]; Ruv1 = Table[Sum[dG[[k, i, k, j]], {k, 4}], {i, 4}, {j, 4}]; Ruv2 = Table[Sum[dG[[k, i, j, k]], {k, 4}], {i, 4}, {j, 4}]; Ruv3 = Table[ Sum[G[[k, i, j]] G[[h, k, h]], {k, 4}, {h, 4}], {i, 4}, {j, 4}]; Ruv4 = Table[ Sum[G[[k, i, h]] G[[h, j, k]], {k, 4}, {h, 4}], {i, 4}, {j, 4}]; Ruv = Ruv1 - Ruv2 - Ruv3 + Ruv4; Result = Append[Result, Simplify[Ruv]]; (*Calculate the Curvature Scalar.*) Rrr = Sum[invg[[i, i]] Ruv[[i, i]], {i, 4}]; Result = Append[Result, Simplify[Rrr]]; (*Return the Result.*) Return[Result] ] End[] EndPackage[] 计算结果应该如下  |
回复此楼» 猜你喜欢
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
要不要辞职读博？
已经有6人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
蝶舞天涯ljw
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1367.5
- 散金: 678
- 红花: 5
- 帖子: 192
- 在线: 104.7小时
- 虫号: 2713327
- 注册: 2013-10-10
- 专业: 星系和类星体
2楼2016-04-22 08:47:28
蝶舞天涯ljw
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1367.5
- 散金: 678
- 红花: 5
- 帖子: 192
- 在线: 104.7小时
- 虫号: 2713327
- 注册: 2013-10-10
- 专业: 星系和类星体
3楼2016-04-22 08:47:43
蝶舞天涯ljw
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1367.5
- 散金: 678
- 红花: 5
- 帖子: 192
- 在线: 104.7小时
- 虫号: 2713327
- 注册: 2013-10-10
- 专业: 星系和类星体
4楼2016-04-22 08:47:54
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
蝶舞天涯ljw: 金币+10, ★★★很有帮助 2016-04-22 12:14:08
感谢参与，应助指数 +1
蝶舞天涯ljw: 金币+10, ★★★很有帮助 2016-04-22 12:14:08
|
开始那里，Result={}前面，大括号后面，少了一个逗号。Block的语法是Block[{...}, ...] 发自小木虫IOS客户端 |
5楼2016-04-22 09:04:19