| 查看: 2836 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有16人回复
-
初秋的晨风
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有11人回复
-
Ei源刊怎么投
已经有4人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有6人回复
-
两块石头
已经有5人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
-
找导师
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 博后能申请到牛顿申请奖学金
已经有13人回复
- 牛顿-爱因斯坦永动机机理方面的研究能获得自然科学基金资助吗
已经有7人回复
- 【求助】用牛顿环测量球面的曲率半径
已经有3人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
jjdg:把方程具体的值给出来 2010-04-20 12:52
|
方程就是f = z ^ 3 + xl * z ^ 2 + xm * z + xn=0 其中xl,xm,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 最后得到的结果只有一个根。 |
2楼2010-04-20 10:05:18
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
6楼2010-04-20 23:20:50
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
9楼2010-04-21 11:06:08
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
12楼2010-04-22 17:50:36
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
13楼2010-04-22 17:55:18
liangzidou
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 371
- 红花: 1
- 帖子: 68
- 在线: 15.8小时
- 虫号: 720428
- 注册: 2009-03-11
- 性别: GG
- 专业: 能源化工
15楼2010-04-24 23:23:15