lass="nb">int</span><span class="p">(</span><span class="n">halfm</span><span class="o">*</span><span class="nb">datan</span><span class="p">(</span><span class="mf">1.d0</span><span class="p">)</span><span class="o">/</span><span class="mf">8.d0</span><span class="p">)</span><span class="w"> </span><span class="o">+</span><span class="w"> </span><span class="mi">5</span><span class="w">          </span><span class="c">! compute multiplier and increment for linear congruential method</span>
<span class="w">