# Computing gravitational potential for a point inside the distribution

 P: 1,308 The gravitational potential, U, can be calculated at any point, $$${\rm{\vec r}}$$$, for a mass density distribution, $$${\rm{\rho (r)}}$$$, using the formula: $$${\rm{U = - }}\int_{{\rm{all space}}} {G\frac{{{\rm{\rho }}({\rm{\vec r')}}}}{{\left| {{\rm{ \vec r - \vec r' }}} \right|}}} \,\,\,d^3 {\rm{r'}}$$$. See: http://scienceworld.wolfram.com/phys...Potential.html My question is how is this calculated for points inside the distribution. For points outside the distribution, $$${\rm{\rho (r)}}$$$ is zero, and there is no problem. But inside the distribution where $$${\rm{\rho (r)}}$$$ is not zero, there will be points where $$${{\rm{\vec r - \vec r' }}}$$$ does become zero and send the integrand to infinity. Wouldn't this cause a problem? And how would you work around it? Thanks. Could it be as simple as adding a small constant in the denominator that goes to zero after the integration, such as: $$${\rm{U(\vec r) = }}\mathop {\lim }\limits_{\varepsilon \to 0} {\rm{( - }}\int_{{\rm{all space}}} {G\frac{{{\rm{\rho }}({\rm{\vec r')}}}}{{\left| {{\rm{ \vec r - \vec r' }}} \right| + \varepsilon }}} \,\,\,d^3 {\rm{r')}}$$$
 Emeritus Sci Advisor PF Gold P: 29,238 Er... no. One uses techniques such as residues to do such integration that contains poles. Zz.
P: 1,308
 Quote by ZapperZ Er... no. One uses techniques such as residues to do such integration that contains poles. Zz.
As I recall, residuals have to do with the math of complex numbers, which is a 2D construction resolved with line integrals, when here we have a 3D pole/infinity.

Are you referring to a Stoke's Theorm or a Divergence Theorm to turn integration throughout a volume which has an infinity into a surface integral that surrounds the infinity and so does not integrate through the infinity?

Emeritus
PF Gold
P: 29,238
Computing gravitational potential for a point inside the distribution

 Quote by Mike2 As I recall, residuals have to do with the math of complex numbers, which is a 2D construction resolved with line integrals, when here we have a 3D pole/infinity. Are you referring to a Stoke's Theorm or a Divergence Theorm to turn integration throughout a volume which has an infinity into a surface integral that surrounds the infinity and so does not integrate through the infinity?
Nah-ah. You only make use of the complex plane when you do the contour integrals. Look at your complex analysis text. You'll see many such integrals (real ones) being done using this technique. This integral is very common in E&M.

Zz.
P: 1,308
 Quote by ZapperZ Nah-ah. You only make use of the complex plane when you do the contour integrals. Look at your complex analysis text. You'll see many such integrals (real ones) being done using this technique. This integral is very common in E&M. Zz.
I've seen where they use contour integrals to integrate real integrands. They substitute a complex variable for the real variable and do a contour integral. But I've not seen them do this for integrands of more than one real variable, have you? If you don't feel like giving a detailed explanation, perhaps you'd give a few keywords I can look up. Thanks.

 Related Discussions Classical Physics 0 Earth 1 Classical Physics 32 Classical Physics 3