Let $n\in \mathbb N$ be a natural number, $x_1,\cdots.x_n$ be formal variables. Consider the following $n\times n$-matrix $M_n:=diag\{x_1^2+\cdots+x_n^2,\cdots,x_1^2+\cdots+x_n^2\}$, can we find a solution $A=(a_{ij})$ such that
- (Square root): $A^2=M_n$.
- (Polynomial): $a_{ij}\in \mathbb F[x_1,\cdots,x_n]$ be a polynomial, where $\mathbb F$ is a number field such as $\mathbb Q, \mathbb R$ or $\mathbb C$.
Of course, we do not order $A$ is diagonal. Moreover, we especially care the case of $n=3$, since $n=1,2$ are relatively simple. For example $$ \begin{pmatrix} x_1^2+x_2^2 &\\ &x_1^2+x_2^2 \end{pmatrix}=\begin{pmatrix}\pm x_1 & x_2\\ x_2& \mp x_1\end{pmatrix}^2 .$$