clc; clear
%% Dane wejsciowe
A1 = [1,2;1.1,2];
B1 = [10;10.4];
A2 = [2,5.999999;2,6];
B2 = [8.000001;8];
A3 = [3,-1,-1;1,-3,-2;-1,2,1];
B3 = [4;-1;-1];

MojGauss(A1,B1)
MojGauss(A2,B2)
MojGauss(A3,B3)

%% Wyznaczanie wskaźników uwarunkowania
wsk_1=cond(A1);
wsk_2=cond(A2);
wsk_3=cond(A3);

%% Rozwiązywanie układów
function x = MojGauss(A,B)
    n = size(B,1);
    x = zeros(n,1);
    AB = [A B];
    for i=1:n
        AB(i,:)=AB(i,:)./AB(i,i);
        for j=1:n
            if i ~= j
                AB(j,:) = AB(j,:).*AB(j,i);
            end
        end
    end
    x = AB(:,n+1);
end