Abstract:
An effective general algorithm for the simulation of one-dimensional probability densities defined on a finite interval is presented. The two-side rejection algorithm has been constructed. It uses piecewise linear interval approximations of structurally simple functions. The algorithm has been investigated in detail for the beta-distribution with non-negative parameters. A modification of the algorithm for infinite densities and infinite intervals is proposed.