5.8.3: Plane Discs
( \newcommand{\kernel}{\mathrm{null}\,}\)
Refer to figure V.2A. The potential at P from the elemental disc is
dψ=−GδM(r2+z2)1/2=−2πGσrδr(r2+z2)1/2.
The potential from the whole disc is therefore
ψ=−2πGσ∫a0rdr(r2+z2)1/2.
The integral is trivial after a brilliant substitution such as X=r2+z2 or r=ztanθ, and we arrive at
ψ=−2πGσ(√z2+a2−z).
This increases to zero as z→∞. We can also write this as
ψ=−2πGmπa2⋅[z(1+a2z2)1/2−z],
and, if you expand this binomially, you see that for large z it becomes, as expected, −Gm/z.