| 查看: 532 | 回复: 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[] 计算结果应该如下  |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有8人回复
-
filecode与中标关系的预测
已经有5人回复
-
这样的filecode谁见过
已经有19人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
蝶舞天涯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
100