The gravitational
potential V at a distance
x from a
point mass of mass
M can be defined as the work
W that needs to be done by an external agent to bring a unit mass in from infinity to that point: V(\mathbf{x}) = \frac{W}{m} = \frac{1}{m} \int_{\infty}^{x} \mathbf{F} \left(\mathbf{x}'\right) \cdot d\mathbf{x}' = \frac{1}{m} \int_{\infty}^{x} \frac{G m M}{x'^2} dx' = -\frac{G M}{x}, where
G is the
gravitational constant, and
F is the gravitational force. The product
GM is the
standard gravitational parameter and is often known to higher precision than
G or
M separately. The potential has
units of energy per mass, e.g., J/kg in the
MKS system. By convention, it is always negative where it is defined, and as
x tends to infinity, it approaches zero. The
gravitational field, and thus the acceleration of a small body in the space around the massive object, is the negative
gradient of the gravitational potential. Thus the negative of a negative gradient yields positive acceleration toward a massive object. Because the potential has no angular components, its gradient is \mathbf{a} = -\frac{GM}{x^3} \mathbf{x} = -\frac{GM}{x^2} \hat{\mathbf{x}}, where
x is a vector of length
x pointing from the point mass toward the small body and \hat{\mathbf{x}} is a
unit vector pointing from the point mass toward the small body. The magnitude of the acceleration therefore follows an
inverse square law: \|\mathbf{a}\| = \frac{GM}{x^2}. The potential associated with a
mass distribution is the superposition of the potentials of point masses. If the mass distribution is a finite collection of point masses, and if the point masses are located at the points
x1, ...,
xn and have masses
m1, ...,
mn, then the potential of the distribution at the point
x is V(\mathbf{x}) = \sum_{i=1}^n -\frac{Gm_i}{\|\mathbf{x} - \mathbf{x}_i\|}. If the mass distribution is given as a mass
measure dm on three-dimensional
Euclidean space R3, then the potential is the
convolution of with
dm. In good cases this equals the integral V(\mathbf{x}) = -\int_{\R^3} \frac{G}{\|\mathbf{x} - \mathbf{r}\|}\,dm(\mathbf{r}), where is the
distance between the points
x and
r. If there is a function
ρ(
r) representing the density of the distribution at
r, so that , where
dv(
r) is the Euclidean
volume element, then the gravitational potential is the
volume integral V(\mathbf{x}) = -\int_{\R^3} \frac{G}{\|\mathbf{x}-\mathbf{r}\|}\,\rho(\mathbf{r})dv(\mathbf{r}). If
V is a potential function coming from a continuous mass distribution
ρ(
r), then
ρ can be recovered using the
Laplace operator, : \rho(\mathbf{x}) = \frac{1}{4\pi G}\Delta V(\mathbf{x}). This holds pointwise whenever
ρ is continuous and is zero outside of a bounded set. In general, the mass measure
dm can be recovered in the same way if the Laplace operator is taken in the sense of
distributions. As a consequence, the gravitational potential satisfies
Poisson's equation. See also
Green's function for the three-variable Laplace equation and
Newtonian potential. The integral may be expressed in terms of known transcendental functions for all ellipsoidal shapes, including the symmetrical and degenerate ones. These include the sphere, where the three semi axes are equal; the oblate (see
reference ellipsoid) and prolate spheroids, where two semi axes are equal; the degenerate ones where one semi axes is infinite (the elliptical and circular cylinder) and the unbounded sheet where two semi axes are infinite. All these shapes are widely used in the applications of the gravitational potential integral (apart from the constant
G, with 𝜌 being a constant
charge density) to electromagnetism. ==Spherical symmetry==