The simulations were carried out in Scilab. The code for the sequence of approximate solutions can be download here .
The code generates the following plot: