% wygenerowanie macierzy do rozkładu
A = randi([-2 2], 6)

% algorytm el. Gaussa
[lw, lk] = size(A);
for k = 1 : (lk-1)
  for w = (k+1) : lw
    A(w,:) = A(w,:) - A(w,k)/A(k,k)*A(k,:);
  end
end

% część Jordana
for k = lk:-1:2
  for w = (k-1) : -1 : 1
    A(w,:) = A(w,:) - A(w,k)/A(k,k)*A(k,:);
  end
end
A


  
