This function simulates cases from an exponential distribution. The user may
specify the rate parameter r, the sample size, and the upper bound of the
survival time. If the rate parameter is 0, then this function defaults to the
uniform distribution.
See also
Other simulate:
simulate_gillespie(),
simulate_secondary(),
simulate_uniform_cases()
