[求助]：如何编程解一元三次的状态方程

我爱大自然

求教各位，如何编程解一元三次的状态方程

目不识丁

状态方程是三次方程式，不是解状态方程吧？

WYQ4517

''牛顿迭代法解方程
''注意,迭代的初值离真值太远，可能不收敛而求不出解,
''导地线计算时,初值可设为X=110     (N/mm2)

Private Sub Command1_Click() '
    Dim y1, x1, e, Y, A, B As Single
    X = 11                                           '令X = 11作为迭代的初值
A = 45: B = 27                              '系数A,B已知,此行值举例而已
    e = 1
    While Abs(e) > 0.001               '设求出的解的精度
        Y = X ^ 3 + A * X * X - B        '状态方程X ^ 3 + A * X * X - B=0
        y1 = 3 * X * X - 2 * A * X        '一阶导数
   x1 = X - Y / y1                        '横坐标为X点的切线交X轴于X1
        e = (X - x1) / x1                      'X、X1的接近程度
   X = x1                                      '让X = x1,循环,即再作切线,找横坐标为X1点的切线交X轴的点
  Wend
    Print X
End Sub

我爱大自然

wyq4517，谢谢你啦，昨天我问的关于如何编程求临界档距，你可否给我传一个程序，我总也做不对，很着急，麻烦你了，我的qq 273667431，邮箱 weijinlibb@163.com

WYQ4517

你将解一元三次方程求应力的子程序弄对了再说吧。

invisual

甭费那老劲了，随便用个数学软件matlab什么的就轻松搞定了，我自己用的就是mathematica计算应力弧垂曲线的

WYQ4517

最初由 我爱大自然 发布
[B]wyq4517，谢谢你啦，昨天我问的关于如何编程求临界档距，你可否给我传一个程序，我总也做不对，很着急，麻烦你了，我的qq 273667431，邮箱 weijinlibb@163.com [/B]

我已在你的另一贴中（见链接）给你传去了。不过忘传排序子程序，在此补上。
http://www.xdcad.net/forum/showthread.php?s=&postid=2278727#post2278727

Private Sub paixu()  ''排序
   If aa > p Then
       p = aa
       Y = i
   End If
End Sub

chiu

可用EXCEL"工具"的"单变量求解"

xuesong709

这些很简单的!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

WYQ4517

最初由 我爱大自然 发布
[B]wyq4517，谢谢你啦，昨天我问的关于如何编程求临界档距，你可否给我传一个程序，我总也做不对，很着急，麻烦你了，我的qq 273667431，邮箱 weijinlibb@163.com [/B]

看来1楼并不着急！
他的本贴子，还有个求助临界档距的计算及判别的贴子，我费了很大的劲传了原代码，他却走人了。