+
And maximize this volume.

Input := 


sols = Solve[rv'[a]==0,a]//N

Output =


{{a -> 1.47247}, {a -> 4.52753}}