function[a0,ro,PQ,ENS,ENQ,EN,EWS,EWQ,EW]=neomezeny_n(n,lambda,mi);
% [a0,ro,PQ,ENS,ENQ,EN,EWS,EWQ,EW]=neomezeny_n(n,lambda,mi)
%
% Vypočítá prvek a0 stacionárního rozložení, intenzitu provozu  
% a charakteristiky systému hromadné obsluhy M|M|n|Inf|FIFO.
%
% Vstupní parametry:
% n ......... počet linek obsluhy
% lambda .... parametr vstupního proudu
% mi ........ parametr obsluhy
%
% Výstupní parametry:
% a0 ........ pst, ľe v systému nebude ľádný zákazník
% ro ........ intenzita provozu (vyuľití systému)
% PQ ........ pst, že přicházející zákazník bude čekat ve frontě
% ENS ....... střední hodnota počtu obsluhovaných zákazníků
% ENQ ....... střední hodnota počtu zákazníků ve frontě
% EN ........ střední hodnota počtu zákazníků v systému
% EWS ....... střední hodnota doby, kterou zákazník stráví obsluhou
% EWQ ....... střední hodnota doby, kterou zákazník stráví ve frontě
% EW ........ střední hodnota doby, kterou zákazník stráví v systému

beta=lambda/mi;
ro=beta/n;
if ro<1
  a0=inv(sum(beta.^(0:n-1)./factorial(0:n-1))+n*beta^n/factorial(n)/(n-beta));
  PQ=a0*beta^n/(factorial(n)*(1-ro));
  ENS=n*ro;
  ENQ=PQ*ro/(1-ro);
  EN=ENQ+ENS;
  EWS=1/mi;
  EWQ=PQ*ro/(lambda*(1-ro));
  EW=EWS+EWQ;
else
error('System se nemuze stabilizovat. Intenzita provozu je vetsi nez 1.')
end;