coefa_fixdia function

Fix(Replace) the diagonal values in the matrix