Shevchenko / 2sem / LR7_2

clc
clear
function [DET] = det_Gauss(A,N)
  for K=1:(N-1)
    for I=(K+1):N
      for J=K:N
        if J==K
          aik=A(I,K);
        else
          A(I,J)=A(I,J)-(aik/A(K,K)*A(K,J));
        end
      end
    end
  end
  DET=1;
  for i=1:N
    DET=DET*A(i,i);
  end
end

N=200;
A=rand(N);
tic
da=det_Gauss(A,N)
toc