Inverse of jacobian elliptic function

Hello! I am using the complete elliptic integral of first kind K(k) depending on k.

syms m %m=k^2

Kk=ellipticK(m);

I know that there isn't the inverse function in matlab. Does someone know how to calculate the inverse ( from K(k) to k)?

Thanks

Answer 1

[Barkan Ugurlu]

Hello Luca,

I have implemented the following method.

Say Wy0 = mu*jacobiSN(q0,m). Here, Wy0,mu and m are known, q0 is sought.

eqy = Wy0 == mu*jacobiSN(q0Sym, m);

q0 = double(vpasolve(eqy, q0Sym));

Then you can double check:

KakuninWy = mu*jacobiSN(q0,m) - Wy0

If KakuninWy is zero, then it should be fine.

Barkan