Имитационное моделирование

#include <iostream>
#include <cstdlib>
#include <ctime>
#include <random>
 
int main() {
    const int N = 100;
    const int task_nb = 1000;
    default_random_engine eng;
    eng.seed(time(nullptr));
    uniform_int_distribution<unsigned> income(0, 1);
    auto gen_income = bind(income, eng);
    if (gen_income<0.4)
    {
    }
    else
        if (gen_income>0.4 && gen_income<0.7)
        {
        }
        else
        {
        }
 
    uniform_int_distribution<unsigned> income(1, 3), process(1, 9);
    auto gen_income = bind(income, eng);
    auto gen_process = bind(process, eng);
    double wait_avg_all = 0.0, process_avg_all = 0.0, prob_idle[2] = { 0.0 };
 
 
    for (int j = 0; j < N; ++j) {
        unsigned t_idle[2] = { 0, 0 }, t_wait = 0, t_income = 0, t_finish[2] = { 0, 0 }, t_process_all = 0;
        for (int i = 0; i < task_nb; ++i)
        {
            auto income = gen_income();
            auto t_process = gen_process();
            t_income += income;
            t_process_all += t_process;
            if (t_income >= t_finish[0])
 
            {
                // первый свободен - первый будет работать
                t_idle[0] += t_income - t_finish[0];
                t_finish[0] = t_income + t_process;
            }
            else if (t_income >= t_finish[1])
            {
                // первый занят, второй свободен - второй будет работать
                t_idle[1] += t_income - t_finish[1];
                t_finish[1] = t_income + t_process;
            }
            else
            {
                if (t_finish[0] <= t_finish[1])
                {
                    // первый освободится
                    t_wait += t_finish[0] - t_income;
                    t_finish[0] += t_process;
                }
                else
                {
                    // второй освободится
                    t_wait += t_finish[1] - t_income;
                    t_finish[1] += t_process;
                }
            }
        }
        wait_avg_all += double(t_wait) / double(task_nb);
        process_avg_all += double(t_process_all + t_wait) / task_nb;
        prob_idle[0] += double(t_idle[0]) / t_finish[0];
        prob_idle[1] += double(t_idle[1]) / t_finish[1];
    }
    std::cout << wait_avg_all / double(N) << " " << process_avg_all / double(N) << " "
        << prob_idle[0] / N << " " << prob_idle[1] / N << std::endl;
    getchar();
}