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

        [Graphics:Images/PicardIterationMod_gr_211.gif].  

Solution 10.

First define the function  [Graphics:../Images/PicardIterationMod_gr_212.gif], and  [Graphics:../Images/PicardIterationMod_gr_213.gif] and the initial conditions  [Graphics:../Images/PicardIterationMod_gr_214.gif]    [Graphics:../Images/PicardIterationMod_gr_215.gif]  by typing:

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

[Graphics:../Images/PicardIterationMod_gr_217.gif]
[Graphics:../Images/PicardIterationMod_gr_218.gif]
[Graphics:../Images/PicardIterationMod_gr_219.gif]

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

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

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

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

 

We are done.

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

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

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


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

Observe that the first six terms of the Picard iteration agrees with the first six terms of the Maclaurin series.  

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

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005