function[ro, a, EN, EW]=neomezeny_1_netrpelivy(lambda,mi,b)
% [ro, a, EN, EW]=neomezeny_1_netrpelivy(lambda,mi,b)
% M/M/1/Inf/FIFO s netrpelivymi zakazniky
% vstupni parametry:
  % lambda ... intenzita vstupu
  % mi ... intenzita obsluhy
  % b ... radkovy vektor pravdepodobnosti chovani zak. v zavislosti na obsazenosti
% vystupni parametry:
  % ro ... intenzita provozu
  % a ... stacionarni vektor
  % EN ... stredni hodnota poctu zakazniku
  % EW ... stredni hodnota celkove doby

ro = lambda/mi;
if ro > 1 
    error('System se nemuze stabilizovat. Intenzita provozu je vetsi než 1.')
end
c = [1 cumprod(b)];
j = 0:(length(c)-1);
a = (diag(c' * ro.^j) * inv(ro.^j*c'));
EN = j*a;
EW = EN/lambda;