Re: extracting fractional powers of series expansion?

*To*: mathgroup at smc.vnet.net*Subject*: [mg78294] Re: [mg78225] extracting fractional powers of series expansion?*From*: DrMajorBob <drmajorbob at bigfoot.com>*Date*: Wed, 27 Jun 2007 05:41:31 -0400 (EDT)*References*: <26733444.1182852184761.JavaMail.root@m35>*Reply-to*: drmajorbob at bigfoot.com

For example: Series[Sqrt[Sin[x]], {x, 0, 10}] % // FullForm Sqrt[x]-x^(5/2)/12+x^(9/2)/1440-x^(13/2)/24192-(67 x^(17/2))/29030400+O[x]^(21/2) SeriesData[x,0,List[1,0,0,0,Rational[-1,12],0,0,0,Rational[1,1440],0,0,0,Rational[-1,24192],0,0,0,Rational[-67,29030400]],1,21,2] The first argument to SeriesData is the variable (x in this case) that's raised to powers in the sum. The last argument, 2, is the denominator of all the powers... so we're really talking about powers of Sqrt[x]. Those powers range from 1 to 21 (the fourth and fifth arguments). The third argument is the list of coefficients. You can puzzle things out from there, if necessary, but this also works: SeriesCoefficient[Sqrt[Sin[x]], {x, 0, 21/2}] -1/5677056 Here's a more complicated example: Series[x!, {x, Infinity, 5}] % // FullForm E^SeriesData[$CellContext`x, Infinity, {-1 - Log[$CellContext`x^(-1)]}, -1, 6, 1]*SeriesData[$CellContext`x, Infinity, {Sqrt[2*Pi], 0, Rational[1, 6]*Sqrt[Rational[1, 2]*Pi], 0, Rational[1, 144]*Sqrt[Rational[1, 2]*Pi], 0, Rational[-139, 25920]*Sqrt[Rational[1, 2]*Pi], 0, Rational[-571, 1244160]*Sqrt[Rational[1, 2]*Pi], 0, Rational[163879, 104509440]*Sqrt[Rational[1, 2]*Pi]}, -1, 11, 2] Times[Power[E,SeriesData[x,DirectedInfinity[1],List[Plus[-1,Times[-1,Log[Power[x,-1]]]]],-1,6,1]],SeriesData[x,DirectedInfinity[1],List[Power[Times[2,Pi],Rational[1,2]],0,Times[Rational[1,6],Power[Times[Rational[1,2],Pi],Rational[1,2]]],0,Times[Rational[1,144],Power[Times[Rational[1,2],Pi],Rational[1,2]]],0,Times[Rational[-139,25920],Power[Times[Rational[1,2],Pi],Rational[1,2]]],0,Times[Rational[-571,1244160],Power[Times[Rational[1,2],Pi],Rational[1,2]]],0,Times[Rational[163879,104509440],Power[Times[Rational[1,2],Pi],Rational[1,2]]]],-1,11,2]] The first argument of Times is, of course, a function to be multiplied by the second argument: E^SeriesData[x, Infinity, {-1 - Log[x^(-1)]}, -1, 6, 1] \[ExponentialE]^(-1-Log[1/x])/(1/x)+O[1/x]^6 The second argument of Times is a series with fractional powers, much like the first example, but now the simple extraction method doesn't work, since SeriesCoefficient[x!, {x, Infinity, 9/2}] returns unevaluated. I'm not sure there's any method that works in every case, but you can see, I hope, how you'd have to proceed in an individual case. Good luck! Bobby On Tue, 26 Jun 2007 03:32:16 -0500, Tanim Islam <tanim.islam at gmail.com> wrote: > Hi: > > How do I extract a given fractional power term in a series expansion? For > instance, I do a series expansion of a function about infinity. How do I > extract the -1/2 power term with, say, SeriesCoefficient. I am using > Mathematica 6. > > Tanim Islam > > -- DrMajorBob at bigfoot.com