plot (-((Sqrt[1 - a])/a) + ArcSin[Sqrt[a]]/a^(3/2))/2 for 0 < a < 1