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);