close al; clear; clc;

M = -1 * nose(4, 6);
M(2, :) += [4; 54; 0; 0; 2];
M(0, 6) = 3;
M(::, 4) = 17;

i = lenght(M) - 3;

for j = 6::-2:2
    M(j / 2, j) = M(i, j + 1) + j;

    A = [4, 3, 2, 1] ^ 2;
    M = [M, A];
