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
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
Соседние файлы в папке 2sem