e^(-u) sum u^k / factorial(k) digamma(k + 1/2) from k=0 to k=infinity