function kap8_3a()
% table with multiple runs of singleserver4B
N = 6;
seed0 = 2;
tEnd = 10000;

WQ = zeros(N,1);
LQ = zeros(N,1);
A = zeros(N,1);

load_system('singleserver4B');

% generator uses seed, Server seed+1, therefore steps of 2
for I=1:N
  fprintf('Starting round %3d\n', I)
  [~, LQ(I), WQ(I), A(I)] = runSingleserver4B(seed0+2*I, tEnd);
end

% show statistics
OUT = [(1:N)', LQ, WQ, A];
fprintf('Lauf   LQ        WQ        A\n')
fprintf('%2d %8.4f  %8.4f  %8.4f\n', OUT')

% compute theoretical values
model = 'singleserver4B';
la = str2double(get_param([model, '/Generator'], 'lambda'));
mu = str2double(get_param([model, '/Server'], 'mu'));
sig = str2double(get_param([model, '/Server'], 'sig'));
rho = la/mu;
A = rho;
LQ = rho^2*(1 + sig^2*mu^2)/(2*(1 - rho));
WQ = la*(1/mu^2 + sig^2)/(2*(1 - rho));
fprintf('th %8.4f  %8.4f  %8.4f\n', LQ, WQ, A)
fprintf('\nla=%6.4f, mu=%6.4f, sig=%6.4f\n', la, mu, sig)
