In

Introduction to Electrodynamics

'. 4th ed., Pearson, 2013. Inner front cover. . $$\backslash nabla^2\; f\; =\; \backslash frac\; +\; \backslash frac\; +\; \backslash frac\; =\; 0.$$ In cylindrical coordinates, $$\backslash nabla^2\; f=\backslash frac\; \backslash frac\; \backslash left(\; r\; \backslash frac\; \backslash right)\; +\; \backslash frac\; \backslash frac\; +\; \backslash frac\; =\; 0.$$ In

_{n}'' and ''G_{n}'' denote normal derivatives on . In view of the conditions satisfied by and , this result simplifies to
$$u(x\text{'},y\text{'},z\text{'})\; =\; \backslash iiint\_V\; G\; f\; \backslash ,\; dV\; +\; \backslash iint\_S\; G\_n\; g\; \backslash ,\; dS.\; \backslash ,$$
Thus the Green's function describes the influence at of the data and . For the case of the interior of a sphere of radius , the Green's function may be obtained by means of a reflection : the source point at distance from the center of the sphere is reflected along its radial line to a point ''P that is at a distance
$$\backslash rho\text{'}\; =\; \backslash frac.\; \backslash ,$$
Note that if is inside the sphere, then ''P′'' will be outside the sphere. The Green's function is then given by
$$\backslash frac\; -\; \backslash frac,\; \backslash ,$$
where denotes the distance to the source point and denotes the distance to the reflected point ''P''′. A consequence of this expression for the Green's function is the Poisson integral formula. Let , , and be

Laplace Equation (particular solutions and boundary value problems)

at EqWorld: The World of Mathematical Equations.

Example initial-boundary value problems

using Laplace's equation from exampleproblems.com. * {{MathWorld , urlname= LaplacesEquation , title= Laplace's Equation

Elliptic partial differential equations Harmonic functions Equations Fourier analysis

mathematics
Mathematics is an area of knowledge that includes the topics of numbers, formulas and related structures, shapes and the spaces in which they are contained, and quantities and their changes. These topics are represented in modern mathematics ...

and physics
Physics is the natural science that studies matter, its fundamental constituents, its motion and behavior through space and time, and the related entities of energy and force. "Physical science is that department of knowledge which rela ...

, Laplace's equation is a second-order partial differential equation
In mathematics, a partial differential equation (PDE) is an equation which imposes relations between the various partial derivatives of a multivariable function.
The function is often thought of as an "unknown" to be solved for, similarly to h ...

named after Pierre-Simon Laplace
Pierre-Simon, marquis de Laplace (; ; 23 March 1749 – 5 March 1827) was a French scholar and polymath whose work was important to the development of engineering, mathematics, statistics, physics, astronomy, and philosophy. He summarized ...

, who first studied its properties. This is often written as
$$\backslash nabla^2\backslash !\; f\; =\; 0$$ or $$\backslash Delta\; f\; =\; 0,$$
where $\backslash Delta\; =\; \backslash nabla\; \backslash cdot\; \backslash nabla\; =\; \backslash nabla^2$ is the Laplace operator,The delta symbol, Δ, is also commonly used to represent a finite change in some quantity, for example, $\backslash Delta\; x\; =\; x\_1\; -\; x\_2$. Its use to represent the Laplacian should not be confused with this use. $\backslash nabla\; \backslash cdot$ is the divergence
In vector calculus, divergence is a vector operator that operates on a vector field, producing a scalar field giving the quantity of the vector field's source at each point. More technically, the divergence represents the volume density of th ...

operator (also symbolized "div"), $\backslash nabla$ is the gradient
In vector calculus, the gradient of a scalar-valued differentiable function of several variables is the vector field (or vector-valued function) \nabla f whose value at a point p is the "direction and rate of fastest increase". If the grad ...

operator (also symbolized "grad"), and $f\; (x,\; y,\; z)$ is a twice-differentiable real-valued function. The Laplace operator therefore maps a scalar function to another scalar function.
If the right-hand side is specified as a given function, $h(x,\; y,\; z)$, we have
$$\backslash Delta\; f\; =\; h.$$
This is called Poisson's equation, a generalization of Laplace's equation. Laplace's equation and Poisson's equation are the simplest examples of elliptic partial differential equations. Laplace's equation is also a special case of the Helmholtz equation
In mathematics, the eigenvalue problem for the Laplace operator is known as the Helmholtz equation. It corresponds to the linear partial differential equation
\nabla^2 f = -k^2 f,
where is the Laplace operator (or "Laplacian"), is the eigenvalu ...

.
The general theory of solutions to Laplace's equation is known as potential theory. The twice continuously differentiable solutions of Laplace's equation are the harmonic functions, which are important in multiple branches of physics, notably electrostatics, gravitation, and fluid dynamics
In physics and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids— liquids and gases. It has several subdisciplines, including ''aerodynamics'' (the study of air and other gases in motion) ...

. In the study of heat conduction, the Laplace equation is the steady-state
In systems theory, a system or a process is in a steady state if the variables (called state variables) which define the behavior of the system or the process are unchanging in time. In continuous time, this means that for those properties ''p' ...

heat equation. In general, Laplace's equation describes situations of equilibrium, or those that do not depend explicitly on time.
Forms in different coordinate systems

In rectangular coordinates,Griffiths, David J.Introduction to Electrodynamics

'. 4th ed., Pearson, 2013. Inner front cover. . $$\backslash nabla^2\; f\; =\; \backslash frac\; +\; \backslash frac\; +\; \backslash frac\; =\; 0.$$ In cylindrical coordinates, $$\backslash nabla^2\; f=\backslash frac\; \backslash frac\; \backslash left(\; r\; \backslash frac\; \backslash right)\; +\; \backslash frac\; \backslash frac\; +\; \backslash frac\; =\; 0.$$ In

spherical coordinates
In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a point is specified by three numbers: the ''radial distance'' of that point from a fixed origin, its ''polar angle'' measu ...

, using the $(r,\; \backslash theta,\; \backslash varphi)$ convention,
$$\backslash nabla^2\; f\; =\; \backslash frac\backslash frac\; \backslash left(r^2\; \backslash frac\backslash right)\; +\; \backslash frac\; \backslash frac\; \backslash left(\backslash sin\backslash theta\; \backslash frac\backslash right)\; +\; \backslash frac\; \backslash frac\; =0.$$
More generally, in arbitrary curvilinear coordinates ,
$$\backslash nabla^2\; f\; =\backslash frac\backslash left(\backslash fracg^\backslash right)\; +\; \backslash frac\; g^\backslash Gamma^n\_\; =0,$$
or
$$\backslash nabla^2\; f\; =\; \backslash frac\; \backslash frac\backslash !\backslash left(\backslash sqrtg^\; \backslash frac\backslash right)\; =0,\; \backslash qquad\; (g=\backslash det\backslash )$$
where is the Euclidean metric tensor
In the mathematical field of differential geometry, a metric tensor (or simply metric) is an additional structure on a manifold (such as a surface) that allows defining distances and angles, just as the inner product on a Euclidean space allows ...

relative to the new coordinates and denotes its Christoffel symbols.
Boundary conditions

The Dirichlet problem for Laplace's equation consists of finding a solution on some domain such that on the boundary of is equal to some given function. Since the Laplace operator appears in the heat equation, one physical interpretation of this problem is as follows: fix the temperature on the boundary of the domain according to the given specification of the boundary condition. Allow heat to flow until a stationary state is reached in which the temperature at each point on the domain doesn't change anymore. The temperature distribution in the interior will then be given by the solution to the corresponding Dirichlet problem. TheNeumann boundary condition
In mathematics, the Neumann (or second-type) boundary condition is a type of boundary condition, named after Carl Neumann.
When imposed on an ordinary or a partial differential equation, the condition specifies the values of the derivative appl ...

s for Laplace's equation specify not the function itself on the boundary of but its normal derivative. Physically, this corresponds to the construction of a potential for a vector field whose effect is known at the boundary of alone. For the example of the heat equation it amounts to prescribing the heat flux through the boundary. In particular, at an adiabatic boundary, the normal derivative of is zero.
Solutions of Laplace's equation are called harmonic functions; they are all analytic within the domain where the equation is satisfied. If any two functions are solutions to Laplace's equation (or any linear homogeneous differential equation), their sum (or any linear combination) is also a solution. This property, called the principle of superposition, is very useful. For example, solutions to complex problems can be constructed by summing simple solutions.
In two dimensions

Laplace's equation in two independent variables in rectangular coordinates has the form $$\backslash frac\; +\; \backslash frac\; \backslash equiv\; \backslash psi\_\; +\; \backslash psi\_\; =\; 0.$$Analytic functions

The real and imaginary parts of a complexanalytic function
In mathematics, an analytic function is a function that is locally given by a convergent power series. There exist both real analytic functions and complex analytic functions. Functions of each type are infinitely differentiable, but complex a ...

both satisfy the Laplace equation. That is, if , and if
$$f(z)\; =\; u(x,y)\; +\; iv(x,y),$$
then the necessary condition that be analytic is that and be differentiable and that the Cauchy–Riemann equations be satisfied:
$$u\_x\; =\; v\_y,\; \backslash quad\; v\_x\; =\; -u\_y.$$
where is the first partial derivative of with respect to .
It follows that
$$u\_\; =\; (-v\_x)\_y\; =\; -(v\_y)\_x\; =\; -(u\_x)\_x.$$
Therefore satisfies the Laplace equation. A similar calculation shows that also satisfies the Laplace equation.
Conversely, given a harmonic function, it is the real part of an analytic function, (at least locally). If a trial form is
$$f(z)\; =\; \backslash varphi(x,y)\; +\; i\; \backslash psi(x,y),$$
then the Cauchy–Riemann equations will be satisfied if we set
$$\backslash psi\_x\; =\; -\backslash varphi\_y,\; \backslash quad\; \backslash psi\_y\; =\; \backslash varphi\_x.$$
This relation does not determine , but only its increments:
$$d\; \backslash psi\; =\; -\backslash varphi\_y\backslash ,\; dx\; +\; \backslash varphi\_x\backslash ,\; dy.$$
The Laplace equation for implies that the integrability condition for is satisfied:
$$\backslash psi\_\; =\; \backslash psi\_,$$
and thus may be defined by a line integral. The integrability condition and Stokes' theorem implies that the value of the line integral connecting two points is independent of the path. The resulting pair of solutions of the Laplace equation are called conjugate harmonic functions. This construction is only valid locally, or provided that the path does not loop around a singularity. For example, if and are polar coordinates and
$$\backslash varphi\; =\; \backslash log\; r,$$
then a corresponding analytic function is
$$f(z)\; =\; \backslash log\; z\; =\; \backslash log\; r\; +\; i\backslash theta.$$
However, the angle is single-valued only in a region that does not enclose the origin.
The close connection between the Laplace equation and analytic functions implies that any solution of the Laplace equation has derivatives of all orders, and can be expanded in a power series, at least inside a circle that does not enclose a singularity. This is in sharp contrast to solutions of the wave equation
The (two-way) wave equation is a second-order linear partial differential equation for the description of waves or standing wave fields — as they occur in classical physics — such as mechanical waves (e.g. water waves, sound waves and sei ...

, which generally have less regularity.
There is an intimate connection between power series and Fourier series
A Fourier series () is a summation of harmonically related sinusoidal functions, also known as components or harmonics. The result of the summation is a periodic function whose functional form is determined by the choices of cycle length (or ''p ...

. If we expand a function in a power series inside a circle of radius , this means that
$$f(z)\; =\; \backslash sum\_^\backslash infty\; c\_n\; z^n,$$
with suitably defined coefficients whose real and imaginary parts are given by
$$c\_n\; =\; a\_n\; +\; i\; b\_n.$$
Therefore
$$f(z)\; =\; \backslash sum\_^\backslash infty\; \backslash left;\; href="/html/ALL/s/a\_n\_r^n\_\backslash cos\_n\_\backslash theta\_-\_b\_n\_r^n\_\backslash sin\_n\_\backslash theta\backslash right.html"\; ;"title="a\_n\; r^n\; \backslash cos\; n\; \backslash theta\; -\; b\_n\; r^n\; \backslash sin\; n\; \backslash theta\backslash right">a\_n\; r^n\; \backslash cos\; n\; \backslash theta\; -\; b\_n\; r^n\; \backslash sin\; n\; \backslash theta\backslash right$$
which is a Fourier series for . These trigonometric functions can themselves be expanded, using multiple angle formulae.
Fluid flow

Let the quantities and be the horizontal and vertical components of the velocity field of a steady incompressible, irrotational flow in two dimensions. The continuity condition for an incompressible flow is that $$u\_x\; +\; v\_y=0,$$ and the condition that the flow be irrotational is that $$\backslash nabla\; \backslash times\; \backslash mathbf\; =\; v\_x\; -\; u\_y\; =\; 0.$$ If we define the differential of a function by $$d\; \backslash psi\; =\; v\; \backslash ,\; dx\; -\; u\; \backslash ,\; dy,$$ then the continuity condition is the integrability condition for this differential: the resulting function is called the stream function because it is constant along flow lines. The first derivatives of are given by $$\backslash psi\_x\; =\; v,\; \backslash quad\; \backslash psi\_y=-u,$$ and the irrotationality condition implies that satisfies the Laplace equation. The harmonic function that is conjugate to is called the velocity potential. The Cauchy–Riemann equations imply that $$\backslash varphi\_x=-u,\; \backslash quad\; \backslash varphi\_y=-v.$$ Thus every analytic function corresponds to a steady incompressible, irrotational, inviscid fluid flow in the plane. The real part is the velocity potential, and the imaginary part is the stream function.Electrostatics

According toMaxwell's equations
Maxwell's equations, or Maxwell–Heaviside equations, are a set of coupled partial differential equations that, together with the Lorentz force law, form the foundation of classical electromagnetism, classical optics, and electric circuits.
...

, an electric field in two space dimensions that is independent of time satisfies
$$\backslash nabla\; \backslash times\; (u,v,0)\; =\; (v\_x\; -u\_y)\backslash hat\; =\; \backslash mathbf,$$
and
$$\backslash nabla\; \backslash cdot\; (u,v)\; =\; \backslash rho,$$
where is the charge density. The first Maxwell equation is the integrability condition for the differential
$$d\; \backslash varphi\; =\; -u\backslash ,\; dx\; -v\backslash ,\; dy,$$
so the electric potential may be constructed to satisfy
$$\backslash varphi\_x\; =\; -u,\; \backslash quad\; \backslash varphi\_y\; =\; -v.$$
The second of Maxwell's equations then implies that
$$\backslash varphi\_\; +\; \backslash varphi\_\; =\; -\backslash rho,$$
which is the Poisson equation. The Laplace equation can be used in three-dimensional problems in electrostatics and fluid flow just as in two dimensions.
In three dimensions

Fundamental solution

A fundamental solution of Laplace's equation satisfies $$\backslash Delta\; u\; =\; u\_\; +\; u\_\; +\; u\_\; =\; -\backslash delta(x-x\text{'},y-y\text{'},z-z\text{'}),$$ where theDirac delta function
In mathematics, the Dirac delta distribution ( distribution), also known as the unit impulse, is a generalized function or distribution over the real numbers, whose value is zero everywhere except at zero, and whose integral over the entir ...

denotes a unit source concentrated at the point . No function has this property: in fact it is a distribution rather than a function; but it can be thought of as a limit of functions whose integrals over space are unity, and whose support (the region where the function is non-zero) shrinks to a point (see weak solution
In mathematics, a weak solution (also called a generalized solution) to an ordinary or partial differential equation is a function for which the derivatives may not all exist but which is nonetheless deemed to satisfy the equation in some precise ...

). It is common to take a different sign convention for this equation than one typically does when defining fundamental solutions. This choice of sign is often convenient to work with because −Δ is a positive operator. The definition of the fundamental solution thus implies that, if the Laplacian of is integrated over any volume that encloses the source point, then
$$\backslash iiint\_V\; \backslash nabla\; \backslash cdot\; \backslash nabla\; u\; \backslash ,\; dV\; =-1.$$
The Laplace equation is unchanged under a rotation of coordinates, and hence we can expect that a fundamental solution may be obtained among solutions that only depend upon the distance from the source point. If we choose the volume to be a ball of radius around the source point, then Gauss' divergence theorem implies that
$$-1=\; \backslash iiint\_V\; \backslash nabla\; \backslash cdot\; \backslash nabla\; u\; \backslash ,\; dV\; =\; \backslash iint\_S\; \backslash frac\; \backslash ,\; dS\; =\; \backslash left.4\backslash pi\; a^2\; \backslash frac\backslash \_.$$
It follows that
$$\backslash frac\; =\; -\backslash frac,$$
on a sphere of radius that is centered on the source point, and hence
$$u\; =\; \backslash frac.$$
Note that, with the opposite sign convention (used in physics
Physics is the natural science that studies matter, its fundamental constituents, its motion and behavior through space and time, and the related entities of energy and force. "Physical science is that department of knowledge which rela ...

), this is the potential
Potential generally refers to a currently unrealized ability. The term is used in a wide variety of fields, from physics to the social sciences to indicate things that are in a state where they are able to change in ways ranging from the simple r ...

generated by a point particle, for an inverse-square law
In science, an inverse-square law is any scientific law stating that a specified physical quantity is inversely proportional to the square of the distance from the source of that physical quantity. The fundamental cause for this can be underst ...

force, arising in the solution of Poisson equation. A similar argument shows that in two dimensions
$$u\; =\; -\backslash frac.$$
where denotes the natural logarithm
The natural logarithm of a number is its logarithm to the base of the mathematical constant , which is an irrational and transcendental number approximately equal to . The natural logarithm of is generally written as , , or sometimes, if ...

. Note that, with the opposite sign convention, this is the potential
Potential generally refers to a currently unrealized ability. The term is used in a wide variety of fields, from physics to the social sciences to indicate things that are in a state where they are able to change in ways ranging from the simple r ...

generated by a pointlike sink
A sink is a bowl-shaped plumbing fixture for washing hands, dishwashing, and other purposes. Sinks have a tap (faucet) that supply hot and cold water and may include a spray feature to be used for faster rinsing. They also include a drain t ...

(see point particle), which is the solution of the Euler equations in two-dimensional incompressible flow.
Green's function

AGreen's function
In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.
This means that if \operatorname is the linear different ...

is a fundamental solution that also satisfies a suitable condition on the boundary of a volume . For instance,
$$G(x,y,z;x\text{'},y\text{'},z\text{'})$$
may satisfy
$$\backslash nabla\; \backslash cdot\; \backslash nabla\; G\; =\; -\backslash delta(x-x\text{'},y-y\text{'},z-z\text{'})\; \backslash qquad\; \backslash text\; V,$$
$$G\; =\; 0\; \backslash quad\; \backslash text\; \backslash quad\; (x,y,z)\; \backslash qquad\; \backslash text\; S.$$
Now if is any solution of the Poisson equation in :
$$\backslash nabla\; \backslash cdot\; \backslash nabla\; u\; =\; -f,$$
and assumes the boundary values on , then we may apply Green's identity, (a consequence of the divergence theorem) which states that
$$\backslash iiint\_V\; \backslash left;\; href="/html/ALL/s/G\_\backslash ,\_\backslash nabla\_\backslash cdot\_\backslash nabla\_u\_-\_u\_\backslash ,\_\backslash nabla\_\backslash cdot\_\backslash nabla\_G\_\backslash right.html"\; ;"title="G\; \backslash ,\; \backslash nabla\; \backslash cdot\; \backslash nabla\; u\; -\; u\; \backslash ,\; \backslash nabla\; \backslash cdot\; \backslash nabla\; G\; \backslash right">G\; \backslash ,\; \backslash nabla\; \backslash cdot\; \backslash nabla\; u\; -\; u\; \backslash ,\; \backslash nabla\; \backslash cdot\; \backslash nabla\; G\; \backslash right$$
The notations ''uspherical coordinates
In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a point is specified by three numbers: the ''radial distance'' of that point from a fixed origin, its ''polar angle'' measu ...

for the source point . Here denotes the angle with the vertical axis, which is contrary to the usual American mathematical notation, but agrees with standard European and physical practice. Then the solution of the Laplace equation with Dirichlet boundary values inside the sphere is given by
$$u(P)\; =\backslash frac\; a^3\backslash left(1-\backslash frac\backslash right)\; \backslash int\_0^\backslash int\_0^\; \backslash frac\; d\backslash theta\text{'}\; \backslash ,\; d\backslash varphi\text{'}$$
where
$$\backslash cos\; \backslash Theta\; =\; \backslash cos\; \backslash theta\; \backslash cos\; \backslash theta\text{'}\; +\; \backslash sin\backslash theta\; \backslash sin\backslash theta\text{'}\backslash cos(\backslash varphi\; -\backslash varphi\text{'})$$
is the cosine of the angle between and . A simple consequence of this formula is that if is a harmonic function, then the value of at the center of the sphere is the mean value of its values on the sphere. This mean value property immediately implies that a non-constant harmonic function cannot assume its maximum value at an interior point.
Laplace's spherical harmonics

Laplace's equation inspherical coordinates
In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a point is specified by three numbers: the ''radial distance'' of that point from a fixed origin, its ''polar angle'' measu ...

is:
$$\backslash nabla^2\; f\; =\; \backslash frac\; \backslash frac\backslash left(r^2\; \backslash frac\backslash right)\; +\; \backslash frac\; \backslash frac\backslash left(\backslash sin\backslash theta\; \backslash frac\backslash right)\; +\; \backslash frac\; \backslash frac\; =\; 0.$$
Consider the problem of finding solutions of the form . By separation of variables, two differential equations result by imposing Laplace's equation:
$$\backslash frac\backslash frac\backslash left(r^2\backslash frac\backslash right)\; =\; \backslash lambda,\backslash qquad\; \backslash frac\backslash frac\backslash frac\backslash left(\backslash sin\backslash theta\; \backslash frac\backslash right)\; +\; \backslash frac\backslash frac\backslash frac\; =\; -\backslash lambda.$$
The second equation can be simplified under the assumption that has the form . Applying separation of variables again to the second equation gives way to the pair of differential equations
$$\backslash frac\; \backslash frac\; =\; -m^2$$
$$\backslash lambda\backslash sin^2\backslash theta\; +\; \backslash frac\; \backslash frac\; \backslash left(\backslash sin\backslash theta\; \backslash frac\backslash right)\; =\; m^2$$
for some number . A priori, is a complex constant, but because must be a periodic function
A periodic function is a function that repeats its values at regular intervals. For example, the trigonometric functions, which repeat at intervals of 2\pi radians, are periodic functions. Periodic functions are used throughout science to des ...

whose period evenly divides , is necessarily an integer and is a linear combination of the complex exponentials . The solution function is regular at the poles of the sphere, where . Imposing this regularity in the solution of the second equation at the boundary points of the domain is a Sturm–Liouville problem that forces the parameter to be of the form for some non-negative integer with ; this is also explained below in terms of the orbital angular momentum. Furthermore, a change of variables transforms this equation into the Legendre equation, whose solution is a multiple of the associated Legendre polynomial . Finally, the equation for has solutions of the form ; requiring the solution to be regular throughout forces .Physical applications often take the solution that vanishes at infinity, making . This does not affect the angular portion of the spherical harmonics.
Here the solution was assumed to have the special form . For a given value of , there are independent solutions of this form, one for each integer with . These angular solutions are a product of trigonometric function
In mathematics, the trigonometric functions (also called circular functions, angle functions or goniometric functions) are real functions which relate an angle of a right-angled triangle to ratios of two side lengths. They are widely used in al ...

s, here represented as a complex exponential, and associated Legendre polynomials:
$$Y\_\backslash ell^m\; (\backslash theta,\; \backslash varphi\; )\; =\; N\; e^\; P\_\backslash ell^m\; (\backslash cos\; )$$
which fulfill
$$r^2\backslash nabla^2\; Y\_\backslash ell^m\; (\backslash theta,\; \backslash varphi\; )\; =\; -\backslash ell\; (\backslash ell\; +\; 1\; )\; Y\_\backslash ell^m\; (\backslash theta,\; \backslash varphi\; ).$$
Here is called a spherical harmonic function of degree and order , is an associated Legendre polynomial, is a normalization constant, and and represent colatitude and longitude, respectively. In particular, the colatitude
In a spherical coordinate system, a colatitude is the complementary angle of a given latitude, i.e. the difference between a right angle and the latitude. Here Southern latitudes are defined to be negative, and as a result the colatitude is a non ...

, or polar angle, ranges from at the North Pole, to at the Equator, to at the South Pole, and the longitude
Longitude (, ) is a geographic coordinate that specifies the east– west position of a point on the surface of the Earth, or another celestial body. It is an angular measurement, usually expressed in degrees and denoted by the Greek lett ...

, or azimuth
An azimuth (; from ar, اَلسُّمُوت, as-sumūt, the directions) is an angular measurement in a spherical coordinate system. More specifically, it is the horizontal angle from a cardinal direction, most commonly north.
Mathematicall ...

, may assume all values with . For a fixed integer , every solution of the eigenvalue problem
$$r^2\backslash nabla^2\; Y\; =\; -\backslash ell\; (\backslash ell\; +\; 1\; )\; Y$$
is a linear combination of . In fact, for any such solution, is the expression in spherical coordinates of a homogeneous polynomial
In mathematics, a homogeneous polynomial, sometimes called quantic in older texts, is a polynomial whose nonzero terms all have the same degree. For example, x^5 + 2 x^3 y^2 + 9 x y^4 is a homogeneous polynomial of degree 5, in two variables; ...

that is harmonic (see below), and so counting dimensions shows that there are linearly independent such polynomials.
The general solution to Laplace's equation in a ball centered at the origin is a linear combination of the spherical harmonic functions multiplied by the appropriate scale factor ,
$$f(r,\; \backslash theta,\; \backslash varphi)\; =\; \backslash sum\_^\backslash infty\; \backslash sum\_^\backslash ell\; f\_\backslash ell^m\; r^\backslash ell\; Y\_\backslash ell^m\; (\backslash theta,\; \backslash varphi\; ),$$
where the are constants and the factors are known as solid harmonics. Such an expansion is valid in the ball
A ball is a round object (usually spherical, but can sometimes be ovoid) with several uses. It is used in ball games, where the play of the game follows the state of the ball as it is hit, kicked or thrown by players. Balls can also be used ...

$$r\; <\; R\; =\; \backslash frac.$$
For $r\; >\; R$, the solid harmonics with negative powers of $r$ are chosen instead. In that case, one needs to expand the solution of known regions in Laurent series (about $r=\backslash infty$), instead of Taylor series
In mathematics, the Taylor series or Taylor expansion of a function is an infinite sum of terms that are expressed in terms of the function's derivatives at a single point. For most common functions, the function and the sum of its Taylor seri ...

(about $r\; =\; 0$), to match the terms and find $f^m\_\backslash ell$.
Electrostatics

Let $\backslash mathbf$ be the electric field, $\backslash rho$ be the electric charge density, and $\backslash varepsilon\_0$ be the permittivity of free space. Then Gauss's law for electricity (Maxwell's first equation) in differential form statesGriffiths, David J. ''Introduction to Electrodynamics''. Fourth ed., Pearson, 2013. Chapter 2: Electrostatics. p. 83-4. . $$\backslash nabla\; \backslash cdot\; \backslash mathbf\; =\; \backslash frac.$$ Now, the electric field can be expressed as the negative gradient of the electric potential $V$, $$\backslash mathbf\; E=-\backslash nabla\; V,$$ if the field is irrotational, $\backslash nabla\; \backslash times\; \backslash mathbf\; =\; \backslash mathbf$. The irrotationality of $\backslash mathbf$ is also known as the electrostatic condition. $$\backslash nabla\backslash cdot\backslash mathbf\; E\; =\; \backslash nabla\backslash cdot(-\backslash nabla\; V)=-\backslash nabla^2\; V$$ $$\backslash nabla^2\; V\; =\; -\backslash nabla\backslash cdot\backslash mathbf\; E$$ Plugging this relation into Gauss's law, we obtain Poisson's equation for electricity, $$\backslash nabla^2\; V\; =\; -\backslash frac.$$ In the particular case of a source-free region, $\backslash rho\; =\; 0$ and Poisson's equation reduces to Laplace's equation for the electric potential. If the electrostatic potential $V$ is specified on the boundary of a region $\backslash mathcal$, then it is uniquely determined. If $\backslash mathcal$ is surrounded by a conducting material with a specified charge density $\backslash rho$, and if the total charge $Q$ is known, then $V$ is also unique.Griffiths, David J. ''Introduction to Electrodynamics''. Fourth ed., Pearson, 2013. Chapter 3: Potentials. p. 119-121. . A potential that doesn't satisfy Laplace's equation together with the boundary condition is an invalid electrostatic potential.Gravitation

Let $\backslash mathbf$ be the gravitational field, $\backslash rho$ the mass density, and $G$ the gravitational constant. Then Gauss's law for gravitation in differential form is $$\backslash nabla\backslash cdot\backslash mathbf\; g=-4\backslash pi\; G\backslash rho.$$ The gravitational field is conservative and can therefore be expressed as the negative gradient of the gravitational potential: $$\backslash begin\; \backslash mathbf\; g\; \&=\; -\backslash nabla\; V,\; \backslash \backslash \; \backslash nabla\backslash cdot\backslash mathbf\; g\; \&=\; \backslash nabla\backslash cdot(-\backslash nabla\; V)\; =\; -\backslash nabla^2\; V,\; \backslash \backslash \; \backslash implies\backslash nabla^2\; V\; \&=\; -\backslash nabla\backslash cdot\backslash mathbf\; g.\; \backslash end$$ Using the differential form of Gauss's law of gravitation, we have $$\backslash nabla^2\; V\; =\; 4\backslash pi\; G\backslash rho,$$ which is Poisson's equation for gravitational fields. In empty space, $\backslash rho=0$ and we have $$\backslash nabla^2\; V\; =\; 0,$$ which is Laplace's equation for gravitational fields.In the Schwarzschild metric

S. Persides solved the Laplace equation in Schwarzschild spacetime on hypersurfaces of constant . Using the canonical variables , , the solution is $$\backslash Psi(r,\backslash theta,\backslash varphi)\; =\; R(r)Y\_l(\backslash theta,\backslash varphi),$$ where is a spherical harmonic function, and $$R(r)\; =\; (-1)^l\backslash fracP\_l\backslash left(1-\backslash frac\backslash right)+(-1)^\backslash fracQ\_l\backslash left(1-\backslash frac\backslash right).$$ Here and are Legendre functions of the first and second kind, respectively, while is the Schwarzschild radius. The parameter is an arbitrary non-negative integer.See also

* 6-sphere coordinates, a coordinate system under which Laplace's equation becomes ''R''-separable *Helmholtz equation
In mathematics, the eigenvalue problem for the Laplace operator is known as the Helmholtz equation. It corresponds to the linear partial differential equation
\nabla^2 f = -k^2 f,
where is the Laplace operator (or "Laplacian"), is the eigenvalu ...

, a general case of Laplace's equation.
* Spherical harmonic
In mathematics and physical science, spherical harmonics are special functions defined on the surface of a sphere. They are often employed in solving partial differential equations in many scientific fields.
Since the spherical harmonics form ...

* Quadrature domains
* Potential theory
* Potential flow
In fluid dynamics, potential flow (or ideal flow) describes the velocity field as the gradient of a scalar function: the velocity potential. As a result, a potential flow is characterized by an irrotational velocity field, which is a valid appr ...

* Bateman transform
* Earnshaw's theorem uses the Laplace equation to show that stable static ferromagnetic suspension is impossible
* Vector Laplacian
In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols \nabla\cdot\nabla, \nabla^2 (where \nabla is the ...

* Fundamental solution
Notes

References

Further reading

* * * * * *External links

*Laplace Equation (particular solutions and boundary value problems)

at EqWorld: The World of Mathematical Equations.

Example initial-boundary value problems

using Laplace's equation from exampleproblems.com. * {{MathWorld , urlname= LaplacesEquation , title= Laplace's Equation

Elliptic partial differential equations Harmonic functions Equations Fourier analysis

Equation
In mathematics, an equation is a formula that expresses the equality of two expressions, by connecting them with the equals sign . The word ''equation'' and its cognates in other languages may have subtly different meanings; for example, in ...