function c=cofactor(a)
% cofactor matrix
if length(a)==2,
    c=[a(2,2) -a(1,2)
       -a(2,1) a(1,1)]
else
       c=inv(a)*det(a);
end