KineticMonteCarlo

  private int kinmc(int endtime, double eps) {
            double m = 10;
            double c = 0.01;
            double time = 0;
            int state = 0, nstate = 0;
            double r01 = c * m;
            double r02 = c;
            double r1 = Math.Exp(-eps / 2);
            double r2 = Math.Exp(-eps);

            double[] rtot = { r01 + r02, r1, r2 };

            double[] swtime = new double[3];
            while(time <= endtime) {
                state = nstate;
                if(state == 0)
                    nstate = rng.NextDouble() * rtot[state] > r01 ? 2 : 1;
                else
                    nstate = 0;
                time += expverteilung(rtot[state]);
            }
            return state;
        }