There are some errors to fix: when run for (I believe) 1416 iterations,
the code encounters a convergence issue when trying to run the next
iteration. Besides that, everything is done, I think.