Función de densidad de probabilidad no uniforme

[asuma]

Hola.

Recurro a ustedes por lo siguiente:

Quiero programar una función que me reciba como entrada una función de densidad de probabilidad. La función es discreta, por ejemplo:

p(1) = 1/2
p(2) = 1/4
p(3) = 1/4

Obviamente la notación es diferente... pero la idea es esa: quiero hacer una función que me devuelva "1" con una probabilidad de 1/2. Que me devuelva "2" con una probabilidad de 1/4, y que me devuelva "3" con una probabilidad de 1/4 (esto es un ejemplo).

Requiero de hacer la función de tal manera que los datos p(x) sean una entrada....

Estaba pensando en llenar un arreglo con (para este ejemplo):

dos "1"s
un "2"
un "3"

y utilizar Mat.Random para seleccionar un elemento de dicho arreglo... Pero me gustaría escuchar su opinión, o si conocen algún método se los agradecería mucho.

Saludos!

[asuma]

Uuuuhhhh ya encontré una solución más simple:

No es necesario utilizar un arreglo, pues basta con asignar los intervalos de valores:

[1,2] al 1
[3] al 2
[4] al 3

Utilizamos Math.Random para elegir aleatoriamente valores entre 1 y 4. Vemos en que intervalo se encuentra el valor generado y entregamos el número asociado a dicho intervalo...

Incluso podemos normalizar el tamaño total de los intervalos a, por ejemplo, 1 (o cualquier otro entero).

Saludos!!