      </span><span class="n">ic</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">2</span><span class="o">*</span><span class="nb">int</span><span class="p">(</span><span class="n">halfm</span><span class="o">*</span><span class="p">(</span><span class="mf">0.5d0</span><span class="o">-</span><span class="nb">dsqrt</span><span class="p">(</span><span class="mf">3.d0</span><span class="p">)</span><span class="o">/</span><span class="mf">6.d0</span><span class="p">