Re: roots of BesselJPrimes
- To: mathgroup at smc.vnet.net
- Subject: [mg14165] Re: [mg14152] roots of BesselJPrimes
- From: BobHanlon at aol.com
- Date: Wed, 30 Sep 1998 02:04:14 -0400
- Sender: owner-wri-mathgroup at wolfram.com
Christian,
$Version
"Power Macintosh 3.0 (May 6, 1997)"
Needs["NumericalMath`BesselZeros`"]
BesselJPrimeZeros[1, {1, 20}]
{1.84118,5.33144,8.53632,11.706,14.8636,18.0155,21.1644,24.3113,27.4571,
30.6019,33.7462,36.89,40.0334,43.1766,46.3196,49.4624,52.605,55.7476,58.89,
62.0323}
Alternatively, if your implementation does not work with BesselJZeros:
g[z_] := Evaluate[D[BesselJ[1, z], z]]; z /. Table[FindRoot[g[z] == 0,
{z, n Pi}], {n, 20}]
{1.84118,5.33144,8.53632,11.706,14.8636,18.0155,21.1644,24.3113,27.4571,
30.6019,33.7462,36.89,40.0334,43.1766,46.3196,49.4624,52.605,55.7476,58.89,
62.0323}
Bob Hanlon
In a message dated 9/28/98 11:31:26 PM, meier at tomo.uni-bremen.de wrote:
>I'm trying to find the first 20 roots of the derivation of a
>Besselfunction of the first kind,
>first order. The BesselJPrimeZeros-function doesn't work while for
>example
>BesselYPrimeZeros does.
>I'm working with mathematica3.0 on a sparc 10 computer.
>
>Example for in- and outputs:
>
>In[31]:=
>BesselYPrimeZeros[1,{1,3}]
>
>Out[31]=
>{3.68302,6.9415,10.1234}
>
>but:
>
>In[32]:=
>BesselJPrimeZeros[1,{1,3}]
>
>Out[32]=
>BesselJPrimeZeros[1,{1,3}]