integrate sqrt(x^2+p*x+q)dx