function [probabilidade] = prob(p,n,k,N)
%probabilidade e o return output
%ver se queremos < ou >
  % p -> probabilidade
  % n -> numero de tries
  % k -> numero de caras pretendidas
  % N -> numero de experiencias
  
    lancamentos = (rand(n,N) < p);
    sucessos= sum(lancamentos)==k;
    probabilidade= sum(sucessos)/N;

end