I have an equation that describes the exact solution to a non-linear 2nd order ODE that contains Guass' Hypergeometric2F1 function.
m := n + 1
z := c/c0
x := 1/Phi*Sqrt[m/2/c0^(m - 2)]*2/(2 - m)*(z^(1 - m/2)*Hypergeometric2F1[1/2 - 1/m, 1/2, 3/2 - 1/m, z^-m] - Hypergeometric2F1[1/2 - 1/m, 1/2, 3/2 - 1/m, 1])
I would like to generate a series expansion of that equation and then take the InverseSeries to generate the final solution. This is described in a publication by Eugen Magyari in some detail using Mathematica. I cannot reproduce his results despite talking to the author and getting his detailed Mathematica file. Our results match until I use the InverseSeries command where I get:
SeriesData::scmp: The constant term of series c^2 + O[c]^34 does not match the expansion point of {Series generated from above}
The author used Mathematica v4 and did not get this error. Is there a bug or has something changed with HyperGeometric2F1 or with InverseSeries that keeps it from generating a solution as it did in v4? I've attached the output file from v4.
Attachments: