function d= mydet(a)
% Функция для вычисления определителя
% с помощью рекурсии
[m,n]= size(a);
if m~=n
    error('Матрица должна быть квадратной!');
end 
if m==1
    d= a;
else
    d= 0;
    for i= 1:m
        min= a; min(1,:)= []; min(:,i)= [];
        d= d+a(1,i)*(-1)^(1+i)*mydet(min);
    end
end
