$$M=\begin{pmatrix} 1+x_1^2 & x_1x_2 &...&x_1x_n \\ x_2x_1 & 1+x_2^2 &...&x_2x_n \\...&...& &...& \\x_nx_1 & x_nx_2& ...&1+x_n^2&\end{pmatrix}.$$
So I noticed that $M$ is a symmetric matrix and $ M=I+(x_1, x_2,...,x_n)^T(x_1, x_2,...,x_n)$. That's all I can get.
You can continue from your observation using the matrix-determinant lemma (which is not hard to prove, by the way).
For a column vector $x = (x_1,...,x_n)$, and an $n\times n $ identity matrix $I$, we have
$$ M = I + x x^T $$
hence, in view of the aforementioned lemma,
$$ \det M = (1 + x^T I x)\det I = 1 + \|x\|^2 = 1 + x_1^2 + \cdots +x_n^2. $$