integrate 1/sqrt((x-r_1)(x-r_2)(x-r_3)(x-r_4))