3. Gauss z częściowym wyborem elementu głównego.
A = orgA;
B = orgB;

najwieksza_wartosc = max(A,[],2);
najwiekszy_x = 0;
if (najwieksza_wartosc(1) >= najwieksza_wartosc(2)) && (najwieksza_wartosc(1) >= najwieksza_wartosc(3)) && (najwieksza_wartosc(1) >= najwieksza_wartosc(4))
    najwiekszy_x = 1;
elseif (najwieksza_wartosc(2) >= najwieksza_wartosc(1)) && (najwieksza_wartosc(2) >= najwieksza_wartosc(3)) && (najwieksza_wartosc(2) >= najwieksza_wartosc(4))
    najwiekszy_x = 2;
elseif (najwieksza_wartosc(3) >= najwieksza_wartosc(1)) && (najwieksza_wartosc(3) >= najwieksza_wartosc(2)) && (najwieksza_wartosc(3) >= najwieksza_wartosc(4))
    najwiekszy_x = 3;
elseif (najwieksza_wartosc(4) >= najwieksza_wartosc(1)) && (najwieksza_wartosc(4) >= najwieksza_wartosc(2)) && (najwieksza_wartosc(4) >= najwieksza_wartosc(3))
    najwiekszy_x = 4;
end

temp1_najwiekszy_a = A(najwiekszy_x,1:end);
temp2_najwiekszy_a = A(1,1:end);

A(najwiekszy_x,1:end) = temp2_najwiekszy_a;
A(1,1:end) = temp1_najwiekszy_a;

temp1_najwiekszy_b = B(najwiekszy_x,1:end);
temp2_najwiekszy_b = B(1,1:end);

B(najwiekszy_x,1:end) = temp2_najwiekszy_b;
B(1,1:end) = temp1_najwiekszy_b;

i = 1;
while (i <= length(A)-1)
    j = i +1;
    while (j <= length(A))
        z=A(j,i)/A(i,i);
        A(j,i:end)=A(j,i:end)-(A(i,i:end)*(A(j,i)/A(i,i)));
        B(j)=B(j)-B(i)*z;
        j = j + 1;
    end
    i = i + 1;
end
disp('Ax=b')
x=inv(A)*B;
disp(x);