Least Squares using Cholesky Decomposition

This function computes the (linear) least squares solution using Cholesky decomposition

This function computes the n-dimensional column vector x that minimizes norm(b-A*x),
where A is an m-by-n coefficient matrix and b is the m-dimensional right side column vector ( m >> n ).