%Wyznacznik
%Funckja tworzy macierz trojkatna, a nastepnie liczy wyznaczik, ktory
%jest iloczynem elementów na przekątnej


function wyzn = wyznacznik(B)
A = B;


n_temp = size(A);
n = n_temp(1);

x=zeros(n,1);

for k = 1:n
    for i = k:n
        if i+1 <= n

            if (A(i+1, k) ~= 0)
                temp = A(i+1, k) / A(k,k);
                A(i+1, k:n) = A(i+1, k:n) - temp .* A(k, k:n);
            end
        end
    end
end

wyzn = 1;
for i = 1:n
    wyzn = wyzn*A(i,i);
end