| 查看: 628 | 回复: 0 | |||
ldoop
|
[交流]
【求助】mathematica 求助
|
|
大家好，请帮忙看看我的这个小程序： n = Table[{140, 80, 57, 41, 31, 24, 20, 16, 13, 11, 10, 8, 7, 6, 5, 4, 4, 4, 3, 3, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1}]; f[j_, r_, x_, a_] := n[[j]] Log[r x] + LogGamma[a + 1] + LogGamma[n[[j]]] + Log[Hypergeometric2F1[n[[j]], n[[j]], n[[j]] + a, x - r x]] - Log[x] - Log[n[[j]]] - LogGamma[n[[j]] + a] - Log[HypergeometricPFQ[{1, 1, 1}, {2, 1 + a}, x] - (1 - r) HypergeometricPFQ[{1, 1, 1}, {2, 1 + a}, x - r x]]; L[r_, x_, a_] := Sum[N[f[j, r, x, a]], {j, 1, Length[n]}]; NMaximize[{L[x, r, a], 0 < x, x < 1, 0 < r, r < 0.2, a > 0}, {x, r, a}] 为什么计算不出来值呢？ 如果限制条件是r<1，可以得到值，再改成其他的就得到SystemException["MemoryAllocationFailure", 请大家帮忙看看原因？ 另外 为什么当n中含有29，30，会产生复数，也就是n=29 or 30 计算的L值是复数，应如何处理？或者避免？ 谢谢 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有15人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有13人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有20人回复
20