integrate 1/sqrt((x-p1)*(x-p2)*(x-p3)*(x-p4))