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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » mathematic中使用while循环求方程，出现超过 1024 的递归深度的问题。
41/1返回列表
查看: 2601  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕求调剂 希望是a区院校 +19 好好休息好不好 2026-04-09 22/1100 2026-04-10 11:01 by luosha500
[考研] 一志愿北理工298英一数二求调剂 +8 Reframe 2026-04-10 9/450 2026-04-10 09:59 by lbsjt
[考研] 080100力学316求调剂 +7 L_Hairui 2026-04-07 7/350 2026-04-10 09:45 by vgtyfty
[考研] 301求调剂 +6 静静想想 2026-04-05 6/300 2026-04-10 09:15 by Delta2012
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 材料295 +15 小英11 2026-04-03 16/800 2026-04-09 15:44 by only周
[考研] 生物学308分求调剂（一志愿华东师大） +13 相信必会光芒万� 2026-04-06 16/800 2026-04-09 13:54 by 徐良白眉大侠
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭