from math import sqrt
a=1e-8 b=10 c=1e-8 x1 = ((-b)-sqrt((b**2)-(4*a*c)))/(2*a) x2 = ((-b)+sqrt((b**2)-(4*a*c)))/(2*a) print 'x1 = {}'.format(x1) print 'x2 = {}'.format(x2) print (4*a*c) print (sqrt(b**2-4*a*c)) print b**2 print 2*a
当我运行程序时，将返回：
x1 = -1e+09 x2 = 0.0 4e-16 10.0 100.0 2e-08
我需要的是x2等于-1e-9。
问题似乎出在
sqrt((b**2)-(4*a*c))
因为它给出的结果是10，显然是因为4 （10 ^ -8）（10 ^ -8）几乎等于0，并且被python认为是0。
结果是：
sqrt((b**2)-(4*a*c)) = sqrt(b**2) = sqrt(10**2) = 10
任何帮助将不胜感激
使用十进制模块：
from decimal import Decimal a = Decimal('1E-8') b = 10 c = Decimal('1E-8') x1 = ((-b)-((b**2)-(4*a*c)).sqrt())/(2*a) x2 = ((-b)+((b**2)-(4*a*c)).sqrt())/(2*a) print 'x1 = {}'.format(x1) print 'x2 = {}'.format(x2)
结果是
x1 = -999999999.999999999000000000 x2 = -1.0000000000E-9
