我是python的新手，正在嘗試以數字方式集成函數。一切似乎都正常，但是我得到的結果與我在Mathematica中得到的結果（我知道是正確的）有很大不同。有人可以幫我弄清楚發生了什么嗎？這是代碼：def integrand(x, d, a, b, l, s, wavelength, y):    return b*(np.sinc((np.pi*a/(wavelength*s))*(y + s*b*x/l))**2)*np.cos((np.pi*d/(wavelength*s))*(y + s*b*x/l))**2def intensity(y):    result, error = si.quad(integrand, -1/2, 1/2, epsrel = 1e-16, epsabs = 1e-16,                            args=(0.0006, 0.000150, 0.000164, 0.8, 1.06, 0.0000006328, y))    return result例如，如果我計算強度（0），我會在 python 中得到 0.0001580120220796804，在 Mathematica 中得到 0.000158898。在 0.5% 以內，所以這看起來沒問題。但是如果我計算強度（0.001），我會在 python 中得到 1.8729902318383768e-05，在 Mathematica 中得到 0.00012034，它們相差近一個數量級。請注意，我已嘗試減少絕對誤差和相對誤差，但這沒有任何效果。任何幫助，將不勝感激。這是Mathematica代碼：NumInt[d_, a_, b_, l_, s_, lambda_, y_] := NIntegrate[b Sinc[(a Pi/(s lambda)) (y - (s*b*      x/l))]^2 Cos[(d Pi/(s lambda)) (y - (s*b*x/l))]^2, {x, -1/2,1/2}]然后NumInt[0.0006, 0.000150, 0.000164, 0.8, 1.06, 0.0000006328, 0.001]