Example 9.  Use Picard iteration to find the solution of the second order I.V.P.  

        [Graphics:Images/PicardIterationMod_gr_194.gif].  

Solution 9.

First define the function  [Graphics:../Images/PicardIterationMod_gr_195.gif], and  [Graphics:../Images/PicardIterationMod_gr_196.gif] and the initial conditions  [Graphics:../Images/PicardIterationMod_gr_197.gif]    [Graphics:../Images/PicardIterationMod_gr_198.gif]  by typing:

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

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

 

Picard iteration for generating the first six approximations is started with the Mathematica command:

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

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

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

 

We are done.

Techniques from calculus can be used to find the solution  [Graphics:../Images/PicardIterationMod_gr_204.gif].

We can express  [Graphics:../Images/PicardIterationMod_gr_205.gif]  as a Maclaurin series and observe that the sequence  [Graphics:../Images/PicardIterationMod_gr_206.gif]  is converging to the solution

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


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

For this example all the terms in the Picard approximation agree with all of the terms in the Maclaurin series.  

Aside.  We can let Mathematica solve the differential equation.  This is just for fun.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005