This should be slightly more robust than the legacy cFIREsim. We should be able to allow users to vary the amount of iterations. If we have a more math-oriented programmer out there, we could start adding in other types of distributions, rather than just the normal distribution.
I can provide the legacy code for this, to anyone willing to take up the task. (It isn't currently public).