% Capacity vs. arrival rate given QoS
a = 0:0.01:1; %arrival rate [p/us]
la = length(a);
z = 10; %longest time to wait [us]
L = 1; % average packet length [kb]
C = 0; % link speed [Gbps or kbpus]
Pz_obj = [0.01 0.1 0.2];
lPz = length(Pz_obj);
LS = zeros(la,lPz); % variable to store C(lambda)
for j=1:lPz
    for l = 1:la
        c = C;
        lambda = a(l);
        Pz = Wait(c,lambda,L,z);
        while Pz > Pz_obj(j)
            c = c + 0.001;
            Pz = Wait(c,lambda,L,z);
        end
        LS(l,j) = c;
    end
end
figure
f = plot(a,LS);
f(1).LineWidth = 1.5;
f(2).LineWidth = 1.5;
f(3).LineWidth = 1.5;
grid on
title('Capacity vs. Arrival rate')
xlabel('Arrival rate (\lambda) [packets/\mus]')
ylabel({'Link speed (C)';'[Gbps]'})
legend({['Pz_{obj} = ',num2str(Pz_obj(1))],['Pz_{obj} = ',num2str(Pz_obj(2))],...
    ['Pz_{obj} = ',num2str(Pz_obj(3))]})