When solving the equation for the geometry in spherical coordinates, the following equation for the polar angle arises
where is the unknown function, is a nonnegative integer, is the separation constant. Let us expand the derivative in the first term
Note that for the equation takes the form of the Legendre equation
This equation has solutions bounded on the interval only for the eigenvalues , ; these solutions are given by the Rodrigues formula
where are the Legendre polynomials. We will look for solutions of equation (E.1 ) for the same eigenvalues — then it can be written as
Let us make the substitution into equation (E.3 ). The first derivative takes the form
To compute the second derivative, we first compute
Then the second derivative can be written as
After collecting like terms it takes the form
Let us substitute everything into equation (E.3 ), cancelling along the way :
Let us group the coefficients of the derivatives:
thus, after the substitution, equation (E.3 ) takes the form
Let us take equation (E.3 ), set and take into account that in this case its solutions are the Legendre polynomials — we obtain the following equation
which we differentiate times using the Leibniz formula, defined as follows
where is the binomial coefficient.
Let us start with the first term of equation (E.5 ), taking
now let us take , we obtain
for the case everything is trivial
The binomial coefficients are respectively equal to
The result of differentiating the first term of equation (E.5 ) equals
for the second term
for the third term
let us add all three terms
The left-hand side of equation (E.5 ) is identically zero, hence so is its -fold derivative:
The resulting equation coincides with equation (E.4 ), whence , which means the solution of equation (E.1 ) has the form
where are the eigenvalues, are the associated Legendre polynomials. For the derivative vanishes, so nontrivial solutions exist only for .