function Z=adjoint(A)
[a,b]=size(A);
for m=1:a
    for n=1:b
        c=1;
        for k=1:a
            if k ~= m
                B(c,:)=[A(k,1:(n-1)), A(k,(n+1):b)];
                c=c+1;
            end
        end
        X(m,n)=((-1)^(m+n))*det(B);
    end
end
Z=transpose(X);
return;