function[C]=opt_neomezeny_n(n,lambda,mi,k1,k2);
% [C]=opt_neomezeny_n(n,lambda,mi,k1,k2)
%
% Vypocita hodnotu kriterialni funkce v systemu hromadne
% obsluhy M|M|n|Inf|FIFO.
%
% Vstupni parametry:
% n ......... pocet linek obsluhy
% lambda .... parametr vstupniho proudu
% mi ........ parametr obsluhy
% k1 ........ naklady na cekajiciho zakaznika za jednotku casu
% k2 ........ naklady na nevyuzitou linku obsluhy za jednotku casu
%
% Vystupni parametry:
% C ......... hodnota kriterialni funkce pro n linek obsluhy
%             (celkove prumerne naklady na provoz systemu)

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 nemuze stabilizovat. Intenzita provozu je vetsi nez 1.')
end;