function[C]=opt_neomezeny_n_alt(n,lambda,mi,c1,c2);
% [C]=opt_neomezeny_n_alt(n,lambda,mi,c1,c2)
%
% 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
% c1 ........ naklady na pobyt zakaznika v SHO za jednotku casu
% c2 ........ naklady na provoz linky 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=c1*(an*ro/(1-ro)^2+n*ro)+c2*n;
else
error('System se nemuze stabilizovat. Intenzita provozu je vetsi nez 1.')
end;