</span></span>
<span class="r-in"><span><span class="co">## Dispersion model</span></span></span>
<span class="r-in"><span><span class="va">sim1</span> <span class="op">&lt;-</span> <span class="kw">function</span><span class="op">(</span><span class="va">nfac</span><span class="op">=</span><span class="fl">40</span>, <span class="va">nt</span><span class="op">=</span><span class="fl">100</span>, <span class="va">facsd</span><span class="op">=</span><span class="fl">0.1</span>, <span class="va">tsd</span><s