Simulation and modeling - Old Questions

The rejection method for obtaining samples of random numbers forms a given non-uniform distribution works by generating uniform random numbers repeatedly and accepting only those numbers that meet certain conditions.

The rejection method is applied when the probability density function, f(x), has a lower and upper limit to its range, 'a' and 'b', respectively, and an upper bound 'c'. The method can be specified as follows:

1. Find the maximum value c of f(x) on a ≤ x ≤ b.

        f(X)  ≤ c ∀ x ∈ [a, b]

2. Compute two values μ₁, μ₂ of the uniformity distributed variables, both defined on [a, b] = [0, 1].

3. Compute x₀ = a + μ₁(b - a)

4. Compute y₀ = c μ₂

5. If y₀ ≤ f(x₀), accept x₀ as desired output; otherwise reject x₀ and repeat the process with two next values μ₁ &  μ_2.