numSingular = 0;
numTests = 50;
matSize = 3;

for i = 1:numTests
  A = rand(matSize);
  if (abs(det(A)) < 10^(-6))
    numSingular = numSingular + 1;
  end
end

numSingular