If u(x,y,z)=(x2+y2+z2)−1/2, then find the value of ∂x2∂2u+∂y2∂2u+∂z2∂2u.
Technique
Write u=r−1 with r=x2+y2+z2 and use the chain rule through ∂x∂r=rx. Compute uxx once and obtain uyy,uzz by the symmetry of r in x,y,z; add and use x2+y2+z2=r2. State the exclusion of the origin.
Solution
Step 0 — Notation and domain.
Put
r=x2+y2+z2(>0),u=r−1=(x2+y2+z2)−1/2.
u is defined and infinitely differentiable only for (x,y,z)=(0,0,0) — at the origin r=0 and u is not even finite, so the whole computation is carried out on R3∖{(0,0,0)}. This exclusion is part of the answer, not a footnote.
From r2=x2+y2+z2, differentiating partially with respect to x: 2rrx=2x, so
∂x∂r=rx,and similarly∂y∂r=ry,∂z∂r=rz.
Step 1 — First partial derivative.
∂x∂u=drd(r−1)⋅∂x∂r=−r−2⋅rx=−r3x.
(Direct check without the r-substitution: ∂x∂u=−21(x2+y2+z2)−3/2⋅2x=−x(x2+y2+z2)−3/2. ✓)
By the complete symmetry of r in x,y,z (the same computation with y or z in place of x):
∂y2∂2u=−r31+r53y2,∂z2∂2u=−r31+r53z2.
Step 3 — Add the three.
∇2u=∂x2∂2u+∂y2∂2u+∂z2∂2u=−r33+r53(x2+y2+z2).
Since x2+y2+z2=r2,
∇2u=−r33+r53r2=−r33+r33=0.
Hence u satisfies Laplace’s equation everywhere except the origin: u=1/r is a harmonic function on R3∖{0} — the Newtonian (gravitational / Coulomb) potential of a unit point source at the origin.
Step 4 — The general result, as a check and for depth.
For u=rn the identical computation gives ux=nxrn−2 and
∂x2∂2u=nrn−2+n(n−2)x2rn−4,
so, adding the three and using x2+y2+z2=r2,
∇2(rn)=3nrn−2+n(n−2)rn−2=n(n+1)rn−2.
Putting n=−1: ∇2(r−1)=(−1)(0)r−3=0 ✓ — consistent with Step 3. (Note n=0 and n=−1 are the only exponents making rn harmonic in R3.)
Spot check at a point. Take (x,y,z)=(1,2,−2), so r=1+4+4=3, r3=27, r5=243:
Remark (depth, one line). The exclusion is essential: over all of R3, in the sense of distributions ∇2(r1)=−4πδ(r), so the origin carries the whole source. Classically, ∇2u is simply undefined there.