

lambda = 0.01                % set values of parameters
 
mu = 0.1

omega = 1.0

Q = zeros(128,128)           % declare the generator matrix

b = zeros(128,1)             % declare the solution vector

pclan4trans                  % read in the entries for Q

Qt = Q'			     % transpose the generator matrix

for i=1:128                  % replace the last equation by
   Qt(128,i) = 1.0           % the normalisation constant in
end                          % both Qt and b
                                          
b(128) = 1.0                             

p = linsolve(Qt,b)           % solve to find the unknowns, the
			     % steady state distribution

numpclan4                    % derive the mean number of data packets at PC1

xputpclan4                   % derive the mean throughput of data packets at PC1

w = N/X - 1/mu               % calculate waiting time as the residence 
                             % time (N/X) minus the transmission time (1/mu)
