Assume Newtonian spherical
hydrostatic equilibrium, neglect
rotation and thermal
density changes, and require continuous
pressure at the layer boundary. Write
C=M1−34πρ2R13,D=34πρ2,
so the mantle
mass profile is
M(r)=C+Dr3. Integrating
dP/dr=−Gρ2M(r)/r2 inward from
P(Rp)=P0 gives the boundary
pressurePb=P0+Gρ2[C(R11−Rp1)+2D(Rp2−R12)].
New to topics? Read the docs here!