I remember seeing a closed form for the series :
∞∑k=0Γ(k+z)ψ(k+1)k!xk|x|<1
But i don't seem to recall what it was . I tried replacing the gamma function with its Euler's integral , and then changing the order of summation and integration, using the fact that :
∞∑k=0ψ(k+1)k!yk=ey(logy+Γ(0,y)
Thus :
∞∑k=0Γ(k+z)ψ(k+1)k!xk=∫∞0yz−1e(x−1)y[(log(xy)+Γ(0,xy))]dy
=x−z∫∞0ωz−1e(1−1x)ω[logω+Γ(0,ω)]dω
Where Γ(⋅,⋅) is the incomplete gamma function.
But i don't know how to do the integral !
Answer
Lemma 1. By the extended binomial theorem, ∑k≥0Γ(k+z)k!xk=Γ(z)(1−x)z.
Lemma 2. The series definition of the digamma function ensures
ψ(k+1)=−γ+∑n≥1(1n−1n+k)=−γ+∑n≥1∫+∞0e−nt(1−e−kt)dt=−γ+∫+∞01−e−ktet−1dt=−γ+∫10uk−1u−1du.
By Lemma 1 and Lemma 2,
∑k≥0Γ(k+z)ψ(k+1)k!xk=−γΓ(z)(1−x)z+Γ(z)∫101u−1(1(1−ux)z−1(1−x)z)du
and by applying integration by parts
∑k≥0Γ(k+z)ψ(k+1)k!xk=−γΓ(z)(1−x)z−Γ(z+1)∫10xlog(1−u)(1−ux)z+1du.
Now the last integral can be seen as the derivative of a hypergeometric function (by substituting u↦(1−u) then exploiting log(u)=ddαuα|α=0+) or expanded as a series by exploiting log(1−u)=−∑m≥1umm.
No comments:
Post a Comment