function [ R ] = aleatorio_com_probabilidade( Ps, Ns_R )
% [ R ] = aleatorio_com_probabilidade( Ps, Ns_R )

if nargin<2
    Ns_R = 1;
end

if Ns_R==1
    R = sum( rand > cumsum( Ps )/sum( Ps ) )+1;
else
    RANDOMICOS = rand( Ns_R, 1 );
    R = ones( Ns_R, 1 );
    Lims_Ps = cumsum( Ps )/sum( Ps ); 
    for qregra = 1:(length(Ps)-1)
        R = R + (RANDOMICOS>Lims_Ps(qregra));
    end
end