RE: BesselK[1/2,x]
- To: mathgroup
- Subject: RE: BesselK[1/2,x]
- From: Mike Monagan <mbmonagan at watcauchy.waterloo.edu>
- Date: Wed, 22 Mar 89 14:23:16 CST
> >From Jim_Wendel at ub.cc.umich.edu Tue Mar 21 06:19:52 1989
> Subject: Re: Chebyshev "bug"
> ...
> Plot[BesselK[1/2,x],{x,20,25}] came out beautifully, values from 0 down to
> -0.0000175 on the vertical axis.
=================================================
Except that BesselK[1/2,x] is positive for x > 0.
Abramowitz & Stegun, pp. 444, give a formula for
BesselK[1/2,x] in terms of Exp[-x] namely
10.2.17. sqrt(Pi/2/x) K(1/2,x) = Pi/2/x exp(-x)
from which you can deduce this.
I tried the Plot on our Mac II (version 1.1 of Mathematica)
and it gives random noise beyond x = 21.
Mike
mbmonagan at watcauchy.waterloo.edu