Example 2.  Use Frobenius series to solve the D. E.   
        [Graphics:Images/FrobeniusSeriesMod_gr_145.gif].  
Solution 2.

Determine the nature of the singularity at [Graphics:../Images/FrobeniusSeriesMod_gr_146.gif].  

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


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

Construct the Indicial Equation.

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

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

 

Find the Roots of the Indicial Equation.

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

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

 

Form the first Frobenius solution corresponding to the root  [Graphics:../Images/FrobeniusSeriesMod_gr_153.gif].

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


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

Form the set of equations to solve and do it.

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


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

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

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

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

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

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

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

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

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

The first Frobenius solution is:

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


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

Form the second Frobenius solution corresponding to the root  [Graphics:../Images/FrobeniusSeriesMod_gr_168.gif].

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


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

Form the set of equations to solve and do it.

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


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

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

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

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

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

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

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

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

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

The second Frobenius solution is:

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


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

After you are done, use Mathematica's DSolve subroutine to get the answer and check out its series expansion.
This will require slight fussing around with the appropriate multiple of  "Sinh".

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


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

At this time we could plot the series approximations and the analytic solutions. To see the difference in the graphs we will reduce the number of terms in the series.

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

[Graphics:../Images/FrobeniusSeriesMod_gr_186.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_187.gif]


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

[Graphics:../Images/FrobeniusSeriesMod_gr_189.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_190.gif]

The recursive formula for the coefficients.

    If we look at the series in more depth we will be able to obtain the analytic solutions as infinite sums.  First find the recursive formula for the coefficients of  [Graphics:../Images/FrobeniusSeriesMod_gr_191.gif].  If you try this be sure to use the  " := "  replacement delayed structure to avoid an infinite recursion.  For this example the trial term  [Graphics:../Images/FrobeniusSeriesMod_gr_192.gif] works with the replacements [Graphics:../Images/FrobeniusSeriesMod_gr_193.gif].  

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


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

If you can't get the above computation to work, then just type in the recursive formula.

 

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

Now look at each series individually.  The first Frobenius corresponds to  [Graphics:../Images/FrobeniusSeriesMod_gr_197.gif].  

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


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

Now look at the second Frobenius series corresponding to  [Graphics:../Images/FrobeniusSeriesMod_gr_200.gif].  

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


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

The explicit formulas for the coefficients for the first Frobenius series are:

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


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

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

[Graphics:../Images/FrobeniusSeriesMod_gr_206.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_207.gif]

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


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

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

[Graphics:../Images/FrobeniusSeriesMod_gr_211.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_212.gif]

Note. We cannot add up the infinite number of terms in this sequence because Mathematica does not know this formula yet.

Remark. However, Mathematica can treat infinite series "like a function."

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

[Graphics:../Images/FrobeniusSeriesMod_gr_214.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_215.gif]

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

[Graphics:../Images/FrobeniusSeriesMod_gr_217.gif]
[Graphics:../Images/FrobeniusSeriesMod_gr_218.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004