Dans le cas le plus général d'une distribution continue de matière décrite par une densité de masse ρ(r), où r représente le rayon vecteur d'un point quelconque de l'espace, l'énergie potentielle gravitationnelle du système est donnée par la somme de tous les travaux nécessaires pour amener chacune de ses parties depuis l'infini jusqu'à leur position finale. Cette énergie s'écrit alors
Ep=−21∬G∣r−r′∣ρ(r)ρ(r′)drdr′.
Le facteur 1/2 peut se comprendre par le fait que l'on considère l'ensemble des énergie potentielles prises entre deux points de la distribution de masse, chaque paire de points étant comptée deux fois d'où la nécessité de rajouter un facteur 1/2 dans le résultat final.
Autre écriture
En fonction du potentiel gravitationnel
Du fait que la formule générale du potentiel gravitationnel s'écrit
Φ(r)=−G∫∣r−r′∣ρ(r′)dr′,
on peut effectuer une des deux intégrations dans la formule précédente, pour obtenir
Ep=21∫ρ(r)Φ(r)dr.
En fonction du champ gravitationnel
Si l'on connaît le champ gravitationnel g généré par la distribution des sources, on peut réexprimer la formule précédente selon
Ep=−8πG1∫∣g(r)∣2dr.
Cette expression est essentiellement similaires aux énergies électrostatique et magnétostatique, faisant l'une et l'autre intervenir l'intégrale sur tout l'espace du carré de la norme du champ correspondant (électrique et magnétique, respectivement), le tout multiplié par la constante appropriée. Par exemple, la constante est ε/2 pour l'énergie potentielle électrostatique et -1/8πG pour l'énergie potentielle gravitationnelle, car les forces électrostatiques et gravitationnelles font intervenir les constantes 1/4πε et G, et que l'une est répulsive pour des charges de même signe, alors que l'autre est attractive.
Cette expression est assez commode dans le cas du calcul de l'énergie potentielle gravitationnelle d'une distribution de matière à symétrie sphérique.