matA = matrix(data=c(1,2,3,4,5,6),nrow=2,ncol=3)
xvec = c(1,2,3)
xvec.mat = as.matrix(xvec)
t(mat) t(xvec)
A %*% B
diag(n)
invA = solve(A)
\[A^{-1} = \frac{1}{det(A)} \begin{bmatrix} a_{22} & -a_{12} \\ -a_{21} & a_{11} \end{bmatrix}, ~~ \text{assuming: } det(A) = a_{11}a_{22} - a_{12}a_{21} \neq 0\]
\[f(x) = \frac{1}{2\pi det(\Sigma)^\frac 12} e^{- \frac 12 (x - \mu)^T \Sigma^{-1} (x-\mu)}\]