function[C]=opt_neomezeny_n(n,lambda,mi,k1,k2);
% [C]=opt_neomezeny_n(n,lambda,mi,k1,k2)
%
% Vypočítá hodnotu kriteriální funkce v systému hromadné
% obsluhy M|M|n|Inf|FIFO.
%
% Vstupní parametry:
% n ......... počet linek obsluhy
% lambda .... parametr vstupního proudu
% mi ........ parametr obsluhy
% k1 ........ náklady na čekajícího zákazníka
% k2 ........ náklady na nevyužitou linku obsluhy
%
% Výstupní parametry:
% C ......... hodnota kriteriální funkce pro n linek obsluhy
%             (celkové průměrné náklady na provoz 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));
  an=beta^n*a0/factorial(n);
  C=k1*an*ro/(1-ro)^2+k2*(n-n*ro);
else
error('System se nemůže stabilizovat. Intenzita provozu je větší než 1.')
end;