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