%Analizo la distribucion de numeros aleatorios


%genero 10 numeros reales al azar entre 0 y 1, con una
%distribucion uniforme


samples = rand(10,1)
%NB el 10,1 en la funcion rand es para que me de un vector
%columna con 10 elementos.

%la funcion hist hace el histograma de los datos, es decir cuenta
%cuantos caen en cada uno de los intervalos(bins) indicados.
hist(samples); legend('10 nros aleatorios');

%son pocas muestras, veamos mas

samples = rand(1e5,1); %el ; al final evita que me muestre 10^5 numeros aleatorios

figure; %para que haga una figura nueva y no borre la vieja
hist(samples,[0:1e-2:1]); legend('10^5 nros aleatorios');

%tambien vemos la media y la desviacion std

%la media debe ser 1/2
media = mean(samples)

%sigma = sqrt(<^2> - <>^2), debe ser 1/(2 Sqrt(3)) ~ 0.288675134594813)
sigma = std(samples)

%si quieren guardar el estado del generador de 
%numeros aleatorios
v=rand("state");



%podemos repetir los mismo con otras distribuciones, por ejemplo
%la binomial con parametros N y p

%en Matlab y Octave se necesitan paquetes especiales que puede
%que no tengas instalado

N=500;
p=0.5;
samples=random("bino",N,p,[1 1e5]);

%deberia ser N p
media = mean(samples)

%deberia ser sqrt(N p (1-p))
sigma =  std(samples)
figure;
hist (samples);legend('Binomial - 10^5 nros aleatorios')