コーシー積は以下の公式で求められます。
sum(A[n]*x^n,n,0,inf)*sum(B[n]*x^n,n,0,inf)=sum(sum(A[j]*B[k-j]*x^k,j,0,k),k,0,inf);
今回は冪級数の2乗ですから$A_n=B_n$として、$A_n$にポッホハマー記号を用いた定義式を与えます。そのために下記の定義式を思い出します。 $${}_{2}F_1\left(a,b;a+b+\frac12;x\right)=\sum_{n=0}^{\infty}{\frac{(a)_n\,(b)_n\,x^n}{(a+b+\frac12)_n\,n!}}$$
A[n]:=pochhammer(a,n)*pochhammer(b,n)/(pochhammer(a+b+1/2,n)*n!);
$C_n$としてコーシー積の公式の右辺の$n$次の係数を計算します。
C[n]:=sum(A[j]*A[n-j],j,0,n);
実際に$C_0, C_1, C_2, C_3$を計算してみます。
CC:[C[0],C[1],C[2],C[3]];
同様に$${}_3F_2\left(2\,a,2\,b,a+b;2\,a+2\,b,a+b+\frac12;x\right)=\sum_{n=0}^{\infty}{\frac{(2\,a)_n\,(2\,b)_n\,(a+b)_n\,x^n}{(2\,a+2\,b)_n\,(a+b+\frac12)_n\,n!}}$$を思い出します。$D_n$をこの冪級数の$n$次の係数とします。
D[n]:=pochhammer(2*a,n)*pochhammer(2*b,n)*pochhammer(a+b,n)/(pochhammer(2*a+2*b,n)*pochhammer(a+b+1/2,n)*n!);
実際に$D_0, D_1, D_2, D_3$を計算してみます。
DD:[D[0],D[1],D[2],D[3]];
$CC=[C_0, C_1, C_2, C_3]$から$DD=[D_0, D_1, D_2, D_3]$を引き算します。リストの引き算では独立に要素ごとに計算されます。
CC-DD,ratsimp;
これで定理の両辺のマクローリン展開の$0\sim 3$次の係数が等しいことがわかり、証明は終了しました。