2楼: Originally posted by mchen10 at 2011-10-06 07:58:39:
不是这个领域的 纯属门外汉 不保证我说的能用呵..

简化起见，把你的方程写成以下形式:

a^(1/3) + a/2600 -5 = 0

每一项除以a^(1/3)变为：
1+a^(2/3) - 5a^(-1/3) = 0

令x = a^(1/3)
有1 + x^2/2600 ...

Python Code:

x=0
while True:
  y=13000.0/(x*x+2600.0)
  if abs(y-x) < 0.001:
    print "y=", y
    print "x=", x
    break
  x=y