%ex1
N = 1e5; %Numero de experiencias
p = 0.5; %Probabilidade de cara
k = 2;   %numero de caras que queremos que calhem
n = 3;   %numero de lancamentos por exeperiencia

probEx1=prob(p,n,k,N)

%ex2
N = 1e5;  %Numero de experiencias
p = 0.5;  %Probabilidade do acontecimento (sair caras)
k = 6;    %Numero de vezes que queremos que saiam caras
n = 15;    %Numero total de lancamentos

probEx2=prob(p,n,k,N)


%ex3
N = 1e5;  %Numero de experiencias
p = 0.5;  %Probabilidade do acontecimento (sair caras)
n = 15;    %Numero total de lancamentos
%k = Numero de vezes que queremos que saiam caras
probEx3=0;
for k= 6:15
    probEx3 = prob(p,n,k,N) + probEx3;
end

probEx3

%ex4a)
disp("")

p=0.5;
n=20;
N=1e5;

for k = 0: n
  fprintf("P(Sair %d caras em %d lançamentos) = %d \n",k,n,prob(p,n,k,N));
end

disp("")
n=40;

for k = 0: n
  fprintf("P(Sair %d caras em %d lançamentos) = %d \n",k,n,prob(p,n,k,N));
end

disp("")
n=100;

for k = 0: n
  fprintf("P(Sair %d caras em %d lançamentos) = %d \n",k,n,prob(p,n,k,N));
end