Re: Hypergeometric1F1 polynomial

*To*: mathgroup at smc.vnet.net*Subject*: [mg91492] Re: [mg91437] Hypergeometric1F1 polynomial*From*: Devendra Kapadia <dkapadia at wolfram.com>*Date*: Sat, 23 Aug 2008 01:44:25 -0400 (EDT)*References*: <200808210956.FAA10365@smc.vnet.net>

On Thu, 21 Aug 2008, Alec Mihailovs wrote: > Mathematica gives the wrong answer to the following sum, > > In[1]:= Sum[Binomial[n, k]/Binomial[2 n, k]/k! (2 x)^k, {k, 0, n}] > > Out[1]= 2^(-(1/2) - n) E^x x^(1/2 + n) > BesselI[1/2 (-1 - 2 n), x] Gamma[1/2 - n] > > The correct answer is 1 for n=0 and Hypergeometric1F1[-n, -2 n, 2 x] for > integer n>0, which would be equal to the expression given by Mathematica if > n was not a positive integer. > > Another form of the correct answer is > > (2 x)^(n+1/2) E^x BesselK[n+1/2,x] n!/(2 n)!/Sqrt[Pi] > > Is there a way to apply some assumptions to get the correct answer? > > Alec > Hello Alec, Thank you for reporting the problem with the answer returned by Sum in the above example. As noted by you, this sum can be evaluated in terms of Hypergeometric1F1. A partial workaround for the problem is to introduce a parameter 'a' as shown in In[2] below, to obtain a Hypergeometric1F1 function. The incorrect behavior in your example is caused by auto-simplification of this hypergeometric function (see In[4] below). ========================= In[2]:= Sum[(Binomial[n, k]/Binomial[a*n, k]/k!)*(2*x)^k, {k, 0, n}] Out[2]= Hypergeometric1F1[-n, -(a n), 2 x] In[3]:= Table[% /. {a -> 2}, {n, 0, 3}] // InputForm Out[3]//InputForm= {1, 1 + x, 1 + x + x^2/3, 1 + x + (2*x^2)/5 + x^3/15} In[4]:= %% /. {a -> 2} // InputForm Out[4]//InputForm= 2^(-1/2 - n)*E^x*x^(1/2 + n)*BesselI[(-1 - 2*n)/2, x]*Gamma[1/2 - n] ========================== I apologize for the confusion caused by this problem. Sincerely, Devendra Kapadia, Wolfram Research, Inc.

**References**:**Hypergeometric1F1 polynomial***From:*"Alec Mihailovs" <alec@mihailovs.com>