Re: Bessel K expansion, large argument?
- To: mathgroup at smc.vnet.net
- Subject: [mg70106] Re: Bessel K expansion, large argument?
- From: dimmechan at yahoo.com
- Date: Wed, 4 Oct 2006 05:58:41 -0400 (EDT)
- References: <eftdhr$6b9$1@smc.vnet.net>
Perhaps,
Normal[Series[z*(BesselK[1, z]/BesselK[0, z]), {z, Infinity, 6}]]
Series::esss :...(3 times)
General::stop : Further output of Series::esss will be suppressed
during this calculation.
1/2 + 103/(32*z^6) - 1073/(1024*z^5) + 13/(32*z^4) - 25/(128*z^3) +
1/(8*z^2) - 1/(8*z) + z
List @@ %
{1/2, 103/(32*z^6), -(1073/(1024*z^5)), 13/(32*z^4), -(25/(128*z^3)),
1/(8*z^2), -(1/(8*z)), z}
%[[{1, 8, 7, 6, 5, 4, 3, 2}]]
{1/2, z, -(1/(8*z)), 1/(8*z^2), -(25/(128*z^3)), 13/(32*z^4),
-(1073/(1024*z^5)), 103/(32*z^6)}
(HoldForm[Plus[##1]] & ) @@ %
HoldForm[1/2 + z - 1/(8*z) + 1/(8*z^2) - 25/(128*z^3) + 13/(32*z^4) -
1073/(1024*z^5) + 103/(32*z^6)]
is what you want.