integrate dx/sqrt((x-a)^2+(x-b)^2)