I'm glad you liked this article Yuri! For the 1 million simulations as well as the rest of the simulations, I just continued with the data from the 1-by-3 matrix (the result of q*p^n) -- only difference is it was iterated for 1 million times.
Mechanical engineering student & writer at : www.ghannami.com