EllepticPi argumet is complex. so can not plot it. How to handle this problem?

inttau[r_]=-(1/Sqrt[0.0345106153943703 - ((37.3042 - r) (-25.578 + r) (62.8822 +   r))/(3000 r)]) 

This is my function of r, now I integrated it w r t r

tauanalytical[r_] = Integrate[inttau[r], r] 

This is what I got after integration, I got this result.

-((2.3094 Sqrt[1. + 62.1447/r] (-31.0723 + r) EllipticPi[3.,ArcSin[0.57735 Sqrt[1. + 62.1447/r]], 1.])/Sqrt[-0.965489 + 20./r + 0.000333333 r^2]) 

later I tried to find the numerical value of;

tauanalytical[30] // N  

I got this, no the exact number , which I was expecting

-126.491 EllipticPi[3., -1.5708 + 0.153761 I, 1.] 

Here the second argument of Ellepticpi comes out to be imaginary , it should be real . Please provide necessary assistance.