Numerical evaluation of HypergeometricPFQ[{a,b,c},{d,e},1]
- To: mathgroup at christensen.cybernetics.net
- Subject: [mg497] Numerical evaluation of HypergeometricPFQ[{a,b,c},{d,e},1]
- From: "NELSON M. BLACHMAN" <blachman at gtewd.mtv.gtegsc.com>
- Date: Wed, 01 Mar 1995 22:28:48 PST
This is a postscript to my 26 February message (reproduced below)
complaining about Mma's not giving a numerical value for
HypergeometricPFQ[{1/2,1,1},{1.45,1.55},1].
I've since found that Mma does evaluate things like
HypergeometricPFQ[{1/2,1,1},{1.4,1.6},0.99999], but it takes
something like an hour on my 486DX33 PC--and it takes longer
and longer as the last argument gets closer and closer to 1.
So it's good that Mma quickly announces its inability to
compute HypergeometricPFQ[{1/2,1,1},{1.4,1.6},1].
A different method is evidently needed when the last argument is
1--provided that the sum of the components of the middle argument
exceeds the sum of the components of the first argument. (If not,
the hypergeometric function seems likely to be infinite.) So I
continue to hope for a formula for 3F2[{a,b,c},{d,e},1] as a
ratio of products of gamma functions of (e, f, and e + f minus
0, a, b, c, a + b, etc.), though I've so far been unable to
devise a satisfactory conjecture of this sort.
Nelson M. Blachman
>From: GTEWD::BLACHMAN "NELSON M. BLACHMAN" 26-FEB-1995 23:54:14.57
>To: MX%"mathgroup at christensen.cybernetics.net"
>CC: BLACHMAN
>Subj: Numerical Evaluation of HypergeometricPFQ
I was pleased to see just now that Mma's able to evaluate
Integrate[Sin[t u]/Sin[t],{t,0,Pi/2}] in terms of HypergeomtricPFQ.
When I asked it to plot the result, however, I found it apparently
unable to determine numerical values for HypergeomtricPFQ.
Maple can sometimes compute numerical values for HypergeomtricPFQ,
but for the particular denominator indices here it complains of
iteration limits' being exceeded. Does anyone know of a way to
get Mma 2.2 to evaluate HypergeomtricPFQ numerically?
I suspect there's a simple expression for HypergeomtricPFQ[{ },{ },1]
in terms of gamma functions when P = Q + 1; there is, anyhow, if
P = 2 and Q = 1. Does anyone know if that's true and, if so, what
it is?
Nelson M. Blachman
GTE Government Systems Corp.
Mountain View, California
Mathematica 2.2 for DOS 387
Copyright 1988-93 Wolfram Research, Inc.
In[1]:= Integrate[Sin[t u]/Sin[t],{t,0,Pi/2}]
3 u Pi u 1 3 - u 3 + u
(- + -) u Cos[----] HypergeometricPFQ[{-, 1, 1}, {-----, -----}, 1]
2 2 2 2 2 2
Out[1]= -------------------------------------------------------------------
3 u 1 u 1 u
4 (-(-) - -) (-(-) - -) (- - -)
2 2 2 2 2 2
In[2]:= f[t_]:= %1 /. u -> t Simplify ELIMINATES THE FIRST FACTOR IN
THE NUMERATOR AND IN THE DENOMINATOR.
In[3]:= f[.1] // N
Out[3]= 0.0997665 HypergeometricPFQ[{0.5, 1., 1.}, {1.45, 1.55}, 1.]