ion
above with the single statement</p>
<pre><code>  return mu + sigma * Phi(uniform_rng(normal_cdf(lb, mu, sigma),
                                      normal_cdf(ub, mu, sigma)));</code></pre>
<p>for readability. The names like <code>p</code> indicate probabilities, and
<code>p_lb</code> and <code>p_ub</code> indicate the probabilities of the
bounds. The variable <code>u</code> is clearly named as a uniform variate,
and <code>y</code> is used to denote the variate being generated itself.</p>

</div>
</di