clc
clear
%N=input('ingrese el valor del numero de repeticiones \n') %ingresamos el valor de las repeticiones
v=rand(1,N) %vector de numeros aleatorios, son las repeticiones de los eventos
X=zeros(1,N); % almacenar el numero de veces que se repiten los eventos
p=1/10; %es la probilidad tomada de la definicion clasica
for i=1:N; %la condicionante del numero de repeticiones
if(v(i)<p)
X(i)=1;
else if (p<=v(i) && v(i)<(2*p)) % dividimos en intervalos por cada evento
X(i)=2;
%que este pudiera suceder N numero de veces
else if ((2*p) && v(i)<(3*p))
X(i)=3;
else if ((3*p) && v(i)<(4*p))
X(i)=4;
else if ((4*p) && v(i)<(5*p))
X(i)=5;
else if ((5*p) && v(i)<(6*p))
X(i)=6;
else if ((6*p) && v(i)<(7*p))
X(i)=7;
else if ((7*p) && v(i)<(8*p))
X(i)=8;
else if ((8*p) && v(i)<(9*p))
X(i)=9;
else if ((9*p) && v(i)<(10*p))
X(i)=10;
end
end
end
end
end
end
end
end
end
end
end
buno=find(X==1);
bdos=find(X==2);
btres=find(X==3); %Escribiendo el cuaquie nombre del numero de bola saldra el valor de veces que se repite
bcuatro=find(X==4);
bcinco=find(X==5);
bseis=find(X==6);
bsiete=find(X==7);
bocho=find(X==8);
bnueve=find(X==9);
bdiez=find(X==10);
b1=length(buno);
b2=length(bdos);
b3=length(btres);
b4=length(bcuatro); %si escribimos cualquier valor de f no dara el numero de veces que ha salido el numero
b5=length(bcinco);
b6=length(bseis);
b7=length(bsiete);
b8=length(bocho);
b9=length(bnueve);
b10=length(bdiez);
x(1)=b1/N;
x(2)=b2/N; % el numero de veces que fue extraida una bola del 1 al 10
x(3)=b3/N; %probabilidad de cada evento apartir de la frecuecia relativa
x(4)=b4/N; %b6 es el numero de veces que ocurrio el evento y N el numero total
x(5)=b5/N;
x(6)=b6/N;
x(7)=b7/N;
x(8)=b8/N;
x(9)=b9/N;
x(10)=b10/N;
y=1:10;
figure(4)
bar(y,x),grid