[sympy] Re: Bessel Function Zeros

2009-05-16 Thread Vinzent Steinberg
On May 15, 9:04 pm, Ondrej Certik ond...@certik.cz wrote: vinzent.steinb...@googlemail.com wrote: Please note that you also could use the brand-new nsolve() to do this without calling mpmath directly and using a simpler syntax. (If besselj () was implemented in sympy.) By the way,

[sympy] Re: Bessel Function Zeros

2009-05-16 Thread Ondrej Certik
On Sat, May 16, 2009 at 1:50 AM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: On May 15, 9:04 pm, Ondrej Certik ond...@certik.cz wrote: vinzent.steinb...@googlemail.com wrote: Please note that you also could use the brand-new nsolve() to do this without calling mpmath directly

[sympy] Re: Bessel Function Zeros

2009-05-15 Thread Ondrej Certik
On Thu, May 14, 2009 at 10:04 PM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Thu, May 14, 2009 at 04:10:47PM -0700, Ondrej Certik wrote: However, as far as I know, scipy doesn't do spherical bessel zeros, so of course I can use scipy's jn() to evaluate the function, but I still

[sympy] Re: Bessel Function Zeros

2009-05-15 Thread Gael Varoquaux
On Thu, May 14, 2009 at 11:35:00PM -0700, Ondrej Certik wrote: In [2]: jn_zeros(4, 5, method=scipy) Out[2]: [8.18256145257, 11.7049071546, 15.0396647076, 18.3012559595, 21.5254177334] E.g. it either uses pure sympy's jn() and findroot() (default), or scipy's sph_jn() and newton(). Scipy's

[sympy] Re: Bessel Function Zeros

2009-05-15 Thread Ondrej Certik
On Fri, May 15, 2009 at 3:57 AM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: On May 15, 12:47 am, Ondrej Certik ond...@certik.cz wrote: Hi, for testing that my FEM calculation of the radial Schroedinger equation is correct, I test the claculated eigenvalues. One example is a

[sympy] Re: Bessel Function Zeros

2009-05-14 Thread Ondrej Certik
On Thu, May 14, 2009 at 4:10 PM, Ondrej Certik ond...@certik.cz wrote: On Thu, May 14, 2009 at 3:53 PM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Thu, May 14, 2009 at 03:47:59PM -0700, Ondrej Certik wrote: Hi, for testing that my FEM calculation of the radial Schroedinger