W1 = ellperiods([1,I]);
W2 = ellperiods([1,exp(2/3*Pi*I)]);
ellwp(W1)
deriv(ellwp(W1))
ellwp(W1,1/2,1)
ellwp(W1,I/2,1)
ellwp(W1,1/2+I/2,1)
ellwp(W2)
deriv(ellwp(W2))
ellwp(W2,1/2,1)
ellwp(W2,exp(2/3*Pi*I)/2,1)
ellwp(W2,1/2+exp(2/3*Pi*I)/2,1)