我写了一个计算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