%Analizo la distribucion de numeros aleatorios


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

samples = rand(10,1)


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

%son pocas muestras, veamos mas

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

hist(samples,[0:1e-2:1])

%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)

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

hist(samples,[0:1e-2:1])
media = mean(samples)

sigma = std(samples)


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