$$
Sigma = begin{pmatrix}
            sigma_x^2 & rho sigma_x sigma_y \
            rho sigma_x sigma_y & sigma_y^2
         end{pmatrix}
$$

M_inverse = M.inv()

M = Matrix(2,3, [1,2,3,4,5,6])