cn⁡(z,k)=cn⁡(z/(1+k1),k1)⁢dn⁡(z/(1+k1),k1)1+k1⁢sn2⁡(z/(1+k1),k1),