acX_to_rho function

Computes the known correlation $\rho$between $\widehat{\theta}$ and $\widehat{\tau}$ from the $p$-vectors $a$ and $c$ and the design matrix $X$