Advanced computing with IPython

Oh, and also the Python version returns the value from n=1e8-1, while the Fortran version goes up to n=1e8, because 'range' takes an exclusive upper bound and 'do' is inclusive. But even after fixing that, it seems there's still a (slightly smaller) difference in their outputs.