I would like to know how can we evaluate the Weierstrass $\wp$-functions. That is, I would like to find $\wp(\theta,\omega,i\omega)$ for some $\theta,\omega\in\mathbb{R}$.
I'm only able to find a function which outputs the Laurent series of the Weierstrass $\wp$-function when an elliptic curve has been entered. Should I evaluate that laurent series?
The function weierstrass_p() of EllipticCurve returns the Laurent expansion of $\wp$ at the origin, hence evaluating it only gives a reasonable approximation near it.
Nothing comes to my mind to do this kind of numerical evaluation straightforwardly in Sage. If you have access to Maple, you could give a shot at the NumGFun package, part of the AlgoLib library http://algo.inria.fr/libraries/. It has support for the numerical evaluation of functions satisfying linear differential equations with polynomial coefficients.
The author of NumGFun gave a talk on it at Sage Days 49 http://www.marc.mezzarobba.net/#expose-sd49. As you can read in the slides, nothing of it is already in Sage, but there are plans for the close future. See, for example, http://trac.sagemath.org/ticket/14996, which will add functionality similar to what you need.
The Jacobi elliptic functions in #14996 are already in Sage, the patch just improves them.
Nice, I didn't know that. Then you may try using the formulas given here: <http://en.wikipedia.org/wiki/Weierstrass%27s_elliptic_functions#Relation_to_Jacobi_elliptic_functions></p<>>
Asked: 2013-08-06 11:39:19 +0200
Seen: 814 times
Last updated: Aug 07 '13
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
