Picard's Method for Numerical Solution to First Order Differential Equations

 

 

Problem Statement:    Find an approximate solution to the first order, ordinary d.e.

 

                    dy/dx   =  f [ (x,y(x) ]     subject to the condition        y(a)  =  b

 

 

 

Strategy:    Use successive approximations (another name for Picard's Method)

                                                                                                c             c

First Integrate    dy   =    f[x,y(x)] dx   from a to c    or         dy   =    ∫ f [x,y(x)] dx    

                                                                                               a              a

where a is the initial value of x and c is an arbitrary value.

 

 

                                                c

Then               y(c)  ˗  y(a)  =  ∫ f [x,y(x)] dx         but  y(a)  =  b      so

                                                a

 

                                                c

                       y(c)  ˗  b  =       ∫ f [x,y(x)] dx        here  x  is a dummy variable

                                                a                          (switch variable of integration to  t)

 

 

                                                       x

for arbitrary x,      y(x)    =    b  +  ∫ f [ t,y(t) ] dt   

                                                      a

 

 

                                                                                     x

Successive steps are given by:     yn+1(x)    =    b  +    ∫ f [ t,yn(t) ] dt    start with yo = b

                                                                                    a

 

 

Click here for an example.

Click here to return to the discussion on approximate methods.

 

 

 




Copyright © 2017 Richard C. Coddington

All rights reserved.