%  rando.m genera una variable aleatoria en 1, 2, ..., n para una
%  distribución dada. Usa la F ingenua.

function [index] = rando(p) % la entrada es un vector de prbabilidades y devuelve la realizacion de una v.a generada con esa ditribución
u = rand; % genera la variable uniforme
i = 1; 
s = p(1);

while ((u > s) & (i < length(p))),
  i=i+1;
  s=s+p(i);
end

index=i;

