Я знаю, как решить в Mathematica с помощью EllipticF
:
EllipticF[ζI[P], k2[P]] - EllipticF[ζr[P, r], k2[P]]
Это в основном то, что я хочу сделать в Python. Я пробовал много вещей, таких как функция scipy.special.ellipj
, но она возвращает список, и мне нужен номер. Пожалуйста, кто-нибудь знает, как решить эту проблему?
Может быть, этот блог-пост Джона Кука пригодится здесь?
Эквивалентом EllipticF Mathematica в SciPy является 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
, поэтому я не могу вам с ними помочь.)
Пожалуйста, предоставьте минимальный воспроизводимый пример, а также текущий и ожидаемый результат.