Реализация нескольких генераторов случайных чисел с разными законами распределения в одном проекте

#include <iostream>
#include <random>
#include <fstream>
 
using namespace std;
 
int main()
{
    int i;
    double mean;    double disp;    int n;//берем из файла
    const int N = 10000;    double m[N]; double r[N];//массив для СЧ
    const int r_kr = 40000; //критическое значение для условия
    
    setlocale(LC_ALL, "rus");   //для кириллицы
    ifstream fin("data_in.txt", ios_base::in); //файл входных данных, режим чтения
    ofstream fout("data_out.txt", ios_base::out, ios_base::trunc); //файл выходных данных, режим записи и очистки
 
    fin >> mean;    fin >> disp; fin >> n;
 
    random_device rd; //задание случайного числа (зерна) запускающего генератор
    mt19937 gen(rd()); //Вихрь Мерсенна
    normal_distribution <double> distribution(mean, disp);//Нормальное распределение
 
    double min = DBL_MAX, max = DBL_MIN;//Анализ данных
    for (int i = 0; i < N; i++)
    {       m[i] = distribution(gen);
                if (m[i] < min) min = m[i];
                    if (m[i] > max) max = m[i];     }
 
    for (int i = 0; i < n; i++)//Выводим данные по ряду в файл
    {       fout << m[i] << "\n";   }
    fin.close();    fout.close();//закрываем входной и выходной файлы
 
    for (int i = 0; i < n; i++)//обозначаем новый массив
    {       r[i] = 0;   }
        
    for (int i = 0; i < n; i++)//работаем с новым массивом 
    {       r[i] = m[i] * 10;   }
 
    int k; k = 0; i = 0;//переменная для счета количества операций удовлетворяющих условию 
    do
    {       if (r[i] < r_kr)
                k++;
                    i++;    }
                        while (i < n);
            cout << k << "\n";//выводим количество итераций удовлетворяющих условию
 
    for (int i = 0; i < n; i++)//цикл для вывода значений ряда удовлетворяющих условию
        {   if (r[i] < r_kr)   
                    cout << r[i] << "\n";       }   
        
    cout << k << "/" << n << "=" << k / static_cast <double> (n) << "\n";//делим количество итераций удовлетворяющих условию на 
                                                                        //общее число значений (задается пользователем в файле входных данных) 
}

#include <iostream>
#include <random>
#include <fstream>
 
using namespace std;
 
int main()
{
    double mean_ro; double disp_ro; int n; double mean_lambda; double disp_lambda; double a; double b;//берем из файла
    const int N = 10000;
    double gauss_ro[N]; double gauss_lambda[N]; double weibull_c[N]; //double geom[g]; //переменные для рядов СЧ
    //double r[N]; //массив для результатов
 
    const int r_kr = 40000; //критическое значение для условия
 
    setlocale(LC_ALL, "rus");   //для кириллицы
    ifstream fin("data_in.txt", ios_base::in); //файл входных данных, режим чтения
    ofstream fout("data_out.txt", ios_base::out, ios_base::trunc); //файл выходных данных, режим записи и очистки
 
    fin >> n; fin >> mean_ro;   fin >> disp_ro;  fin >> mean_lambda; fin >> disp_lambda; fin >> a; fin >> b;
 
    random_device rd; //задание случайного числа (зерна) запускающего генератор
    mt19937 gen(rd()); //Вихрь Мерсенна
 
    normal_distribution <double> distribution_ro(mean_ro, disp_ro);//Нормальное распределение плотность и удельная теплоемкость?
        double min = DBL_MAX, max = DBL_MIN;//Анализ данных Гаусс Плотность
    for (int i = 0; i < N; i++)
    {
        gauss_ro[i] = distribution_ro(gen);
        if (gauss_ro[i] < min) min = gauss_ro[i];
        if (gauss_ro[i] > max) max = gauss_ro[i];
    }
 
    normal_distribution <double> distribution_lambda(mean_lambda, disp_lambda);
            for (int i = 0; i < N; i++)
    {
        gauss_lambda[i] = distribution_lambda(gen);
        if (gauss_lambda[i] < min) min = gauss_lambda[i];
        if (gauss_lambda[i] > max) max = gauss_lambda[i];
    }
 
   weibull_distribution <double> distribution_c(a, b);//Распределение Вейбулла: ФОРМА - МАСШТАБ
        for (int i = 0; i < N; i++)
    {
        weibull_c[i] = distribution_c(gen);
        if (weibull_c[i] < min) min = weibull_c[i];
        if (weibull_c[i] > max) max = weibull_c[i];
    }
 
    for (int i = 0; i < n; i++)//Выводим данные по ряду в файл
    {
        fout << gauss_ro[i] << "\t" << gauss_lambda[i] << "\t"<< weibull_c[i] << "\n";
    }
 
    fin.close();    fout.close();//закрываем входной и выходной файлы
}