Für kugelsymmetrische Massenverteilungen darfst du dir die Masse weiter in einem Punkt konzentriert denken. Gemäß dem Newtonschen Schalentheorem folgt das Gravitationspotential identisch zu dem einer Punktmasse (Eine Verallgemeinerung dieses Satzes gilt sogar in der ART. Gemäß dem Birkhoff-Theorem liefert eine sphärisch symmetrische - ggf. sogar pulsierende - Massenverteilungen im Außenraum exakt die statische Schwarzschild-Geometrie - ohne Gravitationswellen)
Für lediglich näherungsweise kugelsymmetrische Massenverteilungen (z.B. Planeten) verwendet man die sphärische Multipolentwicklung. Diese funktioniert für statische Massen- sowie Ladungsverteilungen exakt identisch.
https://en.m.wikipedia.org/wiki/Multipole_expansion