integrate: (x^2 + 1)^k dx