# Problems with series of generalized hypergeometric functions

I have been trying to do the following series:

Series[HypergeometricPFQ[{1, 4, 4}, {4 - Sqrt[3], 4 + Sqrt[3]}, z],{z,1,0}]


Mathematica says that the result is infinite, without telling me precisely how it diverges with 1-z. I have tried FunctionExpand, FullSimplify and other commands, but so far, with no success. Any help would be most welcome.

Best, Jorge

Probably because it has a rather complicated branch point at z=1. It is actually very easy to derive the explicit expansion:

singular[x_] = Csc[Sqrt[3] Pi] (-I Sqrt[3] Pi^2 + Pi (1 + 1/(Sqrt[3] (-1 + x)))
- 2 Sqrt[3] Pi HarmonicNumber[Sqrt[3]] - Sqrt[3] Pi Log[-1 + x]);

finite[x_] = (3 - 2 Sqrt[3]) HypergeometricPFQ[{1, 4, Sqrt[3]},
{2 + Sqrt[3], 4 + Sqrt[3]}, 1];

suppressed[x_] = (x - 1) 16/ 13 (-69 + 40 Sqrt[3])
HypergeometricPFQ[{2, 5, Sqrt[3]}, {3 + Sqrt[3], 5 + Sqrt[3]}, 1]
+ (x - 1) 1/2 Pi Csc[Sqrt[3] Pi] (-10 - 9 Sqrt[3] + 10 I Sqrt[3] Pi
+ 20 Sqrt[3] HarmonicNumber[Sqrt[3]] + 10 Sqrt[3] Log[-1 + x]);

asymp[x_] = singular[x] + finite[x] + suppressed[x] (* + O((x-1)^2 Log[x-1]) *);


Notice the logarithmic function in the singular part. Here is the comparison with the original function:

y = 10^-8;
N[HypergeometricPFQ[{4, 4, 1}, {4 + Sqrt[3], 4 - Sqrt[3]}, 1 + y] - asymp[1 + y], 16]

(* -1.892181059711013*10^(-13) + 3.39982418261034*10^(-13) I *)


## Derivation

I derived this result by using the integral representation of a generalized hypergeometric function:

integral = HypergeometricPFQ[{a_, b_, c_}, {d_, e_}, x_] :>
Gamma[e]/(Gamma[c] Gamma[e - c]) z^(c - 1) (1 - z)^(e - c - 1) Hypergeometric2F1[a, b, d, x z];


You can easily check it:

Integrate[ Gamma[e]/(Gamma[c] Gamma[e - c]) z^(c - 1) (1 - z)^(e - c - 1)
Hypergeometric2F1[a, b, d, x z], {z, 0, 1}, Assumptions -> Re[c] < Re[e] && Re[c] > 0]

(* HypergeometricPFQ[{a, b, c}, {d, e}, x] *)


Secondly, I applied a simple transformation rule for the integrand:

hyperTransform = HoldPattern[Hypergeometric2F1[a_, b_, c_, z_]] :>
(Gamma[c] Gamma[c - a - b])/(Gamma[c - a] Gamma[c - b])
Hypergeometric2F1[a, b, a + b - c + 1, 1 - z]
+ (1 - z)^(c - a - b) (Gamma[c] Gamma[a + b - c])/(Gamma[a] Gamma[b])
Hypergeometric2F1[c - a, c - b, c - a - b + 1, 1 - z];


Thus, I obtained the following integrand:

Collect[HypergeometricPFQ[{4, 4, 1}, {4 + Sqrt[3], 4 - Sqrt[3]}, x]
/. integral /. hyperTransform, _Hypergeometric2F1, Composition[PowerExpand, FullSimplify]]


which is the sum of two terms:

-Pi x^(-3 + Sqrt[3]) (1 - z)^(-1 + Sqrt[3]) z^Sqrt[3] (1 - x z)^(-1 - Sqrt[3]) Csc[Sqrt[3] Pi]
- (3 + 2 Sqrt[3]) (1 - z)^(-1 + Sqrt[3]) z^3 Hypergeometric2F1[1, 4, 2 + Sqrt[3], 1 - x z]


For the second part, the series expansion about x=1 commutes with the integration with respect to z, hence you can expand it and then integrate, e.g, the leading term is (I made a substitution z-> 1-z)

Integrate[ -(3 + 2 Sqrt[3]) z^(-1 + Sqrt[3]) (1 - z)^3 Hypergeometric2F1[1, 4, 2 + Sqrt[3], z],
{z, 0, 1}]

(* (3 - 2 Sqrt[3]) HypergeometricPFQ[{1, 4, Sqrt[3]}, {2 + Sqrt[3], 4 + Sqrt[3]}, 1] *)


The first term you need to integrate first and then expand:

Integrate[-Pi x^(-3 + Sqrt[3]) (1 - z)^(-1 + Sqrt[3]) z^Sqrt[3]
(1 - x z)^(-1 - Sqrt[3]) Csc[Sqrt[3] Pi], {z, 0, 1}, Assumptions -> Im[x] != 0]

(* -((4^-Sqrt[3] Pi^(3/2) x^(-3 + Sqrt[3]) Csc[Sqrt[3] Pi] Gamma[Sqrt[3]]
Hypergeometric2F1[1 + Sqrt[3], 1 + Sqrt[3], 1 + 2 Sqrt[3], x])/Gamma[1/2 + Sqrt[3]]) *)

FullSimplify[Normal@Series[%, {x, 1, 0}], Assumptions -> x > 1 && Arg[x] == 0]


Thus, you can find any term of the expansion.