%4. PDF EXPONENCIAL.
r3 = rand(10000,1); % vector de 10000 números con cota superior 1.
meanexp=0.5; %Escogemos la media que deseemos
e = -meanexp*log(1-r3);
%Medimos varianza, esperanza y coeficiente de variación
Mexp = mean(e);
Vexp = var(e);
CVexp= std(e)/Mexp; %Tiene que dar 1 o cercano a 1 por que es una exponencial.
%Obtenemos la pdf
figure(4)
[ht3,x3]=hist(e,50);
bar(x3,ht3)
title("Histograma de la exponencial");
area3=(x3(2)-x3(1))*sum(ht3);
pdf3=ht3./area3;
figure(5)
bar(x3,pdf3);
title("PDF obtenida de la exponencial");
%Obtenemosmla exponencial real
ereal = (1/meanexp)*exp(-x3*(1/meanexp));
figure(6)
plot(x3, pdf3, x3, ereal)
title("Compración pdf práctica/teórica");