Определить среднее время ожидания в справочном зале

#include <iostream>
#include <cmath>
 
using namespace std;
 
// Функция для вычисления факториала
double factorial(int n) {
    double result = 1.0;
    for (int i = 2; i <= n; ++i) {
        result *= i;
    }
    return result;
}
 
// Функция для вычисления P0
double calculateP0(double lambda, double mu, int c) {
    double sum = 0.0;
    for (int k = 0; k < c; ++k) {
        sum += pow(lambda / mu, k) / factorial(k);
    }
    double lastTerm = pow(lambda / mu, c) / (factorial(c) * (1 - (lambda / (c * mu))));
    return 1.0 / (sum + lastTerm);
}
 
int main() {
    // Ввод исходных данных
    int M;          // Количество терминалов
    double P;       // Среднее число пользователей в сутки
    double t;       // Среднее время обслуживания в минутах
    double T0;      // Время работы за сутки в часах
 
    // Примерные значения
    M = 5;
    P = 140;      // чел/сутки
    t = 40;       // минуты
    T0 = 24;      // часа
 
    // Расчет интенсивности входящего потока (чел/мин)
    double lambda = P / (T0 * 60.0); // чел/мин
 
    // Среднее время обслуживания (обратная)
    double mu = 1.0 / t; // обслуживание за минуту
 
    // Загрузка системы
    double rho = lambda / (M * mu);
    if (rho >= 1.0) {
        cout << "Система перегружена (rho >= 1), расчет невозможен." << endl;
        return 1;
    }
 
    // Вычисление P0
    double P0 = calculateP0(lambda, mu, M);
 
    // Средняя длина очереди Lq
    double numerator = pow(lambda / mu, M) * rho;
    double denominator = factorial(M) * pow(1 - rho, 2);
    double Lq = P0 * numerator / denominator;
 
    // Среднее время ожидания в очереди
    double Wq = Lq / lambda; // в минутах
 
    cout << "Среднее время ожидания в очереди: " << Wq << " минут" << endl;
 
    return 0;
}