This program solves the linear system AX=B by applying Cholesky factorization to the matrix A, which must be square (M-by-M) and Hermitian positive definite.