[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
besseli: Am I missing something?
From: |
Robert T. Short |
Subject: |
besseli: Am I missing something? |
Date: |
Wed, 11 Jan 2012 20:52:00 -0800 |
User-agent: |
Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.1.4) Gecko/20091017 SeaMonkey/2.0 |
Before I file a bug report, maybe I should check that I am not all wet.
If I recall my Bessel function theory,
besseli(n,x) = besseli(-n,x)
if n is an integer.
In octave 3.4.2 I get
octave:3> besseli(1,1),besseli(-1,1)
ans = 0.565159103992485
ans = 0.565159103992485
Looks good for n=1, but for n=10
octave:4> besseli(10,1),besseli(-10,1)
ans = 2.75294803983687e-10
ans = -1.40610343427994e-07
Not so good!
And for n=100
octave:5> besseli(100,1),besseli(-100,1)
ans = 8.47367400813812e-189
ans = 7.38028373423800e+170
BTW, the numbers for positive n agree to about 9 or 10 decimal places
with my tables. Also, BTW, the relationships for besselj seem to work fine.
Anybody?
Bob
- besseli: Am I missing something?,
Robert T. Short <=
- Re: besseli: Am I missing something?, Jordi GutiƩrrez Hermoso, 2012/01/12
- Re: besseli: Am I missing something?, Michael D Godfrey, 2012/01/12
- Re: besseli: Am I missing something?, Jussi Lehtola, 2012/01/12
- Re: besseli: Am I missing something?, John W. Eaton, 2012/01/12
- Re: besseli: Am I missing something?, Dr. Alexander Klein, 2012/01/12
- Re: besseli: Am I missing something?, Jussi Lehtola, 2012/01/12
- Re: besseli: Am I missing something?, Thomas Weber, 2012/01/12
Re: besseli: Am I missing something?, Michael D Godfrey, 2012/01/12