Using properties of determinants, prove that : $ \begin{bmatrix} a^2+1 & ab & ac \\ ab & b^2+1 & bc \\ ca & cb & c^2+1 \end{bmatrix} = (1+a^2+b^2+c^2) $

Question

This question has appeared in textbook.