Weierstrass's elliptic functions implementation

[GoogleCodeExporter]

Weierstrass's elliptic functions implementation request in terms of elliptic 
integrals and jacobi's elliptic functions

http://en.wikipedia.org/wiki/Weierstrass's_elliptic_functions
http://www.mai.liu.se/~halun/complex/elliptic/
http://mathworld.wolfram.com/WeierstrassEllipticFunction.html

Original issue reported on code.google.com by moiseev....@gmail.com on 10 Jan 2010 at 4:51

[moiseevigor]

As proposed on gitter https://gitter.im/moiseevigor/elliptic?at=5773e91e971717155488a036 it seems that it can be done with this function

function P = WeiP(z,e1,e2,e3)
  m = (e2 - e3)/(e1 - e3);
  w = z*sqrt(e1 - e3);
  sn = ellipj(w,m);
  P = e3 + (e1 - e3)/(sn^2);
end