function [h] = detLaplace(A)

[n,~]=size(A);

if n==1
    h=A;
else
    h=0;
    for i=1:n

        Aij=A([1:i-1,i+1:n],2:n);
        h=h+(-1)^(i+1)*detLaplace(Aij)*A(i,1);

    end
end