Prove that Laplace's equation $\Delta u = 0$ is rotation invariant; that is, if $O$ is an orthogonal $n\times n$ matrix and we define
$$ v(x) := u(Ox) \quad (x \in \mathbb{R}^n)$$ then $\Delta v = 0$.
I wanted to see if what I had below was correct and complete. Any feedback on rigor is greatly appreciated.
In cartesian coordinates, we have $\Delta u = \sum_{i = 1}^n \frac{\partial^2 u}{\partial x_i^2}$. Letting the $i$th row and $j$th column of the matrix $O$ be $a_{ij}$ we have
\begin{equation} \frac{\partial v}{\partial x_j} = \sum_{i =1}^n\frac{\partial u}{\partial x_i}a_{ij} \end{equation}
Thus, \begin{equation} \frac{\partial^2 v}{\partial x_j^2} = \sum_{i =1}^n\sum_{k = 1}^n\frac{\partial^2 u}{\partial x_i\partial x_k}a_{ij}a_{kj} \end{equation}
Therefore: \begin{align} \Delta v &= \sum_{j = 1}^n \frac{\partial^2 v}{\partial x_j^2}\\ &= \sum_{j = 1}^n \left(\sum_{i =1}^n\sum_{k = 1}^n\frac{\partial^2 u}{\partial x_i\partial x_k}a_{ij}a_{kj}\right)\\ &= \sum_{i =1}^n\sum_{k = 1}^n\left(\sum_{j = 1}^n \frac{\partial^2 u}{\partial x_i\partial x_k}a_{ij}a_{kj}\right)\\ &= \sum_{i =1}^n\sum_{k = 1}^n \frac{\partial^2 u}{\partial x_i\partial x_k}\left(\sum_{j = 1}^na_{ij}a_{kj}\right)\\ &=\sum_{i =1}^n\sum_{k = 1}^n \frac{\partial^2 u}{\partial x_i\partial x_k}\delta_{i, k}\\ &= \sum_{i =1}^n\frac{\partial^2 u}{\partial x_i^2} = 0 \end{align}