Soient \(\vx_1,\ldots, \vx_n\) les colonnes de \(A^{-1}\text{.}\) On a alors
\begin{equation*}
A\vx_j = \ve_j.
\end{equation*}
Chaque colonne est donc solution d’une équation qu’on peut résoudre avec la méthode de Cramer. Si on écrit
\begin{equation*}
\vx_j = \left[ \begin{array}{r}x_{1,j}\\ x_{2,j}\\ \vdots \\ x_{n,j} \end{array}\right],
\end{equation*}
on a, d’après la règle de Cramer,
\begin{equation*}
x_{i,j} = \frac{|A_i(\ve_j)|}{|A|}\text{.}
\end{equation*}
Or \(A_i(\ve_j)\) est obtenue de \(A\) en remplaçant sa \(i^{\text{e}}\) colonne par \(\ve_j\text{,}\) et donc dans la \(i^{\text{e}}\) colonne il y a une unique entrée qui est non nulle:
\begin{equation*}
A_i(\ve_j) = \left[\begin{array}{ccccccc}
a_{1,1} \amp \cdots \amp a_{1, i-1} \amp 0 \amp a_{1,i+1}\amp \cdots \amp a_{1,n}\\
a_{2,1} \amp \cdots \amp a_{2, i-1} \amp 0 \amp a_{2,i+1}\amp \cdots \amp a_{2,n}\\
\vdots \amp \amp \amp \vdots \amp \amp \amp \vdots \\
a_{j-1,1} \amp \cdots \amp a_{j-1, i-1} \amp 0 \amp a_{j-1,i+1}\amp \cdots \amp a_{j-1,n}\\
a_{j,1} \amp \cdots \amp a_{j, i-1} \amp 1 \amp a_{j,i+1}\amp \cdots \amp a_{j,n}\\
a_{j+1,1} \amp \cdots \amp a_{j+1, i-1} \amp 0 \amp a_{j+1,i+1}\amp \cdots \amp a_{j+1,n}\\
\vdots \amp \amp \amp \vdots \amp \amp \amp \vdots \\
a_{n,1} \amp \cdots \amp a_{n, i-1} \amp 0 \amp a_{n,i+1}\amp \cdots \amp a_{n,n}\\
\end{array} \right].
\end{equation*}
Le calcul de \(|A_i(\ve_j)|\) le long de la \(i^{\text{e}}\) colonne donne un seul terme. On remarque que la sous-matrice qui importe est précisément \(A_{j,i}\text{,}\) de sorte que
\begin{equation*}
|A_i(\ve_j)| = (-1)^{i+j}|A_{j,i}|
\end{equation*}
qui est le coefficient en position \((i,j)\) de l’adjointe de \(A\text{.}\)