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!