0=d2xdr02+1r0[4+dlnP0dlnr0]dxdr0+[(3−4γg)dlnP0dlnr0]xr02−1Δ[dlnP0dlnr0⋅ρcρ0(σc26γg)]xr02
where: Δ≡Mr4πr03ρ0, σc2≡3ω22πGρc