Моделирование

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <time.h>
using namespace std;
 
int main()
{
    int number_of_experiments, countY, countN, summ = 0, summN = 0; // количество опытов, количество обслуженных заявок, количество потеряных заявок
    double  time_to_work, arrival_intensity, tau;// время на работу, интенсивность поступления, тау(время обслуживания)
    double release_time_line, time_of_admission, orbit_time, orbit_for_while; // время освобождения линии, время поступления, орбитное время, орбитное для цикла while
    cout << "Plese enter the number of experiments: " << endl;
    cin >> number_of_experiments;
    cout << "Please enter time to work: " << endl;
    cin >> time_to_work;
    cout << "Please enter arrival intensity :" << endl;
    cin >> arrival_intensity;
    cout << "Please enter tau(service time): " << endl;
    cin >> tau;
    srand(time(NULL));
    for (int j = 1; j <= number_of_experiments; j++){
        release_time_line = 0; // время освобождения линии
        time_of_admission = 0; // время поступления
        orbit_time = 0; // орбитное время
        orbit_for_while = 0;
        countY = 0; countN = 0; // количество обслуженных заявок, количество потеряных заявок
        while(time_of_admission < time_to_work) {
            if(orbit_time == 0) {
                if(release_time_line <= time_of_admission) {
                    release_time_line = time_of_admission + tau;
                    countY = countY + 1;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                } else {
                    orbit_time = time_of_admission + 2*tau;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                }
            } else {
                if(release_time_line <= time_of_admission) {
                    if((orbit_time - time_of_admission) >= tau) {
                        release_time_line = time_of_admission + tau;
                        countY = countY + 1;
                        time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                    } else {
                        release_time_line = orbit_time + tau;
                        countY = countY + 1;
                        orbit_for_while = orbit_time;
                        orbit_time = 0;
                        while(orbit_for_while > time_of_admission) {
                            countN = countN + 1;
                            time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                        }
                    }
                } else {
                    countN = countN + 1;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                }
            }
        }
        summN = summN + countN;
        summ = summ + countY + countN;
    }
    cout << "result program: " << double(summN)/summ << endl;
    cout << "result formula: " << (2*tau*arrival_intensity*exp(arrival_intensity*tau)+exp(arrival_intensity*tau)+0.5*exp(arrival_intensity*tau)-0.5*exp(-(arrival_intensity*tau))-1-2*arrival_intensity*tau)/(3*exp(arrival_intensity*tau)+2*arrival_intensity*tau*exp(arrival_intensity*tau)-2-2*arrival_intensity*tau) << endl;
    system("pause");
    return 0;
}