我正在寻找python的Jacobi theta函数的导数的实现。我发现了这个
http://mpmath.googlecode.com/svn/trunk/doc/build/functions/elliptic.html#jtheta
然而,导数是针对错误的论点计算的。我想要关于q的导数,而不是z。有人知道得到这个的最好方法吗?仅采用有限差分的数字导数并不能很好地工作(我使用的ODE求解器,当我这样做时会崩溃)。
(实际上,如果可能的话,我也只想要特殊情况下的z=0。对于E4 Eisenstein系列来说,同样的问题也足够了)
谢谢!
发布于 2015-03-31 12:25:04
这种关系应该对您有帮助:http://functions.wolfram.com/EllipticFunctions/EllipticTheta2/13/01/0002/ equation
def jtheta_dq(n, z, q):
# cf. http://functions.wolfram.com/EllipticFunctions/EllipticTheta2/13/01/0002/
return -mpmath.jtheta(n, z, q, 2)/(4*q)https://stackoverflow.com/questions/29364762
复制相似问题