| 查看: 449 | 回复: 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[] 计算结果应该如下  |
回复此楼» 猜你喜欢
考研调剂
已经有3人回复
-
266求调剂
已经有5人回复
-
268求调剂
已经有6人回复
-
085600，专业课化工原理，321分求调剂
已经有6人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有3人回复
-
工科求调剂
已经有10人回复
-
材料295
已经有12人回复
-
材料专业383求调剂
已经有6人回复
-
一志愿江南大学085501机械工程专硕326分，本科佳木斯大学
已经有8人回复
-
总分328生物与医药考数学求调剂
已经有8人回复
蝶舞天涯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
- 专业: 星系和类星体
2楼2016-04-22 08:47:28
蝶舞天涯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
