| 查看: 2292 | 回复: 15 | |||
| 本帖产生 2 个 程序强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
liangzidou银虫 (小有名气)
|
[交流]
【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】 已有5人参与
|
||
|
我在用vb编写中，遇到这样一个问题，解一个一元三次方程只得到一个根，实际上不只一个根。请各位高手帮忙找一下错误。 Private Sub Command1_Click() Dim r, p, t, w, pc, tc, m, a, b, ca, cb, xl, xm, xn As Double r = 8.314 p = 101.325 t = 301 pc = 3384 tc = 460.4 w = 0.227 m = 0.37464 + 1.5422 * w + 0.2699 * w ^ 2 tr = t / tc a = 0.45723553 * r ^ 2 * tc ^ 2 * (1 + m * (1 - Sqr(tr))) ^ 2 / pc b = 0.077796074 * r * tc / pc Print "a="; a, "b="; b ca = a * p / (r ^ 2 * t ^ 2) cb = b * p / (r * t) xl = cb - 1 xm = -3 * cb ^ 2 - 2 * cb + ca xn = cb ^ 3 + cb ^ 2 - ca * cb Print "A="; ca, "B="; cb Print "xl="; xl, "xm="; xm Print "xn="; xn Dim f, z As Single z = 0.01 Do z = z1 f = z ^ 3 + xl * z ^ 2 + xm * z + xn f1 = 3 * z ^ 2 + 2 * xl * z + xm z1 = z - f / f1 Loop While Abs(z - z1) >= 0.0005 Print z1 End Sub [ Last edited by liangzidou on 2010-5-15 at 22:22 ] |
回复此楼» 猜你喜欢
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有7人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有7人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 博后能申请到牛顿申请奖学金
已经有13人回复
- 牛顿-爱因斯坦永动机机理方面的研究能获得自然科学基金资助吗
已经有7人回复
- 【求助】用牛顿环测量球面的曲率半径
已经有3人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
10楼2010-04-21 11:18:07