Write $C=(c_{ij})$ and $(d_1,d_2,d_3)=(0,2,4)$. Because $D$ is a diagonal matrix,Thus $CD=DC$ exactly when $(d_j-d_i)c_{ij}=0$ for every $i,j$. The three eigenvalues are distinct, so every off-diagonal entry of $C$ must be zero; diagonal entries are unrestricted. Consequently the centralizer of a diagonal matrix here isEvery displayed matrix indeed commutes with $D$, since diagonal matrices multiply entry by entry on the diagonal.

For the given real matrix $M$, direct matrix multiplication gives the eigenvectorsPut these eigenvectors into the columns ofSince $\det P=-2\ne 0$, $P$ is an invertible matrix. The three column identities give $MP=PD$, soIn particular, the matrix similarity is established by an explicit change-of-basis matrix.

For any matrix $Y$ and compatible invertible matrix $P$, the linear map $r\mapsto rP$ maps $\mathrm{rowspace}(Y)$ onto $\mathrm{rowspace}(YP)$: a linear combination of rows of $Y$ is sent to the same linear combination of rows of $YP$. Its inverse function is $s\mapsto s{P}^{-1}$, so these row spaces have equal dimensions. Hence right matrix multiplication by $P$ also preserves row rank.

If $X=P^{-1}YP$, ii now givesThis proves the required matrix similarity invariant without needing to assume equality of row rank and column rank.

When $A$ is an invertible matrix, $B=A^{-1}(AB)$. Apply i first to $AB$ and then to $A^{-1}(AB)$:Thus the two row spaces are equal. In particular, invertible left matrix multiplication preserves row rank.

Write $A=(a_{ij})$ as an $m\times r$ matrix and $B$ as an $r\times n$ matrix, with rows $b_1,\dots ,b_r$. The rule for matrix multiplication givesEvery row of $AB$ is therefore a linear combination of rows of $B$. Taking linear spans gives

First let $b>0$. Substitute $x=-y$ in the negative half of the integral. Since $g$ is an odd function,The two halves cancel, soFor $b=0$ both limits coincide. For $b<0$, reverse the limits and apply the result with $-b>0$. Thus the identity holds for every real $b$.

Write $a=2\pi m+r$, where $m$ is an integer and $0\le r<2\pi$. The change of variables $x=y+2\pi m$ and the fact that $f$ is a periodic function giveSplit at $2\pi$, then substitute $y=z+2\pi$ in the second piece:ThereforeThis proves translation invariance of an integral over one period without assuming that $f$ is differentiable.

The core contributes no change in $\phi$, since ${\phi }^{\prime }=0$ there. Integrating the radial gravitational field across the shell, with $u=r/r_0$, givesUsing the supplied antiderivative,HenceThe negative sign agrees with ${\phi }^{\prime }(r)<0$ in the material shell under the paper's gravitational potential sign conventions.

On a sphere of radius $r$, the radial gravitational field has constant magnitude. Thus Gauss's law for gravity becomeswhere $M(r)$ is the enclosed mass. In the hollow core $M(r)=0$. In the material shell, put $u=r/r_0$ and use the spherical coordinates volume element:Outside the shell the enclosed mass is the totalConsequently,At $r=0$ the displayed value refers to the radial derivative of the constant core value. The expressions match at both shell boundaries. There is no surface sheet of mass, so a jump in the volume mass density does not produce a jump in the gravitational field. The core has zero gravitational field, and the exterior gravitational field equals that of a central point mass $M_{\ast }$, as in the shell theorem.

The mass density depends only on $r$, so it has spherical symmetry. For this isolated planet, choose the Newtonian gravitational potential to vanish at infinity, or equivalently choose $\phi =-\Phi$ to do so. There is no imposed external gravitational field. Rotating the mass distribution therefore leaves $\phi$ unchanged: $\phi =\phi (r)$.

The gradient in spherical coordinates consequently has no angular components, and the paper's gravitational potential sign conventions giveAt the centre there is no distinguished radial direction, and spherical symmetry makes $\mathbf{g}(0)=0$.

Integrate the steady state equation over $R$ and use the divergence theorem:The zero heat flux condition makes the boundary term vanish. But $\rho >0$ and $S>0$ throughout a region of positive volume, soa contradiction. Thus no steady-state solution exists. This is the insulated steady heat balance: positive internal heat generation continually increases the stored thermal energy when no heat can escape.

Suppose $T^{(1)}$ and $T^{(2)}$ are two sufficiently regular steady state solutions with the same prescribed temperature on the boundary. Their difference $w=T^{(1)}-T^{(2)}$ satisfiesThe identical PDE source terms cancel. Apply the product rule and the divergence theorem to $wk\nabla w$; this is the weighted form of Green's first identity:Since $k>0$, the nonnegative integrand must vanish, so $\nabla w=0$. Hence $w$ is constant on each connected component of $R$, and its zero Dirichlet boundary condition makes each constant zero. Therefore $\boxed{T^{(1)}=T^{(2)}}$: the solution, if it exists, is unique. This proves Dirichlet uniqueness for variable-conductivity heat conduction; it does not require a constant thermal conductivity.

The trigonometric identityshows that only angular Fourier modes $1$ and $3$ are required. Differentiating the harmonic expansion on a disk at $r=1$ multiplies the coefficient of each nonconstant Fourier mode by $n$. Thus $B_1=3/4$ and $3B_3=-1/4$, givingEach term is a regular harmonic function, and the radial derivative isThe Neumann compatibility condition for Laplace's equation holds becauseas is necessary by the divergence theorem applied to the gradient of a harmonic function.

The additive constant is undetermined. This Neumann boundary condition is unchanged by $T\mapsto T+C$, so the paper's blanket uniqueness assertion requires a normalization in this case. The uniqueness of the Neumann problem gives uniqueness up to a constant; prescribing $T(0)=0$, or zero mean temperature, selects the displayed answer with $C=0$.

The boundary data consist of the third angular Fourier mode. Its regular radial factor in the harmonic expansion on a disk is $r^3$, soFor a direct check of Laplace's equation,Also $T(1,\theta )=\sin (3\theta )$. At the origin the apparent coordinate singularity is harmless: this is the harmonic polynomial $3x^{2}y-y^3$. The allowed uniqueness for the Dirichlet problem completes the identification.

The constant function satisfies Laplace's equation and the prescribed Dirichlet boundary condition. In the harmonic expansion on a disk, its only nonzero coefficient is $A_0=1$. By the allowed uniqueness for the Dirichlet problem,

On the unit circle, the positively oriented unit tangent vector is $\mathbf{u}=(-y,x,0)$, and $|\mathbf{u}|=1$. Hence $d\mathbf{r}=\mathbf{u}ds$ andSince the scalar field $\psi$ has no $z$ dependence, the curl product rule for a scalar field givesOn the unit upper hemisphere $S_1$, the outward unit normal vector is $(x,y,z)$, so its vertical component is $z=\cos \theta$. Applying the Stokes theorem to the vector field $\psi \mathbf{u}$ therefore provesThe smoothness needed here is that $\psi$ extend as a continuously differentiable scalar field to a neighbourhood of the projected unit disk. The left side is a scalar line integral with respect to arc length; the positive boundary orientation is used when identifying it with the vector line integral.

Choose the outward orientation on the upper half of the ellipsoid, giving counterclockwise boundary orientation on $C$ as viewed from $z>0$. Parametrize the boundary circle byThen $\mathbf{u}(\mathbf{c}(t))=(-\sin t,\cos t,0)={\mathbf{c}}^{\prime }(t)$, so the line integral isFor a direct surface integral, use the parametrized surfaceIts outward vector surface element is supplied by the cross productUsing the curl from a,The two integrals agree for every $\alpha >0$, verifying the Stokes theorem. Only the vertical component of the vector surface element enters, which explains the absence of $\alpha$ in the answer.

For $\alpha =3$, the divergence vanishes at every point with $r>0$, so its ordinary volume integral, ignoring the undefined single point, is zero. Nevertheless, the outward flux integral through the unit sphere is still $4\pi$:The hypotheses of the divergence theorem fail at the origin. Indeed $|\mathbf{u}|=r^{-2}$ there, so no continuous extension is possible. On the punctured region $\epsilon <r<1$, the divergence theorem remains valid: the outer flux is $4\pi$ and the inner flux is $-4\pi$. Unlike iii, this inner contribution does not tend to zero. The discrepancy therefore comes from the persistent flux at the singular point, rather than a failure of the divergence theorem under its hypotheses.

On the unit sphere, the outward unit normal vector is $\mathbf{n}=\mathbf{r}$ and $\mathbf{u}=\mathbf{r}$. Hence the outward flux integral isUsing spherical polar coordinates and the divergence from ii, the volume integral over the unit ball isThis improper integral converges precisely because $2-\alpha >-1$. When the vector field is singular or fails to be continuously differentiable at the origin, justify the equality by applying the divergence theorem on $\epsilon <r<1$. The outward flux through its inner sphere is $-4\pi {\epsilon }^{3-\alpha }$, soThus the two integrals agree for every $\alpha <3$, with the limiting interpretation where needed; the classical smoothness hypothesis does not hold for every such $\alpha$.

For $r>0$, the chain rule gives $\partial r/\partial x_i=x_i/r$. Applying the product rule to the components of the vector field yieldsThe calculation is pointwise away from the origin; any extension of the vector field to the origin must be considered separately.