Моделирование многоканальных СМО методом последовательной проводки заявок

#include "stdafx.h" 
#include "conio.h"
#include "locale.h"
#include "stdlib.h"
#include "time.h"
#include "math.h"
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(0,"RUSSIAN");
    srand(time(NULL));
 
    const int N=12;
    int n, i,j,k, m=2, nk; 
    float Mn=5,sigma=0.5, Mn1=10,sigma1=0.5, a,b,a1,b1, fi, min, Sum_og, Sum_pr, x[N],
        Tsv[2]; //время освобождения канала
 
    printf("Количество заявок = ");
    scanf("%d", &n);
    float *t=new float  [n];    //момент прихода очередной заявки
    float *T_n=new float[n];    //момента начала обслуживания заявки
    float *tau=new float[n];    //время через которое приходит заявка
    float *T_obs=new float[n];  //время обслуживания заявки
    float *T_pr=new float[n];    //время простоя первого канала
    //float *T_pr2=new float[n];    //время простоя второго канала
    float *T_og=new float [n];      //время ожидания заявки
    
    //Высчитываем границы (время через которое прийдет след. заявка)
    a=(Mn/N - sqrtf(3)*sigma/N);
    b=(Mn/N + sqrtf(3)*sigma/N);
    printf("Границы прихода след.заявки: [%.2f, %.2f]\n", a, b);
 
    //Высчитываем границы (время сколько она будет обслуживаться)
    a1=(Mn1/N - sqrtf(3)*sigma1/N);
    b1=(Mn1/N + sqrtf(3)*sigma1/N);
    printf("Границы обслуживания: [%.2f, %.2f]\n\n", a1, b1);
    
    Tsv[0]=0; Tsv[1]=0; min=10000; Sum_og=0, Sum_pr=0;
    for(i=0; i<n; i++){ //цикл по заявкам
        fi=0+(float)rand()*(1-0)/RAND_MAX;
        tau[i]=a+(b-a)*fi;//время через которое приходит заявка
        if(i==0) t[i]=0;    //момент прихода первой заявки
        else t[i]=t[i-1]+tau[i];
        
        for(j=0; j<m; j++){ //цикл по каналам
            if(Tsv[j]<min){ //определяем какой из каналов освободился раньше
                min=Tsv[j];
                nk=j;
            }
        }
        printf("%d\n", nk);
        
        //проверяем имеются ли свободные каналы
        if(t[i]<min){//все каналы заняты (момент начала обслуживания равен минимальному моменту освобождения)
            T_pr[i]=0;
            T_og[i]=t[i]-Tsv[i]; //время ожидания заявки
            Sum_og+=T_og[i];
            T_n[i]=min;
        } 
        else { //если есть свободные каналы
            T_og[i]=0;
            T_pr[i]=t[i]-Tsv[i];
            Sum_pr+=T_pr[i];
            T_n[i]=t[i];//(момент начала равен моменту прихода заявки)
        }   
        T_obs[i]=0;
        for(k=0; k<N; k++){
            fi=0+(float)rand()*(1-0)/RAND_MAX;
            x[k]=(a1+(b1-a1)*fi);
            T_obs[i]+=x[k];  //нормальный закон распределения
            
        }
        Tsv[nk]=T_n[i]+T_obs[i];
        //printf("%f  %f\n", Tsv[0], Tsv[1]);
    }
    
    delete []t, delete []tau, delete []T_obs, delete []T_n, delete []T_og, delete []T_pr; //delete []T_pr2;
    getch();
    return 0;
}