title: How do I best simulate an arbitrary univariate random variate using its probability function?
body: <p>In R, what's the best way to simulate an arbitrary univariate random variate if only its probability density function is available?</p>