Method of Frobenius.

    This method is attributed to the german mathematician Ferdinand Georg Frobenius (1849-1917 ).  Assume that  [Graphics:Images/FrobeniusSeriesProof_gr_35.gif]  is regular singular point of the differential equation

        [Graphics:Images/FrobeniusSeriesProof_gr_36.gif].  
    
A Frobenius series (generalized Laurent series) of the form  

[Graphics:Images/FrobeniusSeriesProof_gr_37.gif]

where  [Graphics:Images/FrobeniusSeriesProof_gr_38.gif]  can be used to solve the differential equation.  The parameter [Graphics:Images/FrobeniusSeriesProof_gr_39.gif] must be chosen so that when the series is substituted into the D.E. the coefficient of the smallest power of  [Graphics:Images/FrobeniusSeriesProof_gr_40.gif] is zero.  This is called the indicial equation.  Next, a recursive equation for the coefficients is obtained by setting the coefficient of  [Graphics:Images/FrobeniusSeriesProof_gr_41.gif]  equal to zero. Caveat. There are some instances when only one Frobenius solution can be constructed.

Definition (Indicial Equation).  The parameter [Graphics:Images/FrobeniusSeriesProof_gr_42.gif] in the Frobenius series is a root of the indicial equation  

        [Graphics:Images/FrobeniusSeriesProof_gr_43.gif].

Assuming that the singular point is  [Graphics:Images/FrobeniusSeriesProof_gr_44.gif],  we can calculate  [Graphics:Images/FrobeniusSeriesProof_gr_45.gif]  as follows:

        [Graphics:Images/FrobeniusSeriesProof_gr_46.gif]
and
        [Graphics:Images/FrobeniusSeriesProof_gr_47.gif]

Definition of  [Graphics:Images/FrobeniusSeriesProof_gr_48.gif]  We state the following definition of  [Graphics:Images/FrobeniusSeriesProof_gr_49.gif]

        [Graphics:Images/FrobeniusSeriesProof_gr_50.gif].  

The exponents of the singularity are the roots  [Graphics:Images/FrobeniusSeriesProof_gr_51.gif]  of  [Graphics:Images/FrobeniusSeriesProof_gr_52.gif].   

Derivation.

 

The Recursive Formula for  [Graphics:Images/FrobeniusSeriesProof_gr_78.gif]  

    We are now in a position to derive the recursive formula for the sequence of coefficients  [Graphics:Images/FrobeniusSeriesProof_gr_79.gif]  for the Frobenius series solution  

[Graphics:Images/FrobeniusSeriesProof_gr_80.gif]

The recursive formula for computing  [Graphics:Images/FrobeniusSeriesProof_gr_81.gif]  is

[Graphics:Images/FrobeniusSeriesProof_gr_82.gif]

where

[Graphics:Images/FrobeniusSeriesProof_gr_83.gif]

Derivation.

 

Starting with the differential equation  

[Graphics:../Images/FrobeniusSeriesProof_gr_84.gif]

Rewrite it in the form  

[Graphics:../Images/FrobeniusSeriesProof_gr_85.gif]

Multiply each term by the factor  [Graphics:../Images/FrobeniusSeriesProof_gr_86.gif].   

[Graphics:../Images/FrobeniusSeriesProof_gr_87.gif]


Use series for all the terms  [Graphics:../Images/FrobeniusSeriesProof_gr_88.gif].  

[Graphics:../Images/FrobeniusSeriesProof_gr_89.gif]

Make the substitutions  [Graphics:../Images/FrobeniusSeriesProof_gr_90.gif],  [Graphics:../Images/FrobeniusSeriesProof_gr_91.gif] and regroup the second and third terms in the D. E. as follows.  

 

[Graphics:../Images/FrobeniusSeriesProof_gr_92.gif]

Making all the series substitutions we get  

[Graphics:../Images/FrobeniusSeriesProof_gr_93.gif]

 

Move the terms [Graphics:../Images/FrobeniusSeriesProof_gr_94.gif] and [Graphics:../Images/FrobeniusSeriesProof_gr_95.gif]  into the summations where they belong  

[Graphics:../Images/FrobeniusSeriesProof_gr_96.gif]

 

Recall the following Cauchy product form for multiplying two infinite series  

[Graphics:../Images/FrobeniusSeriesProof_gr_97.gif]
[Graphics:../Images/FrobeniusSeriesProof_gr_98.gif]

Multiply each side by of the above Cauchy products by   [Graphics:../Images/FrobeniusSeriesProof_gr_99.gif].  

 

[Graphics:../Images/FrobeniusSeriesProof_gr_100.gif]
[Graphics:../Images/FrobeniusSeriesProof_gr_101.gif]

Now substitute the above series into

[Graphics:../Images/FrobeniusSeriesProof_gr_102.gif]

Get

[Graphics:../Images/FrobeniusSeriesProof_gr_103.gif]

 

Which can be combined as follows

 

[Graphics:../Images/FrobeniusSeriesProof_gr_104.gif]

And simplified

 

[Graphics:../Images/FrobeniusSeriesProof_gr_105.gif]

And further simplified

 

[Graphics:../Images/FrobeniusSeriesProof_gr_106.gif]

From which we get the equations

 

[Graphics:../Images/FrobeniusSeriesProof_gr_107.gif]

for  n = 0,1,2,...

 

Now split off the last term in the sum corresponding to  k=n  and write this as follows

 

[Graphics:../Images/FrobeniusSeriesProof_gr_108.gif]

From which we easily get

 

[Graphics:../Images/FrobeniusSeriesProof_gr_109.gif]

Recall that  [Graphics:../Images/FrobeniusSeriesProof_gr_110.gif]  and   

        [Graphics:../Images/FrobeniusSeriesProof_gr_111.gif]  
        
so the above quantity can be written in the form  

 

[Graphics:../Images/FrobeniusSeriesProof_gr_112.gif]

From this we can easily obtain the desired recursion formula for computing [Graphics:../Images/FrobeniusSeriesProof_gr_113.gif].

 

[Graphics:../Images/FrobeniusSeriesProof_gr_114.gif]

for  n = 1,2,3,...

 

Caveat.  Although the general formula above can easily be used in a recursive program it is slow because of the recursion depth that occurs with each successive term in the sequence.  For many practical problems the series for  [Graphics:../Images/FrobeniusSeriesProof_gr_115.gif]  and   [Graphics:../Images/FrobeniusSeriesProof_gr_116.gif] are short and the traditional methods of pencil and paper will create a simple recursion formula.

 

Caveat.  Assume that the roots of the indicial equation are [Graphics:../Images/FrobeniusSeriesProof_gr_117.gif] and [Graphics:../Images/FrobeniusSeriesProof_gr_118.gif].  If  [Graphics:../Images/FrobeniusSeriesProof_gr_119.gif]  or if  the roots differ by an integer, then only one Frobenius series solution can be found and other methods such as reduction of order must be used to create the second solution.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004