Shevchenko / 2sem / LR7_1

clc
clear
function [DET] = det_Gauss(A, N) %Задание функции
  for K=1:N-1 %Проходим по этементам строки (1)
    for I=K+1:N %Проходим по элементам столбца начиная со 2 строки (2)
      for J=K:N %Проходим по каждому элементу строки (3), начиная с элемента (1)
        if J==K %Если положение в строке элемента (3) равно положению (1)
          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
  S=A %вывод преобразованной матрицы
end

n=5;
a=[1 2 3 4 5; 10 9 8 7 6; 5 9 11 12 13; 20 1 3 17 14; 12 10 4 16 15]
tic
da=det_Gauss(a, n)
toc