Система массового обслуживания

program smo;
var j,m,x, T: integer;
     p, prost, otk, obsl, a, k, y1, y2, y3, y4: real;
     factorial: longint;
     n, i: byte;
begin
writeln ('Введите скорость потока заявок j=');
readln (j);
writeln ('Введите скорость обслуживания одним каналом m=');
readln (m);
writeln ('Введите количество обслуживащих каналов х=');
readln (x);
writeln ('Введите тип СМО T: 0 - с откатами,1 - с неограниченным временем ожидания');
readln (T);
if T=0 then
begin
p:=j/m;
    begin
    write('n = '); readln(n);
    factorial := 1;
    for i:=2 to n do
    factorial := factorial * i;
    readln;
    end;
prost:=1/(p/factorial);
otk:= p*prost/x;
obsl:=1-otk;
A:=j*obsl;
k:=p*obsl/x;
writeln ('Вероятность простоя prost=', prost:2:2);
writeln ('Вероятность отказа в обслуживании otk=', otk:2:2);
writeln ('Вероятность обслуживания obsl=', obsl:2:2);
writeln ('Абсолютная пропускная способность A=', A:2:2);
writeln ('Коэффициент загрузки СМО k=', k:2:2);
readln
end
else
begin
    write('n = '); readln(n);
    factorial := 1;
    for i:=2 to n do
    factorial := factorial * i;
    readln;
    end;
prost:=1/((p/factorial)+p/(x*(x-p)));
y1:= p*prost/((x-1)*(x-p)*(x-p));
y2:=y1+p;
y3:=y1/j;
y4:=y2/j;
k:=p/x;
writeln ('Вероятность простоя prost=', prost:2:2);
writeln ('Среднее число заявок в очереди у1=', y1:2:2);
writeln ('Среднее число заявок в системе y2=', y2:2:2);
writeln ('Среднее время ожидания заявки в очереди у3=', y3:2:2);
writeln ('Среднее время нахождения заявки в системе у4=', y4:2:2);
writeln ('Коэффициент загрузки СМО k=', k:2:2);
readln;
end.