
m=0.5;
K = ellipke(m);
x=linspace(0,4*K,100);
[SN,CN,DN] = ellipj(x,m);
%AM=jacobiAM(x,m);
plot(x,SN,'k',x,CN,'--k',x,DN,':k')
legend('sn','cn','dn','location','southwest')
 axis([0,4*K,-1,1])

 
[SN2,CN2,DN2] = ellipj(x+K,m);
figure(2)
plot(x,SN2,'k',x,CN/DN,'--k')