In mathematics, a fundamental solution for a linear partial differential operator L is a formulation in the language of distribution theory of the older idea of a Green's function (although unlike Green's functions, fundamental solutions do not address boundary conditions). In terms of the Dirac delta function δ(x), a fundamental solution F is a solution of the inhomogeneous equation
Here F is a priori only assumed to be a distribution. This concept has long been utilized for the Laplacian in two and three dimensions. It was investigated for all dimensions for the Laplacian by Marcel Riesz. The existence of a fundamental solution for any operator with constant coefficients — the most important case, directly linked to the possibility of using convolution to solve an arbitrary right hand side — was shown by Bernard Malgrange and Leon Ehrenpreis, and a proof is available in Joel Smoller (1994). In the context of functional analysis, fundamental solutions are usually developed via the Fredholm alternative and explored in Fredholm theory.
Example Consider the following differential equation Lf = sin(x) with
L = d 2 d x 2 . {\displaystyle L={\frac {d^{2}}{dx^{2}}}.}
The fundamental solutions can be obtained by solving LF = δ(x), explicitly,
d 2 d x 2 F ( x ) = δ ( x ) . {\displaystyle {\frac {d^{2}}{dx^{2}}}F(x)=\delta (x)\,.}
Since for the unit step function (also known as the Heaviside function) H we have
d d x H ( x ) = δ ( x ) , {\displaystyle {\frac {d}{dx}}H(x)=\delta (x)\,,}
there is a solution
d d x F ( x ) = H ( x ) + C . {\displaystyle {\frac {d}{dx}}F(x)=H(x)+C\,.}
Here C is an arbitrary constant introduced by the integration. For convenience, set C = −1/2. After integrating d F d x {\displaystyle {\frac {dF}{dx}}} and choosing the new integration constant as zero, one has
F ( x ) = x H ( x ) − 1 2 x = 1 2 | x | . {\displaystyle F(x)=xH(x)-{\frac {1}{2}}x={\frac {1}{2}}|x|~.}
Motivation Once the fundamental solution is found, it is straightforward to find a solution of the original equation, through convolution of the fundamental solution and the desired right hand side. Fundamental solutions also play an important role in the numerical solution of partial differential equations by the boundary element method.
Application to the example Consider the operator L and the differential equation mentioned in the example,
d 2 d x 2 f ( x ) = sin ( x ) . {\displaystyle {\frac {d^{2}}{dx^{2}}}f(x)=\sin(x)\,.}
We can find the solution f ( x ) {\displaystyle f(x)} of the original equation by convolution (denoted by an asterisk) of the right-hand side sin ( x ) {\displaystyle \sin(x)} with the fundamental solution F ( x ) = 1 2 | x | {\textstyle F(x)={\frac {1}{2}}|x|} :
f ( x ) = ( F ∗ sin ) ( x ) := ∫ − ∞ ∞ 1 2 | x − y | sin ( y ) d y . {\displaystyle f(x)=(F*\sin )(x):=\int _{-\infty }^{\infty }{\frac {1}{2}}|x-y|\sin(y)\,dy\,.}
This shows that some care must be taken when working with functions which do not have enough regularity (e.g. compact support, L1 integrability) since, we know that the desired solution is f(x) = −sin(x), while the above integral diverges for all x. The two expressions for f are, however, equal as distributions.
An example that more clearly works
… excerpt ends here. Continue reading the full article.
