Integrate[Sqrt[(x^2-a)^2+b],x]