Integrate[1/Sqrt[(x - r[1])*(x - r[2])*(x - r[3])*(x - r[4])],x]