pfqn_xzgsbup.m

%{
%{
 % @file pfqn_xzgsbup.m
 % @brief Upper asymptotic bound on throughput (Zahorjan-Gittelsohn-Schweitzer-Bryant).
%}
%}
 
%{
%{
 % @brief Upper asymptotic bound on throughput (Zahorjan-Gittelsohn-Schweitzer-Bryant).
 % @fn pfqn_xzgsbup(L, N, Z)
 % @param L Service demand vector.
 % @param N Population.
 % @param Z Think time.
 % @return X Upper bound on throughput.
%}
%}
function X=pfqn_xzgsbup(L,N,Z)
M=length(L);
R=Z+sum(L)+max(L)*(N-1);
for ist=1:M
    if L(ist) < max(L)
        R=R+(L(ist)-max(L))*pfqn_qzgbup(L,N-1,Z,ist);
    end
end
X=2*N/(R+sqrt(R^2-4*Z*max(L)*N));
end