Soit $S \in \mathcal{S}_{n}^{+}(\mathbb{R})$
Montrons que : $\exists ! R \in \mathcal{S}_n^{+}(\mathbb{R}), R^{2} = S$

Existence

$\newline$
$S$ est symétrique réelle donc orthodiagonalisable par le théorème spectral pour les endomorphismes auto-adjoints. Ils existent donc $P \in \mathcal{O}_{n}(\mathbb{R})$ orthogonale et $D = diag(\lambda_{i})_{i\in [\![ 1,n ]\!]}$ diagonale dont ses coefficients sont les valeurs propres de $S$ telles que $S = P^{T} D P = P^{-1}DP$.
$\newline$
Comme $S \in \mathcal{S}_{n}^{+}(\mathbb{R})$, le spectre de $S$ est positif. On peut donc considérer la matrice diagonale $\Delta$ définie par $\Delta = diag(\sqrt{\lambda_{i}})_{i \in [\![1,n]\!]}$.
$\newline$
On pose ensuite $R = P^{T} \Delta P$
$\newline$
Cette matrice est bien une racine carrée de $S$ dans $\mathcal{S}_{n}^{+}(\mathbb{R})$ :

• $R^{2} = (P^{-1} \Delta P)^{2} = P^{-1} \Delta^{2} P = P^{-1}DP = S$
• $R^{T} = (P^{T} \Delta P)^{T} = P^{T} \Delta^{T} (P^{T})^{T} = P^{T} \Delta P = R$, donc $R \in \mathcal{S}_{n}(\mathbb{R})$
• $R$ est à spectre positif (son spectre est $(\sqrt{\lambda_{i}})_{i \in [\![1,n]\!]}$) donc $R \in \mathcal{S}_{n}^{+}(\mathbb{R})$
  $\newline$

Unicité

$\newline$
Soit $R' \in \mathcal{S}_{n}^{+}(\mathbb{R})$ vérifiant $R'^{2} = S$. Montrons que $R = R'$ où $R$ est la même matrice que celle de la phase "Existence" de la preuve. Ce qui permettra de dire que $R$ est la seule racine carrée de $S$ dans $\mathcal{S}_{n}^{+}(\mathbb{R})$.
$\newline$
Par interpolation de Lagrange, il existe un polynôme $Q$ vérifiant : $\forall i \in [\![1,n]\!], Q(\lambda_{i}) = \sqrt{\lambda_{i}}$
Nous avons ainsi la relation polynomiale entre $R$ et $S$ suivante :
$Q(S) = Q(P^{-1} D P) = P^{-1} Q(D)P = P^{-1} diag(Q(\lambda_{i})) P = P^{-1} diag(\sqrt{\lambda_i}) P = P^{-1} \Delta P= R$
$\newline$
La relation $R'^{2} = S$ nous dit que $S$ est un polynôme en $R'$, par conséquent $R'$ commute avec $S$ et plus généralement avec tout polynôme en $S$ donc, en particulier, avec $R$.
$\newline$
De plus, $R$ et $R'$ sont des matrices symétriques réelles donc sont diagonalisables par le théorème spectral.
$\newline$
Donc, comme $R$ et $R'$ sont diagonalisables et commutent, nous pouvons alors les diagonaliser simultanément. Il existe alors une matrice $\tilde{P} \in GL_{n}(\mathbb{R})$ telle que $R = \tilde{P}^{-1} \Delta \tilde{P}$ et $R' = \tilde{P}^{-1} \Delta' \tilde{P}$ où $\Delta'$ est une matrice diagonale contenant les valeurs propres de $R'$.
$\newline$
En partant de l’hypothèse de départ faite sur $R'$ : $S = R^{2} = R'^{2}$
On a : $\tilde{P}^{-1} \Delta^{2} \tilde{P} = \tilde{P}^{-1} (\Delta') ^{2} \tilde{P}$ puis en simplifiant par $\tilde{P}$ et $\tilde{P}^{-1}$ : $\Delta^{2} = (\Delta')^{2}$.
$\newline$
En notant par $(\mu_{i})_{i\in [\![1,n]\!]}$ le spectre de $R'$, la relation $\Delta^{2} = (\Delta')^{2}$ nous dit que $\forall i \in [\![1,n]\!], \lambda_i = \mu_{i} ^2$.
Comme le spectre de $R'$ est positif car $R' \in \mathcal{S}_{n}^{+}(\mathbb{R})$, on a $\forall i \in [\![1,n]\!], \mu_{i} = \sqrt{\lambda_{i}}$.
Ce qui permet de conclure en affirmant que $\Delta = \Delta'$ et ainsi $\fbox{R = R'}$.