mathematica大家帮帮忙看这个程序哪里出错了？

dongyufeng

" O& C. g6 [' e/ \, s
我用二分法解方程组h1[x1,x2]和h2[x1,x2],可是出现的“General::stop: Further output of PolynomialGCD::lrgexp will be suppressed during this calculation. >>”没办法解决，我不知道怎么调试。热心人谁帮帮忙啊，到底哪里出错了。而且解出的x1的最优值竟然会随着它的范围变化而变化（b1），按理来说方程组的解是唯一的呀（我要的是大于0的解），为什么呀？
/ A( O' T# Y0 ^/ t2 w6 C附件里是我的程序，求帮忙啊。 谢谢！！

linanisyugioh

你这个设置太坑了，下载附件居然要扣体力。
是不是说帮你分析的代价是扣我体力，那我多划不来。

dongyufeng

linanisyugioh 发表于 2013-4-19 17:13
2 N1 F. H% B0 }0 K& J你这个设置太坑了，下载附件居然要扣体力。
8 _0 Q+ R; B) m是不是说帮你分析的代价是扣我体力，那我多划不来。

# E' `2 ~& _) J) W2 H- V" p这个我还真没注意，呵呵，不好意思啊

dongyufeng

linanisyugioh 发表于 2013-4-19 17:13
你这个设置太坑了，下载附件居然要扣体力。
是不是说帮你分析的代价是扣我体力，那我多划不来。

第一次用这个网站，不知道附件下载要扣体力，也不懂怎么设置......你能帮我看看不，就是浪费了你2体力了，，，

linanisyugioh

貌似mathematica自带findroot函数可以用来求数值解哦。
  y  F% @9 a" N" y" z$ q6 c4 E好吧，也许这个不能满足你的需求呢。

linanisyugioh

另外，你能不能把你的问题描述详细点，比如说你要解的方程组具体是多少，有哪些参数，具体问题，疑惑在哪里，光给个nb文件，而且你的程序又极有可能是错的，如果程序很长很费解的话，那么很难让人真正明白你想要搞什么，只能去猜。这样很浪费精力的，而且这些精力都是不必要的，只要你能把问题尽可能的描述一翻，就算人家不能帮你找出问题，也有可能给出你一种可行的方案。总之呢，网上流行一篇《提问的艺术》的帖子，可以去看看的。

Valentiner

干嘛要用二分法解呢，用切线法或割线法不是更好更准确吗   用FindRoot命令

Valentiner

干嘛要用二分法解呢，用切线法或割线法不是更好更准确吗   用FindRoot命令

dongyufeng

linanisyugioh 发表于 2013-4-24 18:01
另外，你能不能把你的问题描述详细点，比如说你要解的方程组具体是多少，有哪些参数，具体问题，疑惑在哪里 ...

/ U; K8 }7 i: n) b2 P  n8 j你的意见真的是很实用，谢谢了！《提问的艺术》我看过了，果真是我太二了，呵呵。

linanisyugioh

好像更普遍的译法是《提问的智慧》。我搞错了。抱歉了，大家要搜的话，搜《提问的智慧》吧，可以看到更完全的版本。