Next Article in Journal
Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant
Next Article in Special Issue
Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction
Previous Article in Journal
Efficient Time and Space Representation of Uncertain Event Data
Previous Article in Special Issue
Simple Iterative Method for Generating Targeted Universal Adversarial Perturbations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application

by
Kyle R. Bryenton
1,†,
Andrew R. Cameron
2,†,
Keegan L. A. Kirk
1,†,
Nasser Saad
1,*,†,
Patrick Strongman
1,† and
Nikita Volodin
1,†
1
School of Mathematical and Computational Sciences, University of Prince Edward Island, Charlottetown, PE C1A 4P3, Canada
2
Department of Physics, University of Prince Edward Island, Charlottetown, PE C1A 4P3, Canada
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Algorithms 2020, 13(11), 286; https://doi.org/10.3390/a13110286
Submission received: 24 October 2020 / Revised: 5 November 2020 / Accepted: 6 November 2020 / Published: 9 November 2020

Abstract

:
The analysis of many physical phenomena is reduced to the study of linear differential equations with polynomial coefficients. The present work establishes the necessary and sufficient conditions for the existence of polynomial solutions to linear differential equations with polynomial coefficients of degree n, n 1 , and n 2 respectively. We show that for n 3 the necessary condition is not enough to ensure the existence of the polynomial solutions. Applying Scheffé’s criteria to this differential equation we have extracted n generic equations that are analytically solvable by two-term recurrence formulas. We give the closed-form solutions of these generic equations in terms of the generalized hypergeometric functions. For arbitrary n, three elementary theorems and one algorithm were developed to construct the polynomial solutions explicitly along with the necessary and sufficient conditions. We demonstrate the validity of the algorithm by constructing the polynomial solutions for the case of n = 4 . We also demonstrate the simplicity and applicability of our constructive approach through applications to several important equations in theoretical physics such as Heun and Dirac equations.

1. Introduction

Differential equations with polynomial coefficients have played an important role not only in understanding engineering and physics problems [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], but also as a source of inspiration for some of the most crucial results in special functions and orthogonal polynomials [33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]. The classical differential equation
α 2 r 2 + α 1 r + α 0 y ( r ) + β 1 r + β 0 y ( r ) τ 0 y ( r ) = 0 ,
has been a cornerstone of many fundamental results in special functions since the early work of L. Euler [33], see also References [47,48,49]. More recently, the Heun-type equation [41,42,43],
r α 3 r 2 + α 2 r + α 1 y ( r ) + β 2 r 2 + β 1 r + β 0 y ( r ) τ 1 r + τ 0 y ( r ) = 0 ,
has become a classic equation in mathematical physics. These two equations are members of the general differential equation
P n ( r ) y ( r ) + Q n 1 ( r ) y ( r ) R n 2 ( r ) y ( r ) = 0 , n 2 , P n ( r ) = k = 0 n α k r k , Q n 1 ( r ) = k = 0 n 1 β k r k , R n 2 ( r ) = k = 0 n 2 τ k r k ,
where { α k } k = 0 n , { β k } k = 0 n 1 , { τ k } k = 0 n 2 are real parameters. We assume that the leading term of R n 2 0 and at least one of the leading terms α n and β n 1 is not zero. The purpose of the present work is to provide an answer to the following question:
“Under what conditions on the equation parameters { α k } k = 0 n , { β k } k = 0 n 1 , and { τ k } k = 0 n 2 does the differential Equation (3) have polynomial solutions? If it does have polynomial solutions, how can we construct them?”
With the general theory of linear differential equations as a guide [50,51], the zeros of the leading polynomial P n ( r ) classify the solutions of the Equation (3). The point r = ζ is called an ordinary point if Q n 1 ( r ) / P n ( r ) and R n 2 ( r ) / P n ( r ) are analytic functions at r = ζ , or a regular singular point if Q n 1 ( r ) / P n ( r ) and R n 2 ( r ) / P n ( r ) are not analytic at r = ζ but the products ( r ζ ) Q n 1 ( r ) / P n ( r ) and ( r ζ ) 2 R n 2 ( r ) / P n ( r ) are analytic at that point. Using this characterization of singularities, S. Bochner [48] classified the polynomial solutions of (1) in terms of the classical orthogonal polynomials. A different approach which depends on the parameters of the leading coefficient P n ( r ) , was introduced in Reference [13] to study the polynomial solutions of (1). It was shown that there are seven possible nonzero leading polynomials depending on the combination of the nonzero α i parameters, for i = 0 , 1 , 2 . By analyzing each of these cases, the authors [13] were able to explicitly construct all the polynomial solutions of Equation (1) in terms of hypergeometric functions along with the associated weight functions.
The nth-degree leading polynomial coefficient P n ( r ) of the differential Equation (3) contains ( 2 n + 1 1 ) nonzero polynomials that depend on the nonzero values of α i , i = 0 , 1 , , n . Out of this polynomial set, there are 2 n differential equations for which r = 0 is an ordinary point, ( 2 n 1 + 2 n 2 ) differential equations for which r = 0 is a regular singular point, and the remaining ( 2 n 2 1 ) differential equations with irregular singular points that fall outside of the scope of this present work.
The criteria for polynomial solutions of second-order linear differential Equation (3) was introduced, using the Asymptotic Iteration Method, in References [52,53]:
Theorem 1.
[52,53] Given λ 0 λ 0 ( x ) and s 0 s 0 ( x ) in C , the differential equation y = λ 0 ( x ) y + s 0 ( x ) y has the general solution
y = exp x α ( t ) d t C 2 + C 1 x exp t λ 0 ( τ ) + 2 α ( τ ) d τ d t ,
if, for some positive integer m > 0 ,
s m λ m = s m 1 λ m 1 = α ( x ) , or δ m ( x ) = λ m s m 1 λ m 1 s m = 0 ,
where λ m = λ m 1 + s m 1 + λ 0 λ m , and s m = s m 1 + s 0 λ m . Further, the differential equation has a polynomial solution of degree m, if
δ m ( x ) = λ m s m 1 λ m 1 s m = 0 .
With λ 0 ( r ) = Q n 1 ( r ) / P n ( r ) and s 0 ( r ) = R n 2 ( r ) / P n ( r ) , a simple algorithm based on Theorem 1 can be used to examine the polynomial solutions of (3).
The present work establishes polynomial solutions in simple and constructive forms to (3) along with all necessary and sufficient conditions they require. This work is done independently from that of Theorem 1, although can be verified using the Asymptotic Iteration Method.
The paper is organized as follows:
  • In Section 2, we introduce the concept of a necessary but not sufficient condition to find polynomial solutions, and then demonstrate how the inverse square-root potential [22,23,24,25] is one such example that does not have polynomial solutions yet satisfies the necessary condition.
  • In Section 3, Scheffé criteria is revised to analyze Equation (3) and to generate exactly solvable equations where the coefficients of their power series solutions are easily computed using a two-term recurrence relation. The closed form solutions of these equations are written in terms of the generalized hypergeometric functions
    p F q α 1 , , α p β 1 , β q ; r = k = 0 j = 1 p ( α p ) k j = 1 q ( α q ) k r k k ! ,
    where the Pochhammer symbol ( α ) k is defined in terms of the Gamma function as ( α ) k = α ( α + 1 ) ( α + k 1 ) = Γ ( α + k ) / Γ ( α ) .
  • In Section 4 we present three theorems corresponding to the cases: (1) α 0 0 , (2) α 0 = 0 , α 1 0 , and (3) α 0 = α 1 = β 0 = 0 , α 2 0 which are required to establish the polynomial solutions of Equation (3). We shall show that for the mth-degree polynomial solutions, there are n + m 1 conditions that ultimately assemble these polynomials. At the end of the section we briefly discuss the Mathematica® program used to generate the solutions of these equations.
  • Lastly, Section 5 demonstrates the validity of these constructions through the application of our results to some problems that have appeared in mathematics and physics literature.

2. A Necessary but Not Sufficient Condition

We begin by stating the necessary condition for the existence of polynomial solutions to the general differential equation given by (3):
Theorem 2.
A necessary condition for the differential Equation (3) to have a polynomial solution of degree m is
τ n 2 = α n m ( m 1 ) + β n 1 m , n = 2 , 3 , .
Proof. 
Substitute the polynomial solution y = i = 0 m C i r i into the differential Equation (3). After collecting the common terms of the common r power, the necessary condition for the leading coefficient C m to be nonzero is τ n 2 = α n m ( m 1 ) + β n 1 m for all m 0 follows. □
For n = 2 , Equation (7) represents both the necessary and sufficient conditions for the polynomial solutions of (3). Although the condition (7) is necessary, it is not sufficient for n 3 to guarantee the polynomial solutions of (3). There are precisely n 2 additional conditions that are sufficient to ensure the existence of such polynomials. These additional conditions can be understood as constraints that relate the remaining parameters of R n 2 ( r ) with the coefficients { α k } k = 0 n and { β k } k = 0 n 1 in P n ( r ) and Q n 1 ( r ) , respectively. In later sections, we shall devise a procedure to find these sufficient conditions and provide a method of evaluating the corresponding polynomial solutions explicitly.
The remainder of the section investigates the inverse square-root potential. We shall show that despite satisfying the necessary condition in Theorem 2 that the Schrödinger equation with this potential as mentioned in the literature has no polynomial solutions, thus demonstrating that the condition (7) is necessary but not sufficient.

The Inverse Square-Root Potential

The exact solutions of the radial Schrödinger equation
d 2 ψ ( r ) d r 2 + v r + l ( l + 1 ) r 2 ψ ( r ) = E ψ ( r ) , v > 0 , r > 0 , ψ ( 0 ) = ψ ( ) = 0 ,
were recently studied by a number of authors [21,22,23,24,25]. The differential Equation (8) has two singular points, one regular at r = 0 and another irregular at r = of rank 3 [50]. To analyze these singular points further, we employ a change of variables z = 2 κ r , for some constant κ to be determined shortly. Equation (8) may then be expressed as
z 2 ψ ( z ) z ψ ( z ) + 2 v κ 3 / 2 + E κ 2 z z 3 4 l ( l + 1 ) ψ ( z ) = 0 , ψ ( 0 ) = ψ ( ) = 0 .
As z 0 , the solution to the differential Equation (9) behaves asymptotically as the solution of the Euler equation, z 2 ψ ( z ) z ψ ( z ) 4 l ( l + 1 ) ψ ( z ) = 0 , which has a physically acceptable solution ψ ( z ) z 2 ( l + 1 ) . Meanwhile as z , the solution to the differential equation behaves asymptotically as the solution of
ψ ( z ) + 2 v κ 3 / 2 z z ψ ( z ) = 0 .
To examine the solution of this equation we complete the square to
ψ ( z ) + v 2 κ 3 z 2 ψ ( z ) = v 2 κ 3 2 ψ ( z ) .
The substitution x = ( ( v / 2 κ 3 ) z ) allows us to compare the solution of the equation to the solvable Schröinger equation with the harmonic oscillator potential, as z , via
ψ ( z ) = exp z ( u + λ ) d u = exp z 2 2 + λ z , λ = v 2 κ 3 .
Hence, we may assume the solution of Equation (9) to be
ψ ( z ) = N z 2 ( l + 1 ) exp z 2 2 + λ z f ( z ) , λ = v 2 κ 3 < 0 ,
up to a normalization constant N. Upon substituting the ansatz solution (12) into Equation (9), we obtain a second-order differential equation (an example of a biconfluent Heun equation [41]) for f ( z ) as
z f ( z ) 2 z 2 + 2 λ z 3 4 l f ( z ) + λ 2 4 ( l + 1 ) z ( 3 + 4 l ) λ f ( z ) = 0 ,
subject to κ = E . Clearly, Equation (13) is of the form of Equation (3) with n = 3 with the parameters α 3 = α 2 = 0 , α 1 = 1 , β 2 = 2 , β 1 = 2 λ , β 0 = 3 + 4 l , τ 1 = λ 2 + 4 ( l + 1 ) , and τ 0 = ( 3 + 4 l ) λ . Using Theorem 2, for the existence of an mth-degree polynomial solution it is necessary that
λ 2 4 ( l + 1 ) = 2 m , m = 0 , 1 , 2 , .
Note that this result can be also deduced using the Bethe Ansatz method [12,54]. On the other hand, an application of the Frobenius method establishes the following three-term recurrence relation for the coefficients of the polynomial solutions f m ( r ) = j = 0 m C j r j :
( j + 1 ) j + 3 + 4 l C j + 1 λ 2 j + 3 + 4 l C j 2 j λ 2 + 4 l + 2 C j 1 = 0 ,
where j = 0 , 1 , 2 , , m , m + 1 and C j = 0 for j < 0 . Implementing the necessary condition (14), with λ = 4 + 4 l + 2 m < 0 , Equation (15) then reads
( j + 1 ) j + 3 + 4 l C j + 1 + 4 + 4 l + 2 m 2 j + 3 + 4 l C j + 2 m j + 1 C j 1 = 0 .
The consistency of the linear homogeneous system generated by the three-term recurrence relation (16), for j = 0 , 1 , , m + 1 , enforces the vanishing of the determinant, denoted by Δ ( m + 1 ) , that establish the sufficient condition
Δ m + 1 = S 0 T 1 γ 1 S 1 T 2 γ 2 S 2 T 3 γ m 2 S m 2 T m 1 γ m 1 S m 1 T m γ m S m , where S j = ( 2 j + 3 + 4 l ) 4 + 4 l + 2 m , T j = j ( j + 2 + 4 l ) , γ j = 2 ( m j + 1 ) .
Therefore, if the determinant Δ m + 1 is different than zero, the only possible solution to the linear system is the trivial solution (i.e., the zero solution). Hence any nonzero value of the determinant Δ m + 1 is an indication that a polynomial solution is not possible, which is the case for the coefficients of Equation (13).
Remark 1.
The point here is not solving the Schrödinger equation with the inverse-root potential; this is already done in several manuscripts [21,22,23,24,25]. We stress that the necessary condition is not enough to guarantee polynomial solutions. Also, the necessary and sufficient constraints must have common roots that are physically acceptable.
Remark 2.
In Reference [55], the author claims the existence of polynomial solutions of a differential equation which matches (13) that satisfy the necessary condition (7). The author missed the possibility of the non-vanishing determinant (4), in his work, for some β 0 and β 2 .
Remark 3.
For the applications in physics, there is nothing special about the necessity of the polynomial solutions. However, the problem of the existence of polynomial solutions is important. Studying the problem in its general form and developing mathematical tools to treat it is significant on its own.

3. Scheffé’s Criteria: Two-Term Recurrence Relation

Generally speaking, recurrence relations with more than two terms are difficult to solve explicitly. Differential equations that are known to have a two-term recurrence relation for their power series solutions guarantee the solvability of such equations.
In a paper presented to the American Mathematical Society (1941), H. Scheffé [34] establishes criteria for the necessary and sufficient conditions for differential equations of the form
p 2 ( r ) u ( r ) + p 1 ( r ) u ( r ) + p 0 ( r ) u ( r ) = 0 ,
to have a two-term recurrence relation. Here p j ( r ) 0 , j = 0 , 1 , 2 , are analytic functions (not necessarily polynomials) in some region consisting of all points in an arbitrary neighbourhood of a regular singular point r = r 0 except the point r 0 itself. We adopt this criterion to examine the differential Equation (18) with
p 2 ( r ) = k = 0 n α k r k , p 1 ( r ) = k = 0 n 1 β k r k , p 0 ( r ) = k = 0 n 2 γ k r k , n 2 ,
and provide a formula for the two-term recurrence relation. Without loss of generality, we shall take the ordinary or singular point as r = 0 , otherwise a simple shifting of r is applied first.
Theorem 3.
The necessary and sufficient conditions for the differential Equation (18) to have a two-term recurrence relation between successive coefficients in its series solution, relative to the ordinary or regular singular point r = 0 , is that in the neighbourhood of r = 0 Equation (18) can be written as:
[ q 2 , 0 + q 2 , h r h ] r 2 m u ( r ) + [ q 1 , 0 + q 1 , h r h ] r 1 m u ( r ) + [ q 0 , 0 + q 0 , h r h ] r m u ( r ) = 0 ,
where, for m Z , h Z + ,
q j ( r ) k = 0 q j , k r k = p j ( r ) r m j , j = 0 , 1 , 2 ,
and at least one of the product q j , 0 q i , h , for i , j = 0 , 1 , 2 , is nonzero. The two-term recurrence formula is given by
c k = q 2 , h ( k + λ h ) ( k + λ h 1 ) + q 1 , h ( k + λ h ) + q 0 , h q 2 , 0 ( k + λ ) ( k + λ 1 ) + q 1 , 0 ( k + λ ) + q 0 , 0 c k h , c 0 0 ,
where λ = λ 1 , λ 2 are the roots of the indicial equation
q 2 , 0 λ ( λ 1 ) + q 1 , 0 λ + q 0 , 0 = 0 .
The closed form of the series solution can be written in terms of the generalized hypergeometric function as
u ( r ; λ ) = r λ k = 0 c h k r h k = r λ 3 F 2 1 , 2 λ 1 2 h + q 1 , h 2 h q 2 , h ( q 1 , h q 2 , h ) 2 4 q 0 , h q 2 , h 2 h q 2 , h , 2 λ 1 2 h + q 1 , h 2 h q 2 , h + ( q 1 , h q 2 , h ) 2 4 q 0 , h q 2 , h 2 h q 2 , h ; 1 + 2 λ 1 2 h + q 1 , 0 2 h q 2 , 0 ( q 1 , 0 q 2 , h ) 2 4 q 0 , 0 q 2 , 0 2 h q 2 , 0 , 1 + 2 λ 1 2 h + q 1 , 0 2 h q 2 , 0 + ( q 1 , 0 q 2 , h ) 2 4 q 0 , 0 q 2 , 0 2 h q 2 , 0 ; q 2 , h q 2 , 0 r h .
As an example, we consider the differential equation
( α 2 r 2 + α 1 r + α 0 ) u ( r ) + ( β 1 r + β 0 ) u ( r ) + γ 0 u ( r ) = 0 ,
mentioned in the introduction. We want to extract the possible equations that have two-term recurrence relations between successive coefficients in their series solutions.
Comparing Equation (25) with (20), we find that h + 2 m = 2 and h + 2 m = 2 > 2 m 0 , which leads to h = m , 2 > 2 m 0 . Consequently 0 < m 2 which gives us two cases: ( m , h ) = ( 1 , 1 ) , and ( m , h ) = ( 2 , 2 ) .
  • In the case of m = 1 and h = 1 , Equation (20) reads:
    [ q 2 , 0 + q 2 , 1 r ] r u ( r ) + [ q 1 , 0 + q 1 , 1 r ] u ( r ) + [ q 0 , 0 + q 0 , 1 r ] r 1 u ( r ) = 0 ,
    and the corresponding Equation (25) with q 2 , 0 = α 1 , q 2 , h = 1 = α 2 , q 1 , 0 = β 0 , q 1 , h = 1 = β 1 , q 0 , h = 1 = γ 0 and q 0 , 0 = 0 reads
    ( α 1 r + α 2 r 2 ) u ( r ) + ( β 0 + β 1 r ) u ( r ) + γ 0 u ( r ) = 0 .
    The two-term recurrence formula for the solution of Equation (27) is
    c k = α 2 ( k + λ 1 ) ( k + λ 2 ) + β 1 ( k + λ 1 ) + γ 0 α 1 ( k + λ ) ( k + λ 1 ) + β 0 ( k + λ ) c k 1 , c 0 = 1 ,
    where λ are the zeros of the indicial equation α 1 λ ( λ 1 ) + β 0 λ = 0 , namely, λ = 0 , and λ = 1 β 0 / α 1 .
    The closed form solutions in terms of the Gauss hypergeometric functions, respectively, are:
    u 1 ( r , λ = 0 ) = 2 F 1 1 2 + β 1 2 α 2 ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 , 1 2 + β 1 2 α 2 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 ; β 0 α 1 ; α 2 α 1 r ,
    and
    u 2 ( r , λ = 1 β 0 α 1 ) = r 1 β 0 α 1 × 2 F 1 1 2 + β 1 2 α 2 β 0 α 1 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 1 α 2 , 1 2 + β 1 2 α 2 + β 0 α 1 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 1 α 2 ; 2 β 0 α 1 ; α 2 α 1 r .
  • In this case of m = 2 and h = 2 , Equation (20) reads
    [ q 2 , 0 + q 2 , h r 2 ] u ( r ) + [ q 1 , 0 + q 1 , h r 2 ] r 1 u ( r ) + [ q 0 , 0 + q 0 , h r 2 ] r 2 u ( r ) = 0 ,
    and the corresponding Equation (25) reads
    ( α 0 + α 2 r 2 ) u ( r ) + β 1 r u ( r ) + γ 0 u ( r ) = 0 .
    The two-term recurrence relation for the solutions of Equation (32) reads
    c k = α 2 ( k + λ 2 ) ( k + λ 3 ) + β 1 ( k + λ 2 ) + γ 0 α 0 ( k + λ ) ( k + λ 1 ) c k 2 , k 2 ,
    where, λ are the roots of the indicial equation λ ( λ 1 ) = 0 , α 0 0 , with closed form solutions
    u 1 ( r , λ = 0 ) = 2 F 1 β 1 4 α 2 1 4 ( α 2 β 1 ) 2 4 α 2 γ 0 4 α 2 , β 1 4 α 2 1 4 + ( α 2 β 1 ) 2 4 α 2 γ 0 4 α 2 ; 1 2 ; α 2 α 0 r 2 ,
    and
    u 2 ( r , λ = 1 ) = r 2 F 1 β 1 4 α 2 + 1 4 ( α 2 β 1 ) 2 4 α 2 γ 0 4 α 2 , β 1 4 α 2 + 1 4 + ( α 2 β 1 ) 2 4 α 2 γ 0 4 α 2 ; 3 2 ; α 2 α 0 r 2 .
Remark 4.
For the admissible values of the equation parameters, the recurrence relations (28) and (33) are generic formulas for the two-term recurrence relations for the series solutions of the differential Equations (27) and (32) respectively. That is to say, by assigning admissible values of the parameters, (27) and (32) will generate solvable differential equations.
Remark 5.
Polynomial solutions for the differential Equations (27) and (32) can be easily obtained using suitable values of the equation parameters to terminate the recurrence relations (28) and (33).
As a second example we consider the differential equation
α 3 r 3 + α 2 r 2 + α 1 r + α 0 y ( r ) + β 2 r 2 + β 1 r + β 0 y ( r ) + γ 1 r + γ 0 y ( r ) = 0 .
As before, we want to extract the equations with a two-term recurrence relationship between successive coefficients in its series solution. We note by comparing Equation (36) and Equation (20) that h + 2 m = 3 and 3 > 2 m 0 , so 1 < m 2 . Three cases ( m , h ) = ( 0 , 1 ) , ( m , h ) = ( 1 , 2 ) and ( m , h ) = ( 2 , 3 ) follow from these constraints, and we consider each of these cases separately:
  • In the case where m = 0 and h = 1 , Equation (20) reads, with q 2 , h = 2 = α 3 , q 1 , h = 2 = β 2 , q 0 , h = 2 = γ 1 , q 2 , 0 = α 2 , q 1 , 0 = β 1 , and q 0 , 0 = γ 0
    [ α 2 r 2 + α 3 r 3 ] u ( r ) + [ β 1 r + β 2 r 2 ] u ( r ) + ( γ 0 + γ 1 r ) u ( r ) = 0 , γ 0 0 ,
    For this differential equation, the recurrence relation (22) implies
    c k = α 3 ( k + λ 1 ) ( k + λ 2 ) + β 2 ( k + λ 1 ) + γ 1 α 2 ( k + λ ) ( k + λ 1 ) + β 1 ( k + λ ) + γ 0 c k 1 , c 0 0 ,
    where λ are the roots of the indicial equation α 2 λ ( λ 1 ) + β 1 λ + γ 0 = 0 . The values of λ with (38) gives the following two linearly independent series solutions for the differential Equation (37):
    u 1 ( r ) = r α 2 β 1 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 × 2 F 1 β 2 2 α 3 β 1 2 α 2 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 ( α 3 β 2 ) 2 4 α 3 γ 1 2 α 3 , β 2 2 α 3 β 1 2 α 2 + ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 + ( α 3 β 2 ) 2 4 α 3 γ 1 2 α 3 ; α 2 + ( α 2 β 1 ) 2 4 α 2 γ 0 α 2 ; α 3 α 2 r ,
    and
    u 2 ( r ) = r α 2 β 1 ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 × 2 F 1 β 2 2 α 3 β 1 2 α 2 ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 ( α 3 β 2 ) 2 4 α 3 γ 1 2 α 3 , β 2 2 α 3 β 1 2 α 2 ( α 2 β 1 ) 2 4 α 2 γ 0 2 α 2 + ( α 3 β 2 ) 2 4 α 3 γ 1 2 α 3 ; α 2 ( α 2 β 1 ) 2 4 α 2 γ 0 α 2 ; α 3 α 2 r .
  • In the case where m = 1 and h = 2 , Equation (20) reads,
    [ α 1 r + α 3 r 3 ] u ( r ) + [ β 0 + β 2 r 2 ] u ( r ) + γ 1 r u ( r ) = 0 , β 0 0 ,
    For this differential equation, the recurrence relation (22) implies
    c k = α 3 ( k + λ 2 ) ( k + λ 3 ) + β 2 ( k + λ 2 ) + γ 1 α 1 ( k + λ ) ( k + λ 1 ) + β 0 ( k + λ ) c k 2 ,
    where λ are the roots of the indicial equation a 1 λ ( λ 1 ) + b 0 λ = 0 . The values of λ with (42) gives the following two linearly independent series solutions for the differential Equation (41):
    u 1 ( r ) = 2 F 1 β 2 4 α 3 1 4 ( α 3 β 2 ) 2 4 α 3 γ 1 4 α 3 , β 2 4 α 3 1 4 + ( α 3 β 2 ) 2 4 α 3 γ 1 4 α 3 ; α 1 + β 0 2 α 1 ; α 3 r 2 α 1 ,
    and
    u 2 ( r ) = r 1 β 0 α 1 ×
    2 F 1 1 4 + β 2 4 α 3 β 0 2 α 1 ( α 3 β 2 ) 2 4 α 3 γ 1 4 α 3 , 1 4 + β 2 4 α 3 β 0 2 α 1 + ( α 3 β 2 ) 2 4 α 3 γ 1 4 α 3 ; 3 2 β 0 2 α 1 ; α 3 r 2 α 1 .
  • In the case where m = 2 and h = 3 , Equation (20) reads,
    ( α 0 + α 3 r 3 ) u ( r ) + β 2 r 2 u ( r ) + γ 1 r u ( r ) = 0 , α 0 0 ,
    For this differential equation, the recurrence relation (22) implies
    c k = α 3 ( k + λ 3 ) ( k + λ 4 ) + β 2 ( k + λ 3 ) + γ 1 α 0 ( k + λ ) ( k + λ 1 ) c k 3 , c 0 0 ,
    where λ are the roots of the indicial equation α 0 λ ( λ 1 ) = 0 . The values of λ with (47) gives the following two linearly independent series solutions for the differential Equation (46):
    u 1 ( r ) = 2 F 1 ( α 3 β 2 ) 2 4 α 3 γ 1 + α 3 β 2 6 α 3 , ( α 3 β 2 ) 2 4 α 3 γ 1 α 3 + β 2 6 α 3 ; 2 3 ; α 3 α 0 r 3 ,
    and
    u 2 ( r ) = r 2 F 1 ( α 3 β 2 ) 2 4 α 3 γ 1 + α 3 + β 2 6 α 3 , ( α 3 β 2 ) 2 4 α 3 γ 1 + α 3 + β 2 6 α 3 ; 4 3 ; α 3 α 0 r 3 .
Remark 6.
When the coefficients of the differential Equation (18) are polynomials, there are n generic equations with series solutions explicitly found using the two-term recurrence relation (22). The closed-form solutions of these equations are given by (24). Indeed, for h + 2 m = n , n > 2 m 0 , we have n cases given h = n + m 2 , 2 n < m 2 , n 2 . These n equations are as follows (See Theorem 3):
( α 0 + α n r n ) y + β n 1 r n 1 y + γ n 2 r n 2 y = 0 , α 0 0 ( α 1 r + α n r n ) y + ( β 0 + β n 1 r n 1 ) y + γ n 2 r n 2 y = 0 , β 0 0 ( α 2 r 2 + α n r n ) y + ( β 1 r + β n 1 r n 1 ) y + ( γ 0 + γ n 2 r n 2 ) y = 0 , γ 0 0 ( α 3 r 3 + α n r n ) y + ( β 2 r 2 + β n 1 r n 1 ) y + ( γ 1 r + γ n 2 r n 2 ) y = 0 , ( α 4 r 4 + α n r n ) y + ( β 3 r 3 + β n 1 r n 1 ) y + ( γ 2 r 2 + γ n 2 r n 2 ) y = 0 , ( α n 1 r n 1 + α n r n ) y + ( β n 2 r n 2 + β n 1 r n 1 ) y + ( γ n 3 r n 3 + γ n 2 r n 2 ) y = 0 .

4. Theorems and Algorithm

In this section we present three elementary theorems that classify the polynomial solutions of the general differential Equation (3). Following that, we will give a brief description of the Mathematica® program that was developed to accompany the present work. A link to the complete program, available for direct use, is also provided for direct applications of these theorems.

4.1. Theorems

Theorem 4 applies to 2 n equations of (3) and gives a recurrence relation needed to compute polynomial solutions in the neighbourhood of an ordinary point. Theorems 5 and 6 give the polynomial solutions for the 2 n 1 + 2 n 2 equations from (3) about a regular singular point. The proofs of Theorems 4 and 5 may be found in the Appendix A, while the proof of Theorem 6 follows closely to the proof of Theorem 5 and is therefore omitted.
Theorem 4.
( α 0 0 )
For α 0 0 , the second-order linear differential Equation (3) admits the solution y ( r ) = P m ( r ) = i = 0 m C i r i in the neighbourhood of the ordinary point r = 0 and valid to the nearest nonzero singular point of P m ( r ) = 0 . Here P m ( r ) is an mth-degree polynomial if
j = 0 n α n j ( n + j + 2 ) ( n + j + 1 ) + β n j 1 ( n + j + 2 ) τ n j 2 C n + j + 2 = 0 ,
for all = 0 , 1 , 2 , , m + n 2 , and τ i = β i = 0 if i < 0 . The equation corresponding to the value of = m + n 2 gives the necessary condition (7) for the existence of the mth-degree polynomial solution
τ n 2 = α n m ( m 1 ) + β n 1 m , m = 0 , 1 , 2 , .
The remaining ( m + n 2 ) linear equations consist of ( n 2 ) sufficient conditions that relate { α k } k = 0 n and { β k } k = 0 n 1 with { τ k } k = 0 n 2 , in addition to the m linear equations required to evaluate the coefficients { C i } i = 1 m of the polynomial solution in terms of C 0 0 .
Theorem 5.
( α 0 = 0 , α 1 0 )
For α 0 = 0 and α 1 0 , the second-order linear differential Equation (3) admits the solution y ( r ) = r s P m ( r ) = r s i = 0 m C i r i in the neighbourhood of the regular singular point r = 0 , where P m ( r ) is an mth-degree polynomial if
j = 0 n 1 α n j ( n + j + s + 2 ) ( n + j + s + 1 ) + β n j 1 ( n + j + s + 2 ) τ n j 2 C n + j + 2 = 0 ,
for all = 0 , 1 , 2 , , m + n 2 , and τ i = 0 if i < 0 . Here s is a root of the indicial equation: s ( s 1 ) + ( β 0 / α 1 ) s = 0 , so s 0 , 1 ( β 0 / α 1 ) .
Theorem 6.
( α 0 = α 1 = β 0 = 0 , α 2 0 )
For α 0 = α 1 = β 0 = 0 and α 2 0 , the second-order linear differential Equation (3) admits the solution y ( r ) = r s P m ( r ) = r s i = 0 m C i r i in the neighbourhood of the regular singular point r = 0 , where P m ( r ) is an mth-degree polynomial if
j = 0 n 2 α n j ( n + j + s + 2 ) ( n + j + s + 1 ) + β n j 1 ( n + j + s + 2 ) τ n j 2 C n + j + 2 = 0 ,
for all = 0 , 1 , 2 , , m + n 2 . Here s is a nonzero root of the indicial equation: s ( s 1 ) + ( β 0 / α 1 ) s τ 0 α 2 = 0 , so s α 2 β 1 ± ( α 2 β 1 ) 2 + 4 α 2 τ 0 / 2 α 2 .
Remark 7.
(Cauchy-Euler equation): In the case of n = 2 , Equation (3) reduces to the equation,
α 2 r 2 y ( r ) + β 1 r y ( r ) τ 0 y ( r ) = 0 ,
for which the recurrence relation (54) reduces to:
α 2 ( + s ) ( + s 1 ) + β 1 ( + s ) τ 0 C = 0 , = 0 , 1 , , m .
Clearly when = m and if C m 0 then Equation (56) gives the necessary condition:
α 2 ( m + s ) ( m + s 1 ) + β 1 ( m + s ) τ 0 = 0 , C m 0 ,
as expected for the solutions of the Cauchy-Euler equation.

4.2. The Mathematica® Program

A Mathematica® algorithm and accompanying GUI (Figure 1) were developed to assist with this work. The algorithm uses Cramer’s rule to compute the coefficients of the polynomial solutions { C i } i = 1 m for the Equation (3) in the case where the equation is in the form of Theorems 4–6. In addition to computing the coefficients, the program also generates the sufficient and necessary conditions. A complete Mathematica® program may be obtained by contacting the authors, or via the following Github repository: https://github.com/KyleBryenton/OPSOLDE.

5. Applications

The simplicity of the above theorems lay in their constructive approach to explicitly computing the polynomial solution coefficients and to provide all the necessary and sufficient conditions for the existence of these polynomials. Further, they can be easily implemented using any available symbolic algebra system. In this section, we show how to make use of these theorems.

5.1. General Case ( n = 4 )

To illustrate the constructive approach of Theorem 4, we consider Equation (3) in the case of n = 4 which reads:
( α 4 r 4 + α 3 r 3 + α 2 r 2 + α 1 r + α 0 ) y ( r ) + ( β 3 r 3 + β 2 r 2 + β 1 r + β 0 ) y ( r ) ( τ 2 r 2 + τ 1 r + τ 0 ) y ( r ) = 0 , where α 0 0 , α k , β k , τ k R k .
For the mth-degree polynomial solutions P m ( r ) y ( r ) = i = 0 m C i r i , (51) gives the five term recurrence-relation:
[ α 4 ( 2 ) ( 3 ) + β 3 ( 2 ) τ 2 ] C 2 + [ α 3 ( 1 ) ( 2 ) + β 2 ( 1 ) τ 1 ] C 1 + [ α 2 ( 1 ) + β 1 τ 0 ] C + [ α 1 ( + 1 ) + β 0 ( + 1 ) ] C + 1 + [ α 0 ( + 2 ) ( + 1 ) ] C + 2 = 0 ,
for = 0 , 1 , 2 , , m + 2 and C j = 0 for j < 0 .
For the zero-degree ( m = 0 ) polynomial solution, P 0 ( r ) = 1 , = 0 , 1 , 2 . Noting that C j = 0 , for all j 0 , it follows for = 0 and = 1 that τ 0 = 0 and τ 1 = 0 respectively, while for = 2 , the necessary condition τ 2 = 0 follows.
For the first-degree ( m = 1 ) polynomial solution, P 1 ( r ) = C 0 + C 1 r , = 0 , 1 , 2 , 3 . Noting that C j = 0 , for all j 0 , 1 , we have for = 0 that τ 0 C 0 + β 0 C 1 = 0 , while for = 1 and = 2 we have τ 1 C 0 + ( β 1 τ 0 ) C 1 = 0 and τ 2 C 0 + ( β 2 τ 1 ) C 1 = 0 respectively. For = 3 it follows ( β 3 τ 2 ) C 1 = 0 from which it is necessary that τ 2 = β 3 , while the first equation ( = 0 ) gives C 1 = τ 0 C 0 / β 0 . Using this value for C 1 and letting C 0 = 1 , the second and third equations ( = 1 , 2 ) give τ 1 β 0 + ( β 1 τ 0 ) τ 0 = 0 and τ 2 β 0 + ( β 2 τ 1 ) τ 0 = 0 respectively. We can write these conditions compactly as:
= 0 : C 1 = τ 0 / β 0 , = 1 , 2 : τ 0 β 0 τ 1 β 1 τ 0 = 0 , ( S 1 ) τ 0 β 0 τ 2 β 2 τ 1 = 0 , ( S 2 ) = 3 : τ 2 = β 3 , ( NC )
where (S1) and (S2) refer to the two sufficient conditions and (NC) refers to the necessary condition.
For the second-degree ( m = 2 ) polynomial solution P 2 ( r ) = 1 + C 1 r + C 2 r 2 , noting that C j = 0 , for all j 0 , 1 , 2 , where = 0 , 1 , 2 , 3 , 4 now yields
= 0 , 1 : C 1 = τ 0 2 α 0 τ 1 2 α 1 + β 0 β 0 2 α 0 β 1 τ 0 2 α 1 + β 0 , C 2 = β 0 τ 0 β 1 τ 0 τ 1 β 0 2 α 0 β 1 τ 0 2 α 1 + β 0 , = 2 , 3 : τ 0 β 0 2 α 0 τ 1 β 1 τ 0 2 α 1 + β 0 τ 2 β 2 τ 1 2 α 2 + β 1 τ 0 = 0 , ( S 1 ) τ 0 β 0 2 α 0 τ 1 β 1 τ 0 2 α 1 + β 0 0 β 3 τ 2 2 α 3 + β 2 τ 1 = 0 , ( S 2 ) = 4 : 2 α 4 + 2 β 3 τ 2 = 0 . ( NC )
For the third degree ( m = 3 ) polynomial solution P 3 ( r ) = 1 + C 1 r + C 2 r 2 + C 3 r 3 , noting that C j = 0 , for all j 0 , 1 , 2 , 3 , we have for = 0 , 1 , 2 , 3 , 4 , 5
= 0 , 1 , 2 : C 1 = τ 0 2 α 0 0 τ 1 2 α 1 + β 0 6 α 0 τ 2 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 β 0 2 α 0 0 β 1 τ 0 2 α 1 + β 0 6 α 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 , C 2 = β 0 τ 0 0 β 1 τ 0 τ 1 6 α 0 β 2 τ 1 τ 2 3 2 α 1 + β 0 β 0 2 α 0 0 β 1 τ 0 2 α 1 + β 0 6 α 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 , C 3 = β 0 2 α 0 τ 0 β 1 τ 0 2 α 1 + β 0 τ 1 β 2 τ 1 2 α 2 + β 1 τ 0 τ 2 β 0 2 α 0 0 β 1 τ 0 2 α 1 + β 0 6 α 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 ,
= 3 , 4 : τ 0 β 0 2 α 0 0 τ 1 β 1 τ 0 2 α 1 + β 0 6 α 0 τ 2 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 0 β 3 τ 2 2 α 3 + β 2 τ 1 6 α 2 + 3 β 1 τ 0 = 0 , ( S 1 ) τ 0 β 0 2 α 0 0 τ 1 β 1 τ 0 2 α 1 + β 0 6 α 0 τ 2 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 0 0 2 α 4 + β 3 τ 2 6 α 3 + 3 β 2 τ 1 = 0 , ( S 2 )
= 5 : 6 α 4 + 3 β 3 τ 2 = 0 . ( NC )
Higher-order polynomial solutions may be obtained similarly in this constructive manner.

5.2. Two-Dimensional Dirac Equation

The general covariant ( 2 + 1 ) -dimensional Dirac equation reads
i γ μ ( μ + i e A μ ) Ψ ( t , r ) = M Ψ ( t , r ) ,
where M is the rest mass of the particle, A μ = ( A 0 , A ) is the external electromagnetic field potential, and e is the charge of the particle. Here, γ 0 = σ z , γ 1 = i σ x , γ 2 = i σ y , where σ x , σ y , and σ z are the Pauli spin matrices. Using the generalized momentum operator P k = i k + e A k for k = 1 , 2 and the conventional summation over μ , Equation (63) may be written in the form
( i t σ x P y + σ y P x σ z m e A 0 ) Ψ ( t , r ) = 0 .
Chiang et al. [16] considered an electron moving in a Coulombic field, A 0 = Z e / r , with a constant homogeneous magnetic field described by A 1 = B y / 2 and A 2 = B x / 2 . With this choice of potential and for r > 0 , (64) becomes
i t M e i ϕ ( r ϕ + e B r 2 ) e i ϕ ( r + ϕ e B r 2 ) M Z e 2 r Ψ ( t , r ) = 0 ,
where polar coordinates, the identities P x ± i P y = i e ± i ϕ ( r ± ( ϕ ( e B r / 2 ) ) ) , x = cos ϕ r [ sin ( ϕ ) / r ] ϕ , and y = sin ( ϕ ) r + [ cos ( ϕ ) / r ] ϕ were employed [16]. Using separation of variables
Ψ ( t , r ) = 1 r e i E t F ( r ) e i l ϕ G ( r ) e i ( l + 1 ) ϕ ,
where l Z , the Dirac radial Equation (65) reads
d F ( r ) d r l + 1 2 r + e B 2 r F ( r ) + E + M + Z e 2 r G ( r ) = 0 ,
d G ( r ) d r + l + 1 2 r + e B 2 r G ( r ) E M + Z e 2 r F ( r ) = 0 ,
For a strong magnetic field, using the asymptotic solutions near r = 0 and r = , it is beneficial to assume an ansatz [16]
F ( r ) = r γ e e B r 2 / 4 Q ( r ) , G ( r ) = r γ e e B r 2 / 4 P ( r ) ,
where γ = ( l + 1 / 2 ) 2 ( Z e 2 ) 2 . Substituting (69) into (67) and (68) yields
d Q ( r ) d r + ( γ l 1 2 e B r 2 ) r Q ( r ) + E + M + Z e 2 r P ( r ) = 0 ,
d P ( r ) d r + ( γ + l + 1 2 ) r P ( r ) E M + Z e 2 r Q ( r ) = 0 .
Solving (70) for P ( r ) and substituting the resulting equation into (71) forms a second-order differential equation for Q ( r )
r ( r + r 0 ) Q ( r ) ( e B r 3 + e B r 0 r 2 β 1 r β r 0 ) Q ( r ) + ( ( e B + ε ) r 2 + ( ε r 0 + a ) r + c a r 0 ) Q ( r ) = 0
where r 0 Z e 2 / ( E + M ) , , β 2 γ + 1 , a 2 E Z e 2 , ε E 2 M 2 e B ( l + 5 2 + γ ) , c γ l ( 1 / 2 ) . Equation (72) is an example of differential equations scattered over a vast range of applications in applied mathematics and theoretical physics [16,17,18,19]. The coefficients { C i } i = 0 m of the polynomial solutions Q m ( r ) = i = 0 m C i r i satisfy the following recurrence relation (53): for = 0 , 1 , , m + 2
E 2 M 2 B e l + γ + 1 2 C 2 + Z e 2 ( 3 E M ) B e 3 Z ( l + γ + + 3 2 ) E + M C 1 + ( 2 + 1 ) γ + ( 1 ) l a e 2 Z E + M 1 2 C + e 2 Z ( 1 + ) ( 1 + 2 γ + ) E + M C + 1 = 0 .
From Theorem 5, the necessary condition for the polynomial solutions of (72) is
ε = e B ( m 1 ) E m 2 = e B m + l + γ + 3 2 + M 2 , m = 0 , 1 , .
For the zero-degree polynomial solution Q 0 ( r ) = 1 where C 0 = 1 , m = 0 and therefore = 0 , 1 , 2 , it follows that:
= 0 , 1 : ( 1 + 2 l 2 γ ) ( E + M ) + 4 e 4 Z 2 E = 0 , ( S 1 ) e B l + γ + 5 2 ( 3 E M ) ( E + M ) = 0 , ( S 2 )
= 2 : E 2 = B e l + γ + 3 2 + M 2 . ( NC )
For the first-degree polynomial solution Q 1 ( r ) = 1 + C 1 r where C 0 = 1 and C 1 0 , m = 1 and therefore = 0 , 1 , 2 , 3 , it follows that:
= 0 : C 1 = ( 1 + 2 l 2 γ ) ( E + M ) + 4 e 4 E Z 2 2 e 2 Z ( 1 + 2 γ ) , = 1 , 2 : e 2 Z B e ( 5 + 2 l + 2 γ ) 2 ( 3 E M ) ( E + M ) + ( 1 + 2 l 6 γ ) ( E + M ) + 4 e 4 E Z 2 C 1 = 0 , ( S 1 ) E 2 B e l + γ + 3 2 M 2 + e 2 Z ( 3 E M ) B e 3 Z ( l + γ + 7 2 ) E + M C 1 = 0 , ( S 2 ) = 3 : E 2 = e B l + γ + 5 2 + M 2 . ( NC )
For the second-degree polynomial solution Q 2 ( r ) = 1 + C 1 r + C 2 r 2 where C 0 = 1 and C 2 0 , m = 2 and therefore = 0 , 1 , 2 , 3 , 4 , it follows that:
= 0 , 1 : C 1 = ( 1 + 2 l 2 γ ) ( E + M ) + 4 e 4 E Z 2 2 e 2 ( 1 + 2 γ ) Z , l 3 γ + 1 2 ( E + M ) + 2 e 4 Z 2 E C 1 4 e 2 Z ( 1 + γ ) C 2 + e 2 Z B e l + γ + 5 2 ( 3 E M ) ( E + M ) = 0 ,
= 2 , 3 : ( E + M ) E 2 M 2 B e l + γ + 3 2 C 0 + e 2 Z ( 3 E M ) ( E + M ) B e l + γ + 7 2 C 1 + 3 2 l + 5 γ ( E + M ) 2 e 4 E Z 2 C 2 = 0 , ( S 1 ) E 2 M 2 B e l + γ + 5 2 C 1 e 2 Z B e ( l + γ + 9 2 ) ( 3 E M ) ( E + M ) E + M C 2 = 0 , ( S 2 ) = 4 : E 2 = e B l + γ + 7 2 + M 2 . ( NC )
Higher-order polynomial solutions follows similarly.

5.3. The Heun Equation

Over the past two decades Heun’s equation has attracted considerable attention due to its increasing number of applications in applied mathematics and theoretical physics [3,4,5,6,7,8,41,42,43,55]. The general second-order Heun differential equation can be written as
α 3 r 3 + α 2 r 2 + α 1 r y ( r ) + β 2 r 2 + β 1 r + β 0 y ( r ) τ 1 r + τ 0 y ( r ) = 0 ,
where the canonical form [41] may be obtained using the substitutions given in Table 1.
We immediately see that because α 0 = 0 and α 1 0 that (78) falls under Theorem 5. By Equation (53) we obtain the three-term recurrence relation
[ α 3 ( + s 2 ) ( + s 1 ) + β 2 ( + s 1 ) τ 1 ] C 1 + α 2 ( + s 1 ) ( + s ) + β 1 ( + s ) τ 0 C + α 1 ( + s ) ( + s + 1 ) + β 0 ( + s + 1 ) C + 1 = 0 ,
for = 0 , 1 , 2 , , m + 1 and s 0 , 1 ( β 0 / α 1 ) . We may then input this into the Mathematica® program or work directly with the ( n + 2 ) equations generated by the recurrence relations (79). Using either approach with s = 0 will generate the results, as follows.
For the zero-degree ( m = 0 ) polynomial solution, P 0 ( r ) = 1 , = 0 , 1 . Noting that C j = 0 for all j 0 , it follows for = 0 that τ 0 = 0 , while for = 1 , τ 1 = 0 .
For the first-degree ( m = 1 ) polynomial solution, P 1 ( r ) = C 0 + C 1 r , = 0 , 1 , 2 . It follows that
= 0 : C 1 = τ 0 / β 0 = 1 : τ 0 β 0 τ 1 β 1 τ 0 = 0 , ( SC ) = 2 : τ 1 = β 2 . ( NC )
where (SC) refers to the sufficient condition and (NC) refers to the necessary condition, and we have set C 0 = 1 for convenience.
For the second-degree ( m = 2 ) polynomial solution P 2 ( r ) = 1 + C 1 r + C 2 r 2 , = 0 , 1 , 2 , 3 yields
= 0 , 1 : C 1 = τ 0 0 τ 1 2 α 1 + β 0 β 0 0 β 1 τ 0 2 α 1 + β 0 , C 2 = β 0 τ 0 β 1 τ 0 τ 1 β 0 0 β 1 τ 0 2 α 1 + β 0 , = 2 : τ 0 β 0 0 τ 1 β 1 τ 0 2 α 1 + β 0 0 β 2 τ 1 2 α 2 + β 1 τ 0 = 0 , ( SC )
= 3 : 2 α 3 + 2 β 2 τ 1 = 0 . ( NC )
For the third degree ( m = 3 ) polynomial solution P 3 ( r ) = 1 + C 1 r + C 2 r 2 + C 3 r 3 . Noting that C j = 0 , for all j 0 , 1 , 2 , we have for = 0 , 1 , 2 , 3 , 4
= 0 , 1 , 2 : C 1 = τ 0 0 0 τ 1 2 α 1 + β 0 0 0 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 β 0 0 0 β 1 τ 0 2 α 1 + β 0 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 , C 2 = β 0 τ 0 0 β 1 τ 0 τ 1 0 β 2 τ 1 0 3 2 α 1 + β 0 β 0 0 0 β 1 τ 0 2 α 1 + β 0 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 , C 3 = β 0 0 τ 0 β 1 τ 0 2 α 1 + β 0 τ 1 β 2 τ 1 2 α 2 + β 1 τ 0 0 β 0 0 0 β 1 τ 0 2 α 1 + β 0 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 ,
= 3 : τ 0 β 0 0 0 τ 1 β 1 τ 0 2 α 1 + β 0 0 0 β 2 τ 1 2 α 2 + β 1 τ 0 3 2 α 1 + β 0 0 0 2 α 3 + β 2 τ 1 6 α 2 + 3 β 1 τ 0 = 0 , ( SC ) = 4 : 6 α 3 + 3 β 2 τ 1 = 0 . ( NC )
This approach allows us to construct all the Heun polynomial solutions in a constructive way that is not available in the literature.

6. Conclusions

Given a linear differential equation with polynomial coefficients of order n, n 1 , and n 2 respectively, this work provides a constructive and straightforward approach for finding all possible polynomial solutions along with the necessary and sufficient conditions. A Mathematica program which solves these differential equations for an integer value of n [ 2 , 9 ] is available for download (see Section 4.2), but the algorithms are also provided so that the reader may adopt them in their preferred computer software. Although there may be other sophisticated approaches available in the literature for solving differential equations [46], the method presented here has its simplicity characterized by recurrence relations. Aside from the examples discussed in this present work, the authors explored a large number of linear differential equations which appeared in mathematics and physics literature, and have found exact consistency with the results obtained by other researchers [3,4,5,6,7,8,9,10,43,55].

Author Contributions

Conceptualization, N.S.; data curation, P.S. and N.V.; formal analysis, A.R.C. and K.L.A.K.; investigation, P.S and N.V; methodology, K.L.A.K. and N.S.; project administration, N.S; resources, K.R.B., N.S. and P.S.; software, K.R.B.; validation, K.L.A.K. and N.S.; visualization, A.R.C., P.S. and N.V.; Writing—original draft, N.S.; writing— review and editing, K.R.B. and K.L.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

Partial financial support of this work, under Grant No. GP249507 from the Natural Sciences and Engineering Research Council of Canada, is gratefully acknowledged.

Acknowledgments

We would like to extend thanks to Le Phuong Uyen Dao for her helpful correspondence pertaining to Section 3 of this article. We would also like to thank the reviewers for their helpful comments and suggestions to improve the present work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Theorem 4.
For α 0 0 the rational functions p ( r ) = k = 0 n 1 β k r k / k = 0 n α k r k and q ( r ) = k = 0 n 2 τ k r k / k = 0 n α k r k are analytic at r = 0 . Consequently, r = 0 is an ordinary point of the differential equation given by (3) with a power series solution of the form
y ( r ) = = 0 C r , y ( r ) = = 1 C r 1 , y ( r ) = = 2 C ( 1 ) r 2 .
Substituting y ( r ) , y ( r ) , and y ( r ) from (A1) into (3) yields
k = 0 n = 2 α k C ( 1 ) r k + 2 + k = 0 n 1 = 1 β k C r k + 1 k = 0 n 2 = 0 τ k C r k + = 0 .
We then shift the -summation to arrive at
k = 0 n = 2 α k C ( 1 ) r k + 2 + k = 0 n 1 = 2 β k C 1 ( 1 ) r k + 2 k = 0 n 2 = 2 τ k C 2 r k + 2 = 0 .
Now we group these into common degrees of r
= 2 α 0 C ( 1 ) + β 0 C 1 ( 1 ) τ 0 C 2 r 2 + = 2 α 1 C ( 1 ) + β 1 C 1 ( 1 ) τ 1 C 2 r 1 + = 2 α 2 C ( 1 ) + β 2 C 1 ( 1 ) τ 2 C 2 r + + = 2 α n 2 C ( 1 ) + β n 2 C 1 ( 1 ) τ n 2 C 2 r + n 4 + = 2 α n 1 C ( 1 ) + β n 1 C 1 ( 1 ) r + n 3 + = 2 α n C ( 1 ) r + n 2 = 0 .
Then, shifting the individual sums obtain a common degree of r results in
+ n = 2 α 0 C + n ( + n ) ( + n 1 ) + β 0 C + n 1 ( + n 1 ) τ 0 C + n 2 r + n 2 + + n 1 = 2 α 1 C + n 1 ( + n 1 ) ( + n 2 ) + β 1 C + n 2 ( + n 2 ) τ 1 C + n 3 r + n 2 + + n 2 = 2 α 2 C + n 2 ( + n 2 ) ( + n 3 ) + β 2 C + n 3 ( + n 3 ) τ 2 C + n 4 r + n 2 + + + 2 = 2 α n 2 C + 2 ( + 2 ) ( + 1 ) + β n 2 C + 1 ( + 1 ) τ n 2 C r + n 2 + + 1 = 2 α n 1 C + 1 ( + 1 ) + β n 1 C r + n 2 + = 2 α n C ( 1 ) r + n 2 = 0 .
To enforce polynomial solutions of degree m, we may modify the summations to go to an upper limit of m rather than . We wish to shift our bottom index by the degree of the polynomial on the τ terms, thus a shift of = n + 2 is required. The above formula exposes the ( n + 1 ) -term recurrence relation which may be grouped into common C j as
α 0 ( + 2 ) ( + 1 ) C + 2 + α 1 ( + 1 ) + β 0 ( + 1 ) C + 1 + α 2 ( 1 ) + β 1 τ 0 C + + α n ( n + 2 ) ( n + 1 ) + β n 1 ( n + 2 ) τ n 2 C n + 2 = 0 .
Finally, the mth-degree polynomial solution is given by
j = 0 n α n j ( n + j + 2 ) ( n + j + 1 ) + β n j 1 ( n + j + 2 ) τ n j 2 C n + j + 2 = 0 .
where = 0 , 1 , 2 , , m + n 2 , n 2 , and τ i = β i = 0 if i < 0 . □
Proof of Theorem 5.
For α 0 = 0 and α 1 0 the rational functions p ( r ) = k = 0 n 1 β k r k / k = 1 n α k r k and q ( r ) = k = 0 n 2 τ k r k / k = 1 n α k r k , are non-analytic at r = 0 . We see that
lim r 0 r k = 0 n 1 β k r k k = 1 n α k r k = β 0 α 1 , lim r 0 r 2 k = 0 n 2 τ k r k k = 1 n α k r k = 0 .
Thus the indicial roots are given by: s ( s 1 ) + ( β 0 / α 1 ) s = 0 , so s 0 , 1 ( β 0 / α 1 ) and r = 0 is a regular singular point of the differential Equation (3). The method of Frobenius assumes a formal series solution of (3) of the form
y = = 0 C r + s , y ( r ) = = 0 C ( + s ) r + s 1 , y ( r ) = = 0 C ( + s ) ( + s 1 ) r + s 2 .
Substituting y ( r ) , y ( r ) , and y ( r ) from (A9) into (3) yields
k = 1 n = 0 α k C ( + s ) ( + s 1 ) r k + + s 2 + k = 0 n 1 = 0 β k C ( + s ) r k + + s 1 k = 0 n 2 = 0 τ k C r k + + s = 0 .
We then group these into a common degree of r:
= 0 α 1 C ( + s ) ( + s 1 ) + β 0 C ( + s ) r + s 1 + = 0 α 2 C ( + s ) ( + s 1 ) + β 1 C ( + s ) τ 0 C r + s + + = 0 α n C ( + s ) ( + s 1 ) + β n 1 C ( + s ) τ n 2 C r n + + s 2 = 0 .
Now we shift the -summation to obtain r n + + s 2 in each term
n + 1 = 0 α 1 C n + 1 ( n + + s 1 ) ( n + + s 2 ) + β 0 C n + 1 ( n + + s 1 ) r n + + s 2 + n + 2 = 0 α 2 C n + 2 ( n + + s 2 ) ( n + + s 3 ) + β 1 C n + 2 ( n + + s 2 ) τ 0 C r n + + s 2 + + = 0 α n C ( + s ) ( + s 1 ) + β n 1 C ( + s ) τ n 2 C n + 2 r n + + s 2 = 0 .
To enforce polynomial solutions of degree m, we may modify the summations to go to an upper limit of m rather than . We wish to shift our bottom index by the degree of the polynomial on the τ terms, thus a shift of = n + 2 is required. The above formula exposes the n-term recurrence relation which may be grouped into common C j as
[ α 1 ( + s + 1 ) ( + s ) + β 0 ( + s + 1 ) ] C + 1 + α 2 ( + s ) ( + s 1 ) + β 1 ( + s ) τ 0 C + + α n ( n + s + 2 ) ( n + s + 1 ) + β n 1 ( n + s + 2 ) τ n 2 C n 2 = 0 .
Thus, the mth-degree polynomial solution is given by
j = 0 n 1 α n j ( n + j + s + 2 ) ( n + j + s + 1 ) + β n j 1 ( n + j + s + 2 ) τ n j 2 C n + j + 2 = 0 .
where = 0 , 1 , 2 , , m + n 2 , n 2 , and τ i = 0 if i < 0 . □

References

  1. Primitivo, B. Acosta-Humánez, David Blázquez-Sanz, Henock Venegas-Gómez, Liouvillian solutions for second order linear differential equations with polynomial coefficients. São Paulo J. Math. Sci. 2020. [Google Scholar] [CrossRef]
  2. Raposo, A.P.; Weber, H.J.; Alvarez-Castillo, D.E.; Kirchbach, M. Romanovski polynomials in selected physics problems. Centr. Eur. J. Phys. 2007, 5, 253–284. [Google Scholar] [CrossRef] [Green Version]
  3. Ovsiyuk, E.; Amirfachrian, M.; Veko, O. On Schrödinger equation with potential V(r)=αr−1+βr+kr2 and the biconfluent Heun functions theory. Nonl. Phen. Compl. Sys. 2012, 15, 163. [Google Scholar]
  4. Caruso, F.; Martins, J.; Oguri, V. Solving a two-electron quantum dot model in terms of polynomial solutions of a biconfluent Heun equation. Ann. Phys. 2014, 347, 130–140. [Google Scholar] [CrossRef] [Green Version]
  5. Marcilhacy, G.; Pons, R. The Schrödinger equation for the interaction potential x2 + λx2/(1 + gx2), and the first Heun confluent equation. J. Phys. A Math. Gen. 1985, 18, 2441–2449. [Google Scholar] [CrossRef]
  6. Exton, H. The interaction V(r) = −Ze2/(r + β) and the confluent Heun equation. J. Phys. A Math. Gen. 1991, 24, L329–L330. [Google Scholar] [CrossRef]
  7. Ciftci, H.; Hall, R.L.; Saad, N.; Dogu, E. Physical applications of second-order linear differential equations that admit polynomial solutions. J. Phys. A Math. Theor. 2010, 43, 415206. [Google Scholar] [CrossRef]
  8. Hall, R.H.; Saad, N.; Sen, K.D. Soft-core Coulomb potentials and Heun’s differential equation. J. Math. Phys. 2010, 51, 022107. [Google Scholar] [CrossRef] [Green Version]
  9. Hall, R.L.; Saad, N.; Sen, K.D. Discrete spectra for confined and unconfined −a/r + br2 potentials in d-dimensions. J. Math. Phys. 2011, 52, 092103. [Google Scholar] [CrossRef] [Green Version]
  10. Hall, R.L.; Saad, N.; Sen, K.D. Spectral characteristics for a spherically confined −a/r + br2 potential. J. Phys. A Math. Theor. 2011, 44, 185307. [Google Scholar] [CrossRef] [Green Version]
  11. Hortacsu, M. Heun functions and their uses in physics. arXiv 2013, arXiv:1101.0471v1. [Google Scholar]
  12. Zhang, Y.-Z. Exact polynomial solutions of second order differential equations and their applications. J. Phys. A Math. Theor. 2012, 45, 065206. [Google Scholar] [CrossRef]
  13. Saad, N.; Hall, R.L.; Trenton, V. Polynomial solutions for a class of second-order linear differential equations. Appl. Math. Comput. 2014, 226, 615–634. [Google Scholar] [CrossRef]
  14. Wen, F.-K.; Yang, Z.-Y.; Liu, C.; Yang, W.-L.; Zhang, Y.-Z. Exact polynomial solutions of Schrödinger equation with various hyperbolic potentials. Commun. Theor. Phys. 2014, 61, 153. [Google Scholar] [CrossRef]
  15. Sun, D.; You, Y.; Lu, F.; Chen, C.; Dong, S. The quantum characteristics of a class of complicated double ring-shaped non-central potential. Phys. Scr. 2014, 89, 045002. [Google Scholar] [CrossRef]
  16. Chiang, C.; Ho, C. Planar Dirac electron in Coulomb and magnetic fields: A Bethe ansatz approach. J. Math. Phys. 2002, 43, 43–51. [Google Scholar] [CrossRef] [Green Version]
  17. Chen, C.Y.; Lu, F.L.; Sun, D.S.; Dong, S.H. The origin and mathematical characteristics of the Super-Universal Associated-Legendre polynomials. Commun. Theor. Phys. 2014, 62, 331–337. [Google Scholar] [CrossRef]
  18. Chen, C.Y.; You, Y.; Lu, F.L.; Sun, D.S.; Dong, S.-H. Exact solutions to a class of differential equation and some new mathematical properties for the universal associated-Legendre polynomials. Appl. Math. Lett. 2015, 40, 90–96. [Google Scholar] [CrossRef]
  19. Chen, C.Y.; Lu, F.L.; Sun, D.S.; You, Y.; Dong, S.-H. Spin-orbit interaction for the double ring-shaped oscillator. Ann. Phys. 2016, 371, 183–198. [Google Scholar] [CrossRef]
  20. Dong, S.; Sun, G.; Falaye, B.J.; Dong, S.-H. Semi-exact solutions to position-dependent mass Schrödinger problem with a class of hyperbolic potential V0 tanh(ax). Eur. Phys. J. Plus 2016, 131, 176. [Google Scholar] [CrossRef]
  21. Ishkhanyan, A. Exact solution of the Schrödinger equation for the inverse square root potential. Europhys. Lett. 2015, 112, 10006. [Google Scholar] [CrossRef] [Green Version]
  22. Li, W.; Dai, W. Exact solution of inverse-square-root potential V(x) = −α/r. Ann. Phys. 2016, 373, 207–215. [Google Scholar] [CrossRef] [Green Version]
  23. Fernández, F.M. Comment on: Exact solution of the inverse-square-root potential V(r)=−α/r. Ann. Phys. 2017, 379, 83–85. [Google Scholar] [CrossRef]
  24. Ishkhanyan, A.M. Schrödinger potentials solvable in terms of the general Heun functions. Ann. Phys. 2018, 388, 456–471. [Google Scholar] [CrossRef] [Green Version]
  25. Ishkhanyan, A.M. Exact solution of the Schrödinger equation for a short-range exponential potential with inverse square root singularity. Eur. Phys. J. Plus 2018, 133, 83. [Google Scholar] [CrossRef] [Green Version]
  26. Maiz, F.; Alqahtani, M.M.; Al Sdran, N.; Ghnaim, I. Sextic and decatic anharmonic oscillator potentials: Polynomial solutions. Phys. B 2018, 530, 101–106. [Google Scholar] [CrossRef]
  27. Dong, S.; Dong, Q.; Sun, G.-H.; Femmam, S.; Dong, S.-H. Exact solutions of the Razavy Cosine Type potential. Adv. High Energy Phys. 2018, 2018, 5824271. [Google Scholar]
  28. Dong, Q.; Serrano, F.A.; Sun, G.-H.; Jing, J.; Dong, S.-H. Semi-exact Solutions of the Razavy Potential. Adv. High Energy Phys. 2018, 2018, 9105825. [Google Scholar] [CrossRef]
  29. Dong, Q.; Dong, S.S.; Hernández-Márquez, E.; Silva-Ortigoza, R.; Sun, G.H.; Dong, S.H. Semi-exact Solutions of Konwent Potential. Commun. Theor. Phys. 2019, 71, 231. [Google Scholar] [CrossRef] [Green Version]
  30. Dong, Q.; Sun, G.-H.; Jing, J.; Dong, S.-H. New findings for two new type sine hyperbolic potentials. Phys. Lett. A 2019, 383, 270–275. [Google Scholar] [CrossRef]
  31. Bazighifan, O. Some New Oscillation Results for Fourth-Order Neutral Differential Equations with a Canonical Operator. Math. Probl. Eng. 2020. [Google Scholar] [CrossRef]
  32. Agarwal, R.P.; Zhang, C.; Li, T. Some remarks on oscillation of second order neutral differential equations. Appl. Math. Comput. 2016, 274, 178–181. [Google Scholar]
  33. Euler, L. Institutiones Calculi Integralis, Vol. 2; Opera Omnia I: St. Petersburg, Russia, 1769; pp. 221–230. [Google Scholar]
  34. Scheffé, H. Linear differential equations with two-term recurrence formulas. J. Math. Phys. 1941, 21, 240–249. [Google Scholar] [CrossRef]
  35. Shore, S.D. On the second order differential equation which has orthogonal polynomial solutions. Bull. Calcutta Math. Soc. 1964, 56, 195–198. [Google Scholar]
  36. Hahn, W. On differential equations for orthogonal polynomials. Funkcialaj Ekuacioj 1978, 21, 1–9. [Google Scholar]
  37. Littlejohn, L.; Shore, S. Nonclassical orthogonal polynomials as solutions to second-order differential equations. Can. Math. Bull. 1982, 25, 291–295. [Google Scholar] [CrossRef]
  38. Littlejohn, L. On the classification of differential equations having orthogonal polynomial solutions. Ann. Mat. Pura Appl. 1984, 138, 35–53. [Google Scholar] [CrossRef]
  39. Al-Salam, W.A. Characterization Theorems for orthogonal polynomials. In Orthogonal Polynomials: Theory and Practice; Paul, N., Ed.; Springer: Dordrecht, The Netherlands, 1990. [Google Scholar]
  40. Atakishiyev, N.M.; Rahman, M.; Suslov, S.K. On classical orthogonal polynomials. Constr. Approx. 1995, 11, 181–226. [Google Scholar] [CrossRef]
  41. Ronveaux, A. Heun’s Differential Equation; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  42. Maier, R.S. The 192 solutions of the Heun equation. Math. Comput. 2007, 76, 811–843. [Google Scholar] [CrossRef]
  43. El-Jaick, L.J.; Bartolomeu, D.B.; Figueiredo, D.B. On certain solutions for confluent and double-confluent Heun equations. J. Math. Phys. 2008, 49, 082508. [Google Scholar] [CrossRef] [Green Version]
  44. Kwon, K.H.; Littlejohn, L.L. Sobolev orthogonal polynomials and second-order differential equations. Rocky Mt. J. Math. 1998, 28, 547–594. [Google Scholar] [CrossRef]
  45. Everitt, W.N.; Kwon, K.H.; Littlejohn, L.L.; Wellman, R. Orthogonal polynomial solutions of linear ordinary differential equations. J. Comput. Appl. Math. 2001, 133, 85–109. [Google Scholar] [CrossRef] [Green Version]
  46. Bluman, G.; Dridi, R. New solutions for ordinary differential equations. J. Symb. Comput. 2012, 47, 76–88. [Google Scholar] [CrossRef] [Green Version]
  47. Routh, E.J. On some properties of certain solutions of a differential equation of the second order. Proc. Lond. Math. Soc. 1884, 1, 245–262. [Google Scholar] [CrossRef] [Green Version]
  48. Bochner, S. Über Sturm-Liouvillesche Polynomsysteme. Math. Z. 1929, 29, 730–736. [Google Scholar] [CrossRef]
  49. Slavyanov, S.Y.; Lay, W. Special Functions, A Unified Theory Based on Singularities; Oxford Mathematical Monographs: Oxford, UK, 2000. [Google Scholar]
  50. Ince, E.L. Ordinary Differential Equations; Dover Publications: Mineola, NY, USA, 1956. [Google Scholar]
  51. Forsyth, A.R. A Treatise on Differential Equations, 6th ed.; Macmillan: London, UK, 1933. [Google Scholar]
  52. Ciftci, H.; Hall, R.L.; Saad, N. Asymptotic iteration method for eigenvalue problems. J. Phys. A Math. Gen. 2003, 36, 11807–11816. [Google Scholar] [CrossRef] [Green Version]
  53. Saad, N.; Hall, R.L.; Ciftci, H. Criterion for polynomial solutions to a class of linear differential equations of second order. J. Phys. A Math. Gen. 2005, 38, 1147. [Google Scholar] [CrossRef] [Green Version]
  54. Dong, S.-H. Wave Equations in Higher Dimensions; Springer: Dordrecht, The Netherlands, 2011. [Google Scholar]
  55. Rovder, J. Zeros of the polynomial solutions of the differential equation xy + (β0 + β1x + β2x2)y + (γnβ2x)y = 0. Mat. Căs. 1974, 24, 15. [Google Scholar]
Figure 1. The OPSOLDE Algorithm GUI. Upon inputting an integer value of n [ 2 , 9 ] , the coefficient field will resize appropriately and permit the user to choose the values for each coefficient. A working example is provided: by selecting “Dirac Coefficients” and the program will generate the results discussed in Section 5.2 of this article.
Figure 1. The OPSOLDE Algorithm GUI. Upon inputting an integer value of n [ 2 , 9 ] , the coefficient field will resize appropriately and permit the user to choose the values for each coefficient. A working example is provided: by selecting “Dirac Coefficients” and the program will generate the results discussed in Section 5.2 of this article.
Algorithms 13 00286 g001
Table 1. Tabulated values for parameters α i , β i , and τ i which would transform Equation (78) to the canonical form of the Heun equation.
Table 1. Tabulated values for parameters α i , β i , and τ i which would transform Equation (78) to the canonical form of the Heun equation.
i        3                2                1                0        
α i 1 ( c + 1 ) c0
β i δ + ϵ + γ c ( δ + γ ) + ϵ + γ γ c
τ i α β q
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bryenton, K.R.; Cameron, A.R.; Kirk, K.L.A.; Saad, N.; Strongman, P.; Volodin, N. On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application. Algorithms 2020, 13, 286. https://doi.org/10.3390/a13110286

AMA Style

Bryenton KR, Cameron AR, Kirk KLA, Saad N, Strongman P, Volodin N. On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application. Algorithms. 2020; 13(11):286. https://doi.org/10.3390/a13110286

Chicago/Turabian Style

Bryenton, Kyle R., Andrew R. Cameron, Keegan L. A. Kirk, Nasser Saad, Patrick Strongman, and Nikita Volodin. 2020. "On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application" Algorithms 13, no. 11: 286. https://doi.org/10.3390/a13110286

APA Style

Bryenton, K. R., Cameron, A. R., Kirk, K. L. A., Saad, N., Strongman, P., & Volodin, N. (2020). On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application. Algorithms, 13(11), 286. https://doi.org/10.3390/a13110286

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop