sum(x^k*gamma(y+k),k=0..n)