Picard's Method for Higher Order D.E.'s

   
    It is easy to extend the idea of Picard iteration to higher order problems.  For illustration purposes we mention second order D.E.'s.
    
Theorem 3 (Picard Iteration for Second Order D. E.'s).  Given the second order initial value problem  
        [Graphics:Images/PicardIterationProof_gr_46.gif],   with  
(10)    
        [Graphics:Images/PicardIterationProof_gr_47.gif].  
The solution to the I.V.P in (10) is found by constructing recursively a sequence [Graphics:Images/PicardIterationProof_gr_48.gif] of functions
        [Graphics:Images/PicardIterationProof_gr_49.gif],   and  
(11)
        [Graphics:Images/PicardIterationProof_gr_50.gif]    for   [Graphics:Images/PicardIterationProof_gr_51.gif].  
Then the solution [Graphics:Images/PicardIterationProof_gr_52.gif] to (10) is given by the limit:

(12)        [Graphics:Images/PicardIterationProof_gr_53.gif].
Proof.

Begin by reformulating (10) as an equivalent integral equation.  If  [Graphics:../Images/PicardIterationProof_gr_54.gif] is the solution of (10), then  

(13)        [Graphics:../Images/PicardIterationProof_gr_55.gif]  

The Fundamental Theorem of Calculus is used to integrate the left side of (13), and the result after rearranging terms is;

(14)        [Graphics:../Images/PicardIterationProof_gr_56.gif]

Now integrate (14) one more time and get

        [Graphics:../Images/PicardIterationProof_gr_57.gif]
    then
        [Graphics:../Images/PicardIterationProof_gr_58.gif]

which can be further simplified to produce equivalent integral equation

(15)        [Graphics:../Images/PicardIterationProof_gr_59.gif].  

    If we use equation (15) and input a formula for  [Graphics:../Images/PicardIterationProof_gr_60.gif] and its derivative  [Graphics:../Images/PicardIterationProof_gr_61.gif]  in the integrand  [Graphics:../Images/PicardIterationProof_gr_62.gif],  then the function  [Graphics:../Images/PicardIterationProof_gr_63.gif]  on the left side is considered output.  
    Start the iteration with the initial function

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

then define the next function  [Graphics:../Images/PicardIterationProof_gr_65.gif]  as follows

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

Next  [Graphics:../Images/PicardIterationProof_gr_67.gif]  is used to construct  [Graphics:../Images/PicardIterationProof_gr_68.gif]  as follows  

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

The process is repeated, and once  [Graphics:../Images/PicardIterationProof_gr_70.gif]  has been obtained, the next function is given recursively by  

(16)        [Graphics:../Images/PicardIterationProof_gr_71.gif].

    In a fashion similar to the first order case we must take the limit as  [Graphics:../Images/PicardIterationProof_gr_72.gif]  in (16).  Assume that the limit (12) exists, then [Graphics:../Images/PicardIterationProof_gr_73.gif] and [Graphics:../Images/PicardIterationProof_gr_74.gif] and we write  

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

If there is "no problem" when taking limits on the right side then we might expect the following

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

This is the "intuitive proof"  of equation (15).                                  Q.E.D.

For More Proof.  
    
    A more rigorous proof based on of uniform convergence can be found in the article:

James Fabrey, Picard's Theorem (in Classroom Notes),The American Mathematical Monthly, Vol. 79, No. 9. (Nov., 1972), pp. 1020-1023, Jstor.  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005