Comment by fsh

It's a Monte-Carlo simulation, not a statistical model. These are usually written in C++ for performance reasons.