MM = 10 + 8 * randn(5,7);
SR = mean(MM(:));
MM(MM > SR + 8) = SR + 8;
MM(MM < SR - 8) = SR - 8;
MMC = MM(:);
MMC = sort(MMC);
index = round(length(MMC) / 2);
median_val = MMC(index)
MM1 = log(MM)