Modify theta_epsilon.rvs() for speed