matrix_estimation function

Matrix Estimation