以下函数peval2给出了收敛于1的饱和曲线,这是正确的,例如在枫树。然而,Python计算错误的总和值。谁能给我一个提示我做错了什么?
from numpy import *
from scipy.optimize import leastsq
import matplotlib.pyplot as plt
def peval2(x, p):
b,c = p
err = 0 * x
for ii in range(0,30):
for jj in range(0,30):
for kk in range(0,30):
err += exp(-1.0*x*pi**2.0*(b*(((2.0*ii+1.0)/75.0)**2+((2.0*jj+1)/75.0))**2+c*((2.0*kk+1)/2.0)**2.0))/((2.0*ii+1)*(2.0*jj+1)*(2.0*kk+1))**2.0
print err
return 1.0-((8.0/pi**2)**3*err)
rr = linspace(0,4000,100)
plt.plot(rr,peval2(rr,[5e-7,2e-7]))
plt.show()
答案 0 :(得分:2)
您使用与Maple代码相同的参数吗?我将参数更改为:
plt.plot(rr,peval2(rr,[0.001, 0.002]))
结果: