% Glivenko-Cantelli

clear;clf;
N=input('Nombre de realisations:');
l=input('Parametre de la loi exponentielle:');
xmax=input('Borne sup de l intervalle considere:');
h=input('Taille de la subdivision:');
Z=rexpweib(N,1);
Femp=[];
Ftheo=[];
for i=1:h
   Femp=[Femp sum(Z<=i/h*xmax)/N];
   Ftheo=[Ftheo 1-exp(-l*i/h*xmax)];
end
plot(Femp,'b')
title('Theoreme de Glivenko-Cantelli')
xlabel('Nombre de realisations:')
ylabel('Fonctions de repartition:')
hold on
plot(Ftheo,'r--')
legend('Empirique','Theorique')
hold off



