For $\left(U_1, V_1\right) \sim \mathcal{N}(0,0,1,1, \rho)$. In a standard bivariate normal with zero means and unit variances, the joint probability of both variables being positive equals:
$$ P\left[U_1>0, V_1>0\right]=\frac{1}{4}+\frac{\arcsin (\rho)}{2 \pi} $$
How to show the above result?