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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » mathematic中使用while循环求方程，出现超过 1024 的递归深度的问题。
北京石油化工学院2026年研究生招生接收调剂公告
41/1返回列表
查看: 2591  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +12 哦呦呼o 2026-04-04 12/600 2026-04-06 23:04 by chenzhimin
[考研] 一志愿武理车辆专硕总分 281 求调剂 +3 上岸研究生. 2026-04-02 3/150 2026-04-06 16:46 by 小月儿呀呀
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 求助071001调剂！！！ +4 黄守松 2026-04-05 5/250 2026-04-06 10:55 by 1028907439
[考研] 生物与医药086000调剂一志愿西北农林320分 +3 美美女士 2026-04-03 3/150 2026-04-05 21:55 by 学员8dgXkO
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 292分，材料与化工，申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 309分085801求调剂 +11 MY_angel 2026-03-31 11/550 2026-04-04 19:11 by 蓝云思雨
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 085601一志愿北理325分求调剂 +6 找调剂，， 2026-04-02 6/300 2026-04-03 22:20 by 柟�?
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 296求调剂 +4 汪！？！ 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
信息提示
关闭
请填处理意见
关闭