clear;clf;
n=input('Taille de l echantillon:');
t=input('Entrer la valeur du parametre:');
X=[];
tchap=[];
for i=1:n
   X=[X probadis([1:t],1/t*ones(t,1))];
   tchap=[tchap 2*sum(X)/i-1];
end
plot(tchap,'b')
title('Estimation de theta')
xlabel('Nombre de realisations')
ylabel('Estimateur theta chapeau')
hold on
plot(t*ones(n,1),'r--')
legend('Empirique','Theorique')
hold off