Final version of the π approximation.
I performed [a benchmark](https://gist.github.com/kdungs/10963341) and found out that the loop-based version is
slightly faster (clang-503.0.40) so we go with that one. It is also more
memory efficient.