应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » mathematic中使用while循环求方程，出现超过 1024 的递归深度的问题。
41/1返回列表
查看: 2597  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

c736737270

新虫 (初入文坛)

[求助] mathematic中使用while循环求方程，出现超过 1024 的递归深度的问题。

求大神解决啊，代码如下，复制黏贴即可运

g[t_] := \!\(
\*SubsuperscriptBox[\(\[Integral]\), \(td\), \(t\)]\(\((0.2 +
     0.1 \((v - td)\))\) \[DifferentialD]v\)\)

m[nu_] := 1 - E^(-0.01*nu)

beta[x_] := E^-x

td := 1/24;

nu := 69.3147

newnu := nu; oldnu := newnu - 5

k := 120


c := 20

h := 3

DE := 1000

p := 35

PPi := 5

s := 4

\[CapitalLambda]1[t1_, T_, nu_] := ((p - c)*DE*(t1 + \!\(
\*SubsuperscriptBox[\(\[Integral]\), \(t1\), \(T\)]\(\((beta[
           T - t])\) \[DifferentialD]t\)\)) - k - (c + h*td)*DE*\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(td\), \(t1\)]\(\((\((\(-1\) +
\*SuperscriptBox[\(E\), \(\(-0.008246527777777776`\)\
\*SuperscriptBox[\(E\), \(\(-0.01`\)\ nu\)]\)])\) +
         0.19583333333333333`\
\*SuperscriptBox[\(E\), \(\(-0.008246527777777776`\)\
\*SuperscriptBox[\(E\), \(\(-0.01`\)\ nu\)] - 0.01`\ nu\)]\ t +
         0.019175347222222222`\
\*SuperscriptBox[\(E\), \(\(-0.008246527777777776`\)\
\*SuperscriptBox[\(E\), \(\(-0.01`\)\ nu\)] - 0.03`\ nu\)]\ \((1.`\
\*SuperscriptBox[\(E\), \(0.01`\ nu\)] + 2.607514712539611`\
\*SuperscriptBox[\(E\), \(0.019999999999999997`\ nu\)])\)\
\*SuperscriptBox[\(t\), \(2\)] + 0.009791666666666667`\
\*SuperscriptBox[\(E\), \(\(-0.008246527777777776`\)\
\*SuperscriptBox[\(E\), \(\(-0.01`\)\ nu\)] -
            0.05`\ nu\)]\ \((0.12783564814814813`\
\*SuperscriptBox[\(E\), \(0.020000000000000004`\ nu\)] + 1.`\
\*SuperscriptBox[\(E\), \(0.030000000000000002`\ nu\)])\)\
\*SuperscriptBox[\(t\), \(3\)])\) \[DifferentialD]t\)\) +
    h*DE*td^2/2 - h*DE*td*t1 - h*DE*\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(td\), \(t1\)]\(\((
\*SuperscriptBox[\(E\), \(\(-0.08`\)\ nu\)]\ \((1.`\
\*SuperscriptBox[\(E\), \(0.08`\ nu\)]\ t1 + 0.09791666666666667`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] + 0.006391782407407406`\
\*SuperscriptBox[\(E\), \(0.06000000000000001`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 0.016666666666666666`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 0.0003129310136959876`\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)] + 0.002447916666666667`\
\*SuperscriptBox[\(E\), \(0.06`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)])\) +
\*SuperscriptBox[\(E\), \(\(-0.08`\)\ nu\)]\ \((\(-1.`\)\
\*SuperscriptBox[\(E\), \(0.08`\ nu\)] - 0.19583333333333333`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)]\ t1 - 0.019175347222222222`\
\*SuperscriptBox[\(E\), \(0.060000000000000005`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] - 0.0012517240547839502`\
\*SuperscriptBox[\(E\), \(0.05000000000000001`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] - 0.0032638888888888887`\
\*SuperscriptBox[\(E\), \(0.060000000000000005`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] - 0.00006128232351546424`\
\*SuperscriptBox[\(E\), \(0.04`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)] - 0.0004793836805555556`\
\*SuperscriptBox[\(E\), \(0.049999999999999996`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)])\)\ t +
\*SuperscriptBox[\(E\), \(\(-0.08`\)\ nu\)]\ \((0.09791666666666667`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)] + 0.019175347222222222`\
\*SuperscriptBox[\(E\), \(0.06`\ nu\)]\ t1 - 0.05`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)]\ t1 + 0.0018775860821759259`\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] - 0.004895833333333334`\
\*SuperscriptBox[\(E\), \(0.060000000000000005`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] + 0.00012256464703092846`\
\*SuperscriptBox[\(E\), \(0.04000000000000001`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 5.421010862427522`*^-20\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] - 0.0008333333333333334`\
\*SuperscriptBox[\(E\), \(0.060000000000000005`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 6.000560844222541`*^-6\
\*SuperscriptBox[\(E\), \(0.030000000000000002`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)] + 0.00003129310136959877`\
\*SuperscriptBox[\(E\), \(0.039999999999999994`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)] - 0.00012239583333333335`\
\*SuperscriptBox[\(E\), \(0.049999999999999996`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)])\)\
\*SuperscriptBox[\(t\), \(2\)] +
\*SuperscriptBox[\(E\), \(\(-0.08`\)\ nu\)]\ \((0.`  -
            0.006391782407407405`\
\*SuperscriptBox[\(E\), \(0.06`\ nu\)] + 0.03333333333333334`\
\*SuperscriptBox[\(E\), \(0.07`\ nu\)] - 0.0012517240547839504`\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\ t1 + 0.009791666666666667`\
\*SuperscriptBox[\(E\), \(0.06`\ nu\)]\ t1 - 0.00012256464703092848`\
\*SuperscriptBox[\(E\), \(0.04000000000000001`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] + 0.0009587673611111112`\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\
\*SuperscriptBox[\(t1\), \(2\)] - 8.000747792296719`*^-6\
\*SuperscriptBox[\(E\), \(0.030000000000000013`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 0.00004172413515946501`\
\*SuperscriptBox[\(E\), \(0.04000000000000001`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] + 0.00016319444444444446`\
\*SuperscriptBox[\(E\), \(0.05`\ nu\)]\
\*SuperscriptBox[\(t1\), \(3\)] - 3.9170327733119356`*^-7\
\*SuperscriptBox[\(E\), \(0.020000000000000004`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)] + 0.000023969184027777782`\
\*SuperscriptBox[\(E\), \(0.039999999999999994`\ nu\)]\
\*SuperscriptBox[\(t1\), \(4\)])\)\
\*SuperscriptBox[\(t\), \(3\)])\) \[DifferentialD]t\)\) - s*DE*\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(t1\), \(T\)]\(\((\((T -
           t)\) beta[T - t])\) \[DifferentialD]t\)\) - PPi*DE*\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(t1\), \(T\)]\(\((1 -
         beta[T - t])\) \[DifferentialD]t\)\) - T*nu)/T

While[Abs[newnu - oldnu] > 10^-4,
xx := FindRoot[{D[\[CapitalLambda]1[t1, T, nu], t1] == 0,
    D[\[CapitalLambda]1[t1, T, nu], T] == 0}, {t1, 0.2}, {T, 0.2}];
t1 := xx[[1, 2]]; T := xx[[2, 2]];
yy := FindRoot[D[\[CapitalLambda]1[t1, T, nu], nu] == 0, {nu，80}];
nu := yy[[1, 2]]; oldnu := newnu; newnu := nu]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-02-27 16:52:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
如果真是 递归或 迭代的深度 改变 $RecursionLimit 或 $IterationLimit 值就可以吧
不过 感觉 是代码的问题 不是深度不够的问题
代码太乱 copy过去还是乱
赞一下 回复此楼
2楼2014-02-27 20:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c736737270

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by walk1997 at 2014-02-27 20:28:03
如果真是 递归或 迭代的深度 改变 $RecursionLimit 或 $IterationLimit 值就可以吧
不过 感觉 是代码的问题 不是深度不够的问题
代码太乱 copy过去还是乱

恩 已经解决了 谢谢回复  金币就给你吧  哈哈
赞一下 回复此楼
3楼2014-02-28 18:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c736737270

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by walk1997 at 2014-02-27 20:28:03
如果真是 递归或 迭代的深度 改变 $RecursionLimit 或 $IterationLimit 值就可以吧
不过 感觉 是代码的问题 不是深度不够的问题
代码太乱 copy过去还是乱

怎么不能给金币啊？
赞一下 回复此楼
4楼2014-02-28 18:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c736737270 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854调剂 +6 长弓傲 2026-04-09 6/300 2026-04-09 11:16 by 2451080013
[考研] 化学调剂求助 +14 LULONG1 2026-04-03 19/950 2026-04-09 10:43 by chenxi233
[考研] 一志愿武理车辆 281 求调剂 +3 上岸研究生. 2026-04-07 3/150 2026-04-09 10:26 by liuhuiying09
[考研] 一志愿北京2，材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 302分求调剂 +4 凡语祈愿 2026-04-08 5/250 2026-04-08 22:03 by 土木硕士招生
[考研] 290求调剂085701 +20 1314捧花 2026-04-02 21/1050 2026-04-08 21:19 by cyh—315
[考研] 269求调剂 +6 啊啊我我 2026-04-07 6/300 2026-04-08 20:04 by 我减肥1
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 求调剂 +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[考研] 312求调剂 +18 gtw1 2026-04-06 20/1000 2026-04-07 18:16 by 蓝云思雨
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 生物学学硕求调剂：351分一志愿南京师范大学生物学专业 +6 …～、王…～ 2026-04-06 7/350 2026-04-06 18:54 by macy2011
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 生物学308分求调剂（一志愿华东师大） +8 相信必会光芒万� 2026-04-05 10/500 2026-04-05 12:19 by Hdyxbekcb
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭