First time here? Check out the FAQ!

Ask Your Question
0

Arbitrary precision BesselY

asked 12 years ago

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

Preview: (hide)

2 Answers

Sort by » oldest newest most voted
1

answered 12 years ago

achrzesz gravatar image

Try:

sage: import mpmath  
sage: mpmath.bessely?

(mpmath is multiprecision package)

Preview: (hide)
link

Comments

Also note that when http://trac.sagemath.org/sage_trac/ticket/4102 is finished, mpmath will be the default evaluation mode and so asking for as much precision as you want *should* work.

kcrisman gravatar imagekcrisman ( 12 years ago )
0

answered 12 years ago

FrankSt gravatar image

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

Many thanks Frank

Preview: (hide)
link

Your Answer

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

Add Answer

Question Tools

Stats

Asked: 12 years ago

Seen: 322 times

Last updated: Jan 03 '13