Numeric mathematics bessel-functions « Next Oldest | Next Newest »

 ▼ peacecalc Member Posts: 97 Threads: 9 Joined: Nov 2011 06-30-2012, 02:43 PM Hello numeric fans, does anybody knows a numerical function, which calculates every root of the bessel function first kind. I have only an algorithm for the hp 50g in USER-RPL which is able to found a root with input an estimation of the position for the root. It would be more usefull, when one can input let's say: 1. until 5. root and/or when it is possible to find the last root. Greetings peacecalc ▼ Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 06-30-2012, 03:02 PM Hello, Interesting that you pose this question. I am writing an article for the online HP Solve newsletter that discusses a mutli-root finder that locates roots in a range of values. The method can locate roots that are also minima or maxima (something that most root-seeking algorithms fail to do). Namir ▼ peacecalc Member Posts: 97 Threads: 9 Joined: Nov 2011 06-30-2012, 05:32 PM Hello Namir, When your article will be published? I can't wait reading it! It exists an integral representation for first kind functions, the integrand function can be used to find roots of the first derivate, that's a way to find the last root of the bessel function itself. Greetings peacecalc C.Ret Senior Member Posts: 260 Threads: 0 Joined: Oct 2008 07-01-2012, 01:21 AM The functions Bessel function Jn(z) each have an infinite number of real zeros, all of which are simple with the possible exception of z = 0. For nonnegative n, the kth positive zeros of first kind Bessel functions are denoted jn,k, except that z = 0 is typically counted as the first zero of J'0(z)(Abramowitz and Stegun 1972, p. 370). The first few roots jn,k of the Bessel functions are given in the following table for small nonnegative integer values of n and k. ```k J0(x) J1(x) J2(x) J3(x) J4(x) J5(x) 1 2.4048 3.8317 5.1356 6.3802 7.5883 8.7715 2 5.5201 7.0156 8.4172 9.7610 11.0647 12.3386 3 8.6537 10.1735 11.6198 13.0152 14.3725 15.7002 4 11.7915 13.3237 14.7960 16.2235 17.6160 18.9801 5 14.9309 16.4706 17.9598 19.4094 20.8269 22.2178 ``` Sources : http://mathworld.wolfram.com/BesselFunctionZeros.html http://en.wikipedia.org/wiki/Bessel_function ▼ peacecalc Member Posts: 97 Threads: 9 Joined: Nov 2011 07-01-2012, 02:44 AM Hello C.Ret, thank you for your private lesson in mathematics. I forgot the important fact you mentioned. I threw some things to disarray. The approximation for great values of z is written with a cos function, which has obviously an infinit number of roots. But one problem is left, how to find every roots in a given range. I'm waiting for Namir's article, the suspense is killing me. Greetings peacecalc ▼ C.Ret Senior Member Posts: 260 Threads: 0 Joined: Oct 2008 07-01-2012, 02:57 AM Yuo are wel come. I am also waiting on Namir's article. This lesson is only a copy-paste from InterNet sources about Bessel functions. But it was a good (or needed) refresch for my memories on the Bessel topic too! I am glad to see that sharing this can also help anyone else ! :-) As you, I was a bit disapointed that no algorithm or 'approximation function or code' was available nor advertise on the Web. ▼ Ángel Martin Posting Freak Posts: 1,253 Threads: 117 Joined: Nov 2005 07-01-2012, 03:30 AM I´d suggest you check Jean-Marc Baillard´s authoritative pages. He´s got a few programs worth studying, like the Mac Mahon expansion one, see it here: It is also included in the BESSEL ROM, availabe at TOS - and by extension on the CL of course; where speed is definitely not an issue. ▼ C.Ret Senior Member Posts: 260 Threads: 0 Joined: Oct 2008 07-01-2012, 05:42 AM Thanks a lot. This Mac Mahon expansion expression si exactly the 'code' I was looking for ! peacecalc Member Posts: 97 Threads: 9 Joined: Nov 2011 07-01-2012, 10:06 AM Hello Martin, it looks like a taylor expansion of the asymptotic approximation, is that right? Greetings peacecalc Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 07-01-2012, 04:43 PM Balliard's use of the MacMahon approximation seems to work for roots that are high in their sequence number. Edited: 1 July 2012, 4:44 p.m. ▼ Ángel Martin Posting Freak Posts: 1,253 Threads: 117 Joined: Nov 2005 07-02-2012, 02:28 AM Quote: seems to work for roots that are high in their sequence number. Can you give an example where it doesn´t work? ▼ Jean-Michel Member Posts: 191 Threads: 41 Joined: Jun 2007 07-02-2012, 07:43 AM Hi, MacMahon's formulae are very good for large roots. I've used the first 5 terms of these asymptotic expansions, but better formulas are given in the "Digital Library of Mathematical Functions" up to the terms of order 7 cf http://dlmf.nist.gov/10.21 For the first roots, one can use a root-finding program and ascending series or continued fractions. Unfortunately, I don't know direct formulas in these cases. Hope this helps, JM. ▼ peacecalc Member Posts: 97 Threads: 9 Joined: Nov 2011 07-05-2012, 01:37 PM Hello JM at al., thank you for interesting links and further information. If I understand right, the formulae from Mac Mahon works even for real k parameter values. For integer values is an easy approach possible (and x --> oo). again thank you and Greetings peacecalc

 Possibly Related Threads... Thread Author Replies Views Last Post HP Prime suggestion to avoid Numeric/Symbolic confusion Chris Pem10 4 850 11-19-2013, 05:49 AM Last Post: bluesun08 [41CL] New Extra Functions version Monte Dalrymple 0 515 11-08-2013, 04:32 PM Last Post: Monte Dalrymple Prime as laptop numeric keypad Jose Gonzalez Divasson 1 544 10-31-2013, 02:25 AM Last Post: debrouxl HP Prime: in need of help with defining functions Alberto Candel 14 1,922 10-27-2013, 10:48 AM Last Post: Alberto Candel Using units in Numeric Solver Harold A Climer 1 539 10-13-2013, 10:44 AM Last Post: Tim Wessman HP Prime spreadsheet functions SanS 0 1,058 10-04-2013, 04:23 AM Last Post: SanS Stats functions on the HP34S Nicholas van Stigt 5 901 09-24-2013, 02:45 AM Last Post: Nick_S Trig Functions Howard Owen 11 1,568 09-16-2013, 02:53 PM Last Post: Fred Lusk 50g piecewise functions Kurt Fankhauser 6 993 09-15-2013, 08:01 PM Last Post: Kurt Fankhauser Missing functions on the HP Prime!!!??? :-( Namir 6 976 08-22-2013, 08:40 AM Last Post: Gilles Carpentier

Forum Jump: