cn⁡(z,k)=(1+k2′)⁢(dn2⁡(z/(1+k2′),k2)-k2′)k22⁢dn⁡(z/(1+k2′),k2),