如何在python中解决雅可比椭圆函数

问题描述

我知道如何在 Mathematica 中使用 EllipticF 求解：

EllipticF[ζI[P],k2[P]] - EllipticF[ζr[P,r],k2[P]]

这基本上就是我想用 Python 做的事情。我尝试了很多东西，比如函数 scipy.special.ellipj，但它返回一个列表，我想要一个数字。请问有人知道怎么解决这个问题吗？

解决方法

与 SciPy 中 Mathematica 的 EllipticF 等效的是 scipy.special.ellipkinc。

例如，在 Wolfram Alpha 中，EllipticF(3/4,1/10) 给出数值 0.7564213553633020688976920...。

使用 SciPy，

In [49]: from scipy.special import ellipkinc                                                                                               

In [50]: ellipkinc(0.75,0.1)                                                                                                              
Out[50]: 0.7564213553633021

（我不知道您传递给 EllipticF 的参数是什么，所以我无法帮助您。）