Let $X,Y,Z$ have a standard trivariate normal distribution centered at the origin, with zero means, unit variances, and correlation coefficients $\rho_1,\rho_2,\rho_3$. Show that $$P(X>0,Y>0,Z>0)=\frac{1}{8}+\frac{1}{4\pi}(\sin^{-1}(\rho_1)+\sin^{-1}(\rho_2)+\sin^{-1}(\rho_3))$$
My idea is transform $X,Y,Z$ into i.i.d. normal. First, we have $V=\begin{bmatrix}1 & \rho_1& \rho_2\\ \rho_1 & 1 & \rho_3\\ \rho_2 & \rho_3 & 1\end{bmatrix}$. $V=P^TDP$ where columns of $P$ are eigenvectors and $D$ is a diagonal matrix. Then we have $V^{1/2}V^{1/2}=P^TD^{1/2}D^{1/2}P$. So, $V^{-1/2}=PD^{-1/2}$ since $P^TP=I$. Hence, $$\begin{bmatrix}U\\V\\W\end{bmatrix}=V^{-1/2}\begin{bmatrix}X\\Y\\Z\end{bmatrix}$$ I got stuck at this step, and I don't know how to find $V^{-1/2}$. Can someone give me a hint or suggestion? Thanks