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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » mathematic中使用while循环求方程，出现超过 1024 的递归深度的问题。
北京石油化工学院2026年研究生招生接收调剂公告
41/1返回列表
查看: 2590  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +11 倘若起风了呢 2026-04-05 11/550 2026-04-05 23:21 by 来看流星雨10
[考研] 材料工程302分求调剂 +9 zyx上岸！ 2026-04-04 9/450 2026-04-05 22:08 by 醉翁wl
[考研] 求调剂 +6 chenxrlkx 2026-04-05 6/300 2026-04-05 22:05 by imissbao
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 化学0703-一志愿211-338分求调剂 +7 vants 2026-04-05 7/350 2026-04-05 18:17 by cql1109
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-03-31 5/250 2026-04-04 19:45 by 1753564080
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +9 万事宜臻 2026-04-04 9/450 2026-04-04 11:49 by 啵啵啵0119
[考研] 282求调剂 +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 285求调剂 +14 AZMK 2026-04-02 14/700 2026-04-02 15:54 by 上九天揽月（好�
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业 +7 zju51 2026-03-31 9/450 2026-04-01 19:35 by CFQZAFU
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
信息提示
关闭
请填处理意见
关闭