clear
tstart=clock;
%script per la generazione di una processo di interarrivo
%con con generatore geometrico
R=100; %realizzazioni
N=100;  %arrivi medi per realizzazione
dt=0.01;
tmi=1.2; %tempo medio di interarrivo
t1=0:dt:tmi*N;
p=dt/tmi;
ta=rand(R,length(t1))<p;
tmin=0;
tmax=tmi*N;
for i=1:R
    ta1{i}=t1(ta(i,:)==1);
    x=ta1{i};
    ti{i}=x(2:length(x))-x(1:length(x)-1);
    tmin=max(tmin,min(ta1{i}));
    tmax=min(tmax,max(ta1{i}));
end
t=tmin+dt:dt:tmax-dt;
for i=1:R 
    for j=1:length(t) 
        tp=t(j)-ta1{i};
        tla(i,j)=min(tp(tp>0));
        tna(i,j)=-max(tp(tp<0));
    end     
end
mtla=mean(tla); %tempo medio dall'ultimo arrivo
mtna=mean(tna); %tempo medio dal prossimo arrivo
for i=1:3
    subplot(3,2,2*i)
    plot(t(t<3*tmin),tna(i,t<3*tmin),'k'); 
    set(gca, 'Xlim', [tmin,3*tmin]); 
    set(gca, 'Ylim', [0,max(ceil(max(max(tna(1:3,t<3*tmin)))),ceil(max(max(tla(1:3,t<3*tmin)))))]);
    if i==1
        title('tempi per il prossimo arrivo (vita residua)');
    end
end

for i=1:3
    subplot(3,2,2*i-1)
    plot(t(t<3*tmin),tla(i,t<3*tmin),'k'); 
    set(gca, 'Xlim', [tmin,3*tmin]);
    set(gca, 'Ylim', [0,max(ceil(max(max(tna(1:3,t<3*tmin)))),ceil(max(max(tla(1:3,t<3*tmin)))))]);
    if i==1
        title('tempi per dall"ultimo arrivo (vita)');
    end    
end
etime(clock,tstart)