int (x^2 + 4x + 1)/((x + 1)*(x - 1)*(x + 3)) dx