%
% Simulation d'un processus de renouvellement.
%

clf
clear
l=input('Parametre de la loi exponentielle, l= ? ');
T1=rand(1);
X(1)=0;
X(2)=T1;
for i=2:100
  X(i+1)=X(i)+rexpweib(1,l);
end
N(1:101)=0:100;
XX(1:2:201)=X(1:101);
XX(2:2:200)=X(2:101);
NN(1:2:201)=N(1:101);
NN(2:2:200)=N(1:100);
plot([0.5 0.5+(100/l)],[1,101],'r--')
line(XX,NN)
xlabel('t')
ylabel('Nt')

