M = [2  3  1 -6;
     0 -5 -4  1;
     3  4 -1  0;
     0  4  5  2];
 
display(M);
sumOfDiagonals = trace(M) + trace(fliplr(M)); % сума діагональних елементів
sumOfAllElements = sum(sum(M)); % сума всіх елементів матриці
sumOfNonDiaglonal = sumOfAllElements - sumOfDiagonals;% сумма недіагональних елементів
for i = 1:length(M)
    if i ~= (length(M) - i + 1)
        M(i, i) = M(i, i) + sumOfNonDiaglonal;
        M(i, (length(M) - i + 1)) = M(i, (length(M) - i + 1)) + sumOfNonDiaglonal;
    else
        M(i, i) = M(i, i) + sumOfNonDiaglonal;
    end
end
display(M);