Python のヤコビ シータ関数の導関数の実装を探しています。私はこれを見つけました
http://mpmath.googlecode.com/svn/trunk/doc/build/functions/elliptic.html#jtheta
ただし、導関数は間違った引数に対して計算されます。z ではなく、q に関する微分が必要です。これを取得する最良の方法を知っている人はいますか？微分を有限差分で数値的に取るだけではうまくいきません (これを行うと、私が使用している ODE ソルバーがクラッシュします)。
（実際には、可能であれば、z = 0の特別なケースでもそれが必要です。E4 Eisensteinシリーズの同じ質問でも十分です）
ありがとう！