求助怎么解决小数溢出

hanlinbao

在做那个小练习不知道这个怎么错了。。。我没安numpy，所以想使用62的基础去做63的题。
我看网上说是math。exp的小数太多溢出了，但没教怎么解决

import math
def F(x):
        return (math.exp(abs(x))-1)/(math.exp(abs(x))+1)              
X = list([-0.9, -0.8, -0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0])
Y = list([13.1441,11.5536,10.1961,9.0416,8.0625,7.2336,6.5321,5.9376,5.4321,5.0,4.6281,4.3056,4.0241,3.7776,3.5625,3.3776,3.2241,3.1056,3.0281,3.0])
a=b=c=d=e=f=0
for repeat in range(100):
        for i in range(20):
                c1 = a*(X[i]**5) + b*(X[i]**4) + c*(X[i]**3) + d*(X[i]**2) + e*X[i]
                e1 = Y[i] - f - c1
                d1 = e1 * F(e1)
                a += d1
                b += d1
                c += d1
                d += d1
                e += d1
                c2 = a*(X[i]**5) + b*(X[i]**4) + c*(X[i]**3) + d*(X[i]**2) + e*X[i] + f
                e2 =Y[i]-c2
                d2 =e2 * F(e2)
                f += d2
print(a,b,c,d,e,f)

塔利班

我运行没问题

Sj中国智造

我也运行没问题啊

hanlinbao

Sj中国智造 发表于 2019-3-14 13:09
我也运行没问题啊

惊了，我咋不行

wp231957

e:\>python ex11.py
-0.4352252469667711 -0.4352252469667711 -0.4352252469667711 -0.4352252469667711 -0.4352252469667711 5.299108231460074


不知道结果对不对