function[a0,ro,ENS,ENQ,EN,EWS,EWQ,EW]=neomezeny_1(lambda,mi);
% [a0,ro,ENS,ENQ,EN,EWS,EWQ,EW]=neomezeny_1(lambda,mi)
%
% Vypočítá prvek a0 stacionárního rozložení, intenzitu provozu  
% a charakteristiky systému hromadné obsluhy M|M|1|Inf|FIFO.
%
% Vstupní parametry:
% 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
% 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

ro=lambda/mi;
if ro<1
  a0=1-ro;
  ENS=ro;
  ENQ=lambda^2/(mi*(mi-lambda));
  EN=lambda/(mi-lambda);
  EW=1/(mi-lambda);
  EWS=1/mi;
  EWQ=lambda/(mi*(mi-lambda));
else
error('System se nemůže stabilizovat. Intenzita provozu je větší než 1.')
end;