Pivoted Cholesky decomposition for low-rank approximation

Question

I am trying to implement a pivoted Cholesky decomposition in MatLab, which computes a low-rank approximation for a symmetric positive definite matrix (with exponentially dec