MathGroup Archive 2000

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Mod Bessel function bug ?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg23019] Re: [mg22988] Mod Bessel function bug ?
  • From: BobHanlon at aol.com
  • Date: Tue, 11 Apr 2000 23:18:36 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

Plot[Evaluate[Table[BesselK[n, x], {n, 0, 2, 1/3}]], {x, 0.1, 1.5}];

Series[BesselK[n, x], {x, 0, 5}] == Series[BesselK[-n, x], {x, 0, 5}] // 
    Normal // Simplify

True

FullSimplify[BesselK[n, x] == BesselK[-n, x]]

True

Abramowitz & Stegun 10.2.4 argues against your statement that "the I function 
with negative, noninteger order is the same as the K function"

Sqrt[1/2*Pi/z]*BesselK[n + 1/2, z] == 
    Pi/2 * (-1)^(n + 1)*
      Sqrt[1/2*Pi/z]*(BesselI[n + 1/2, z] - BesselI[-n - 1/2, z]);

Simplify[FunctionExpand[%], Element[n, Integers]]

Sqrt[Pi/2]*Sqrt[1/z]*(-BesselK[-(1/2) - n, z] + BesselK[1/2 + n, z]) == 0

FullSimplify[% /. n -> (m - 1/2)]

True

Bob Hanlon

In a message dated 4/10/2000 10:16:07 AM, jrchaff at nwlink.com writes:

>Well, something is very strange.
>
>Thank both of you for your replies.  I am using Mathematica
>Student Version 4.0; supposedly same as full version
>capabilities, or at least so advertised.
>
>My plots show there is a difference near zero;  however, both
>functions come together (and become large) for large argument,
>precisely opposite to what bessel theory says.  Supposedly,
>the I function with negative, noninteger order is the same as
>the K function (not with Mathematica), and the K function
>goes to zero exponentially with large argument.
>
>Something is definitely wrong.
>


  • Prev by Date: Mathematica Wavelet-Explorer Question: Inverse Wavelet Transform on a subset of the coefficients
  • Next by Date: Chop?
  • Previous by thread: Re: Mod Bessel function bug ?
  • Next by thread: Re: Mod Bessel function bug ?