clear all
% tempo finale
T = 100; % dovrebbe essere infinito
% numero di simulazioni
N = 100;
for i = 1:N
  [S,Q,a] = codagg1(T);
  SS(i) = S(end);
  QQ(i) = Q(end);
end
sum(QQ)/N
a*sum(SS)/N