MN_new

A=[6,3,6,30;2,3,3,17;1,1,1,11];                         %deklaracja macierzy dla metody eliminacji Gaussa-Jordana
[m,n]=size(A);                                          %przypisanie do zmiennej m i n wymiarow macierzy A
%-> przekształcenia macierzy aby otrzymac wartośći rózne od 0 na jej przekatnej i w ostatniej kolumnie
for jj=1:1:m-1
  for ii=2:m
    if A(jj,jj)==0                                      %Jezeli na przekątnej macierzy głównej znajduje sie wartośc 0:
        A([1 ii],:) = A([ii 1],:)                       %zamiana wiersza pierwszego z i-tym
    end
  end
for k=jj+1:1:m
  A(k,:)=A(k,:)-A(jj,:)*(A(k,jj)/A(jj,jj));
end
end
for jj=m:-1:2
  for ii=jj-1:-1:1
    A(ii,:)=A(ii,:)-A(jj,:)*(A(ii,jj)/A(jj,jj));
end
end
%<-
for ii=1:1:m
  %podzielenie macierzy przez stale, tak aby otrzymac minor jednostkowy (po wykresleniu ostatniej kolumny)
  A(ii,:)=A(ii,:)/A(ii,ii);
  x(ii)=A(ii,n);                                        %przypisanie wartosci i-tego wiersza z ostatniej kolumny do zmiennej przechowującej wynik
end
disp('Wynik po obliczeniu metodą Gaussa-Jordana');
x'