expectation of y^2+2y-1 for y having a gamma distribution