Here is the numerical solution

Initialize some variables.
We need to run through the integration once to get a first values for the error.  
Here is the main iteration loop.
Here are some plots of the results


Converted by Mathematica      April 30, 2000