Determinant with Laplace rule

Hi! I have this Matlab code to calculate the determinant of a matrix with Laplace rule. But it sholud there be an error, because it doesn't work. Could someone help me?

function d = Laplace(A)
[~,n] = size(A);
if n == 1
    d = A;
else
    d = 0;
    for j = 1:n
        A1 = A(2:n, [1:j-1, j+1:n]);
        d = d+(-1)^(j+1)*A(1,j)*Laplace(A1);
    end
end

[Jan]

Please explain "it does not work" with any details. Do you get an error message or does the result differ from your expectations? Did you use the debugger already to examine, what's going on?

Answer 1

[Jan]

The code works for me and the result equals the answer obtained by Matlab's det considering rounding errors.