Ask Your Question

Arbitrary precision BesselY

asked 2013-01-03 02:52:57 -0500

FrankSt gravatar image

Hi, I consider using Sage for some calculations which contain Bessel functions of complex arguments. Since I have to mix Bessel functions with very small and very large arguments I require precision higher than 15 digits.

Unfortunately, I recognized that BesselY is not implemented in PARI. But BesselJ and BesselH1 are. My physicists understanding of math tells me, I could just use "(BesselH1-BesselJ)/i". But I am surprised that this has not been discussed before (at least I couldn't find it), since this would allow a quick implementation of BesselY. Am I missing something that's obvious for math experts? Or can I just use above definition to get an arbitrary precision BesselY?

Many thanks Frank

edit retag flag offensive close merge delete

2 answers

Sort by ยป oldest newest most voted

answered 2013-01-03 03:16:23 -0500

achrzesz gravatar image


sage: import mpmath  
sage: mpmath.bessely?

(mpmath is multiprecision package)

edit flag offensive delete link more


Also note that when is finished, mpmath will be the default evaluation mode and so asking for as much precision as you want *should* work.

kcrisman gravatar imagekcrisman ( 2013-01-03 03:24:51 -0500 )edit

answered 2013-01-03 03:39:18 -0500

FrankSt gravatar image

Didn't know about mpmath. It seems to provide all I need.

Many thanks Frank

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools


Asked: 2013-01-03 02:52:57 -0500

Seen: 53 times

Last updated: Jan 03 '13