-1

次の関数 peval2 は、Maple などで正しい 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()
4

1 に答える 1

2

Maple コードと同じパラメータを使用していますか? パラメータを次のように変更しました。

plt.plot(rr,peval2(rr,[0.001, 0.002]))

結果:

ここに画像の説明を入力

于 2013-03-12T00:33:53.433 に答える