This post deals with the familiar (to the physics student) Laplace’s equation. I am solving this equation in the context of physics, instead of a pure mathematical perspective. This problem is considered most extensively in the context of electrostatics. This equation is usually considered in the spherical polar coordinate system. A lot of finer details are also considered in a mathematical physics course where the topic of spherical harmonics is discussed. Assuming spherical-polar coordinates, Laplace’s equation is
Suppose that the function . Furthermore, by the method of separation of variables we suppose that the solution is a product of eigenfunctions of the form . Hence, Laplace’s equation becomes
Furthermore, we can separate further the term into . Rewriting (2) and multiplying by , we get
Bringing the radial and angular component to the other side of the equation and setting the azimuthal component equal to a separation constant , yielding
Solving the right-hand side of the equation we get
Now, we set the azimuthal component equal to and carry out the derivative in the angular component to get the Associated Legendre equation:
where I have let . The solutions to this equation are known as the associated Legendre functions. However, instead of solving this difficult equation by brute force methods (i.e. power series method), we consider the case for which . In this case, Eq.(5) simplifies to Legendre’s differential equation discussed previously. Instead of quoting the explicit form of the Legendre polynomials, we equivalently state the Rodrigues formula for these polynomials
However, the solutions to Eq.(5) are the associated Legendre functions, not polynomials. Therefore, we use the following to determine the associated Legendre functions from the Legendre polynomials:
For real solutions and for complex solutions , thus we can write the solutions as series whose indices run from 0 to l for the real-values and from -l to l for the complex-valued solutions.
Now we turn to the radial part which upon substitution of in place of the angular component, we get
and if we evaluate the derivatives in the first term we get an Euler equation of the form
Let us assume that the solution can be represented as
Taking the necessary derivatives and substituting into Eq.(8) gives
Simplifying the powers of r we find that
Now, since , this means that
We can simplify this further to get
Factoring out , we arrive at
Since this must be true for all values of l, the solution therefore becomes
Thus, we can write the solution to Laplace’s equation as
for real solutions, and
for complex solutions.
3 thoughts on “Solution to Laplace’s Equation”
How did you factor j + l out of j^2 + l^2 + j + l = 0?
(j + l)(j + l + 1) = j^2 + l^2 + j + l + 2jl
LikeLiked by 1 person
Hello, thanks for pointing that out. It looks like I made a few sign mistakes that carried through to the values of j. I have updated my solution with my correction.
for some reason my l’s look like larger 1’s…
LikeLiked by 1 person