integral x^n/(pi * gamma * (1 + ((x-m)/gamma)^2)) dx