function adjA=mnhf_adjoint(A)
%MNHF_ADJOINT calculate adjoint of a (square) matrix
%
% A - the matrix in question
%
% Usage  mnhf_adjoint(A);

[r,c] = size(A);  % Compute the size of A.
if isequal(r,c)
 if r==1 
  fprintf('Matrix has dimension 1-by-1.')
 else
  
  % Initialize array.
  adjA = zeros(r,c);
  
  % Evaluate minors.  
  for ii=1:r
   for jj=1:c   
    % Split out cases where we know the row (column) number is 1 or r 
    % (c). Use : notation to facilitate matrix index referencing.
    if ii==1
     if jj==1
      adjA(ii,jj) = mnhf_determinant(A(2:r,2:c));
     elseif (jj>1 && jj<r)
      adjA(ii,jj) = mnhf_determinant(A(2:r,[1:jj-1,jj+1:c]));
     else
      adjA(ii,jj) = mnhf_determinant(A(2:r,1:c-1)); 
     end
    elseif (ii>1 && ii<r)
     if jj==1      
      adjA(ii,jj) = mnhf_determinant(A([1:ii-1,ii+1:r],2:c));
     elseif (jj>1 && jj<r)
      adjA(ii,jj) = mnhf_determinant(A([1:ii-1,ii+1:r],[1:jj-1,jj+1:c]));
     else
      adjA(ii,jj) = mnhf_determinant(A([1:ii-1,ii+1:r],1:c-1)); 
     end
    else
     if jj==1
      adjA(ii,jj) = mnhf_determinant(A(1:r-1,2:c));
     elseif (jj>1 && jj<r)      
      adjA(ii,jj) = mnhf_determinant(A(1:r-1,[1:jj-1,jj+1:c]));
     else
      adjA(ii,jj) = mnhf_determinant(A(1:r-1,1:c-1)); 
     end
    end
   end
  end
  
  % Turn minors into cofactors then compute matrix transpose. 
  for ii=1:r
   for jj=1:c
    adjA(ii,jj) = (-1.0)^(ii+jj)*adjA(ii,jj);
   end
  end
  adjA = adjA';
  
 end
else
 fprintf('Can''t evaluate adjoint of non-square A.\n\n')
end