使用matlab编写矩阵行列式

function f= finddet(A)
[m,n]=size(A);
if m==1
    f=A(1);
else
    cur=0;
    for i=1:n
        B=[];
        for j=1:n
            if j==i
                continue;
            else
                B=[B,A(2:n,j)];
            end
        end
        cur=cur+A(1,i)*(-1)^(i+1)*finddet(B);
    end
    f=cur;
end
end