We use Monte Carlo methods to estimate the expectation over the posterior distribution q(z|x) and q(y|x) inside the objective function 14.
we draw Gumbel noise and Gaussian noise one at a time to compute the latent variables y and z.