Group Abstract Group Abstract

Message Boards Message Boards

1
|
5.7K Views
|
4 Replies
|
1 Total Like
View groups...
Share
Share this post:

Obtain series expansions using Frobenius Method

Posted 6 years ago
POSTED BY: Brad Klee
4 Replies
Posted 6 years ago

Hi Daniel,

I think this looks okay. The other "hacker" way to do the calculation is:

Normal@Series[ x Exp[   FullSimplify@CoefficientList[
       Normal@  Series[- 2 Pi  /Sqrt[3] Hypergeometric2F1[1/3, 2/3, 1, 1 - x]
           /Hypergeometric2F1[1/3, 2/3, 1, x], {x, 0, 10}], Log[x]][[1]]], {x, 0, 3}]

Outputs:

Out[]:=(1/27)*x + (5 /243)*x^2 + (31/2187)* x^3

And the integerization is as above. There is another question as to whether the output should be a series with rational coefficients or lots of Gamma, PolyGamma whatever. I guess it doesn't matter if FullSimplify then works, nice. You are the expert, so I guess I will just leave it up to you from here, thanks though. This is much better than the original call.

--Brad

POSTED BY: Brad Klee

Thank you for pointing out this weak spot. I am now making some changes to better support these exponentials. So in a future release the first example will give the result below.

x^(2/Sqrt[3])*SeriesData[x, 0, {2^(-8/Sqrt[3]),
    1/(2^(8/Sqrt[3])*Sqrt[3]), 
   (2^(-5 - 8/Sqrt[3])*(16 + 13*Sqrt[3]))/3, 
   (2^(-5 - 8/Sqrt[3])*(117 + 85*Sqrt[3]))/27, 
   (2^(-14 - 8/Sqrt[3])*(49544 + 34293*Sqrt[3]))/27, 
   (2^(-14 - 8/Sqrt[3])*(629835 + 424286*Sqrt[3]))/405}, 0, 6, 1]
POSTED BY: Daniel Lichtblau
Posted 6 years ago

Hi Daniel,

Oh no! There was a copy-paste error in the first expansion, it was supposed to be:

Series[Exp[-Pi Hypergeometric2F1[1/2, 1/2, 1, 1 - x]/Hypergeometric2F1[1/2, 1/2, 1, x]], {x, 0, 5}]

and you can test the following by series expansion:

Exp[-Pi Hypergeometric2F1[1/2, 1/2, 1, 1 - x]/Hypergeometric2F1[1/2, 1/2, 1, x]] == EllipticNomeQ[x]

So that first one actually seems to work okay. The other three expansions are more of a concern, they come from Ramanujan's theory, and appear as examples in Bruce Berndt's "Ramanujan's Notebooks Vol. II" ( on p.80-82). Does your fix cover those examples as well?

I may be around the office tomorrow afternoon if you want a better explanation, or to see a few more examples I've cooked up.

--Brad

POSTED BY: Brad Klee

All four of the original examples will give more "expanded" expansions in future. Here is the second one, pared back for brevity.

Series[Exp[-2 Pi/Sqrt[3] Hypergeometric2F1[1/3, 2/3, 1, 1 - x]/Hypergeo
metric2F1[1/3, 2/3, 1, x]], {x, 0, 1}] // InputForm                             

(* Out[20]//InputForm= 
x^((2*Pi)/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]))*
 SeriesData[x, 0, {E^((4*EulerGamma*Pi)/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]) + 
     (2*Pi*PolyGamma[0, 1/3])/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]) + 
     (2*Pi*PolyGamma[0, 2/3])/(Sqrt[3]*Gamma[1/3]*Gamma[2/3])), 
   (-4*E^((4*EulerGamma*Pi)/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]) + 
       (2*Pi*PolyGamma[0, 1/3])/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]) + 
       (2*Pi*PolyGamma[0, 2/3])/(Sqrt[3]*Gamma[1/3]*Gamma[2/3]))*Pi*
     (2 + PolyGamma[0, 1/3] + PolyGamma[0, 2/3] - PolyGamma[0, 4/3] - 
      PolyGamma[0, 5/3]))/(9*Sqrt[3]*Gamma[1/3]*Gamma[2/3])}, 0, 2, 1] *)
POSTED BY: Daniel Lichtblau
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard