Implemented functions for rate constants and jump rates
First implementation of cluster simulation. It seems to work although output is somewhat unexpected. Also model is way too simple