我寫了一個計算Riemann-Siegel Z(t)函數的基本程序。我很好奇,如果有更好的方法來近似剩餘期限。我現在使用的方法使用了來自Haselgrove的可怕表格近似值。 Further information about the Riemann Siegel formula。這可能稍微先進一點,但有關這方面的更多詳細信息,請參見this thesis。另外,在Edwards book。 我知
我試圖實現euler's method來近似py的值e。這是我到目前爲止有: def Euler(f, t0, y0, h, N):
t = t0 + arange(N+1)*h
y = zeros(N+1)
y[0] = y0
for n in range(N):
y[n+1] = y[n] + h*f(t[n], y[n])
f