On 7/1/2010 4:45 PM, Jaasiel Ornelas wrote:
I'm working on modeling drums with bessel functions. However, I need
to calculate the zeroes of the bessel functions. Is there any way to
do that?
I'm not sure if this is the best way, but you could try, e.g.,
sage: import scipy.special
sage: [x for x in dir(scipy.special) if 'zero' in x]
['ai_zeros', 'bei_zeros', 'beip_zeros', 'ber_zeros', 'berp_zeros',
'bi_zeros', 'erf_zeros', 'flatnonzero', 'fresnel_zeros',
'fresnelc_zeros', 'fresnels_zeros', 'jn_zeros', 'jnjnp_zeros',
'jnp_zeros', 'jnyn_zeros', 'kei_zeros', 'keip_zeros', 'kelvin_zeros',
'ker_zeros', 'kerp_zeros', 'nonzero', 'trim_zeros', 'y0_zeros',
'y1_zeros', 'y1p_zeros', 'yn_zeros', 'ynp_zeros', 'zeros', 'zeros_like']
sage: scipy.special.jn_zeros?
sage: scipy.special.jn_zeros(int(0), int(3)).tolist()
[2.4048255576957724, 5.5200781102863106, 8.6537279129110125]
List of scipy.special functions:
http://www.scipy.org/SciPyPackages/Special
Possibly naive questions:
* Are there analogous arbitrary-precision routines in or wrapped by Sage?
* Is it possible to call these GSL functions easily from Sage:
http://www.gnu.org/software/gsl/manual/html_node/Zeros-of-Regular-Bessel-Functions.html
?
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
