init
chegada=newevent
chegar = 'negexp(inter_arrival_time,seed)'; 
inter_arrival_time=5;
t_final=100;
seed = 13; 
t=0:t_final;
%Recolher o nš de eventos em 1000 runs;

for i=1:1000;

tt=0;
while tt < t_final
    [inter_time,seed]=eval(chegar);
    [tt]=[tt,tt(length(tt))+inter_time];
end;

n_events(i)=length(tt)-1;
end;