sn2(12z,k)=1-cn(z,k)1+dn(z,k)=1-dn(z,k)k2(1+cn(z,k))=dn(z,k)-k2cn(z,k)-k2k2(dn(z,k)-cn(z,k)),