How to solve Jacobian elliptical function in python

Question

I know how to solve in Mathematica with "EllipticF": EllipticF[ζI[P], k2[P]] - EllipticF[ζr[P, r], k2[P]]. That\'s basically what i want to do in python. I\'ve tri