Le cofacteur d'indice (j,k) d'une matrice carrée A est (-1)^{k+j}\det(A_{kj}) où A_{kj} désigne la matrice qu'on obtient en enlevant de A la k-ième ligne et la j-ième colonne. Autrement dit, si A est de format nxn alors A_{kj} est la matrice suivante de format (n-1)x(n-1)

A_{kj}= \begin{pmatrix}a_{1,1} & \dots & a_{1,j-1}& a_{1,j+1}& \dots & a_{1,n} \\\vdots & & \vdots & \vdots& &\vdots\\ a_{k-1,1} & \dots & a_{k-1,j-1}& a_{k-1,j+1}& \dots & a_{k-1,n} \\ a_{k+1,1} & \dots & a_{k+1,j-1}& a_{k+1,j+1}& \dots & a_{k+1,n} \\ \vdots & & \vdots & \vdots &&\vdots\\ a_{n,1} & \dots & a_{n,j-1}& a_{n,j+1}& \dots & a_{n,n}\end{pmatrix}\;.

La matrice des cofacteurs de A, s'appelle la comatrice de A, notée com(A). En résumé,

\text{com}(A) = \left((-1)^{k+j}\det(A_{kj})\right)_{1\leq k,j\leq n}

Petit exercice :  la fonction qui à une matrice associe sa comatrice est-elle un difféomorphisme du groupe linéaire GL(n,\mathbb{R}) sur lui-même ? Et de GL(n,\mathbb{C}) sur lui-même ?