自作した確率密度による確率分布に従うランダムなデータをn件生成するには？

ある連続型の確率分布を定める、任意の確率密度をPythonで実装した関数 f(x) が与えられるものとします。たとえば標準正規分布だったら以下のようなものです。

python

1from math import pi 2from math import exp 3 4def f(x):5 return (1/(2*pi)**0.5) * exp(-0.5*(x**2))

上記のような何らかの確率密度を実装した関数 f(x) が与えられたときに、この f とデータ件数を表す n をパラメータとして受け取り、f(x) による確率分布に従う乱数をn個生成して、それらを要素として含むリストを返すような関数を実装するにはどうしたらよいでしょうか?

コードでなくても、考え方や参照すべきサイトなどをご教示頂ければ助かります。