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