span class="p">)</span>    <span class="kd">::</span> <span class="n">relax</span>            <span class="c">!! *relaxation parameter*</span>
   <span class="kt">real</span><span class="p">(</span><span class="nb">kind</span><span class="o">=</span><span class="n">R8</span><span class="p">)</span>    <span class="kd">::</span> <span class="n">eps</span>              <span class="c">!! *error for convergence*</span>
   <span class="kt">integer</span><span class="p">(</span><span class="nb">kind</span><span 