Y = function(N, mu, sigma){ sample = exp(rnorm(N, mu, sigma)) return(sample) } mu_range = seq(0, 5, by = 0.01) for (i in mu_range){ sample_mu = Y(10, i, 2)