2013-03-11 2 views
-1

다음 함수 peval2는 1로 수렴 곡선을 수렴합니다. 메이플에서. 그러나 파이썬은 잘못된 합계 값을 계산합니다. 누구든지 내가 뭘 잘못하고 있는지 힌트를 줄 수 있습니까?이 합계가 파이썬에서 수렴하지 않는 이유

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

여기 생각입니다 : 운영자를위한 우선 순위 규칙은 파이썬과 메이플에 따라 변경 될 수 있습니다을 나는대로 매개 변수를 변경했습니다. – duffymo

+0

반올림 오류 일 수 있습니까? – Cairnarvon

+1

올바른 코드를 작성했는지 어떻게 알 수 있습니까? 그것이해야 할 일이 무엇인지, 아니면 더 나은 방법으로 우리에게 정확한 메이플 코드를 보여주십시오. –

답변

2

Maple 코드와 동일한 매개 변수를 사용하고 있습니까?

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

결과 :

enter image description here

관련 문제