a>01170           dl(i - lp1, j - lp1) = dl(j - lp1, i - lp1) * sgn
<a name="l01171"></a>01171           sgn = sgn * (- 1)
<a name="l01172"></a>01172         <span class="keyword">end do</span>
<a name="l01173"></a>01173       <span class="keyword">end do</span>
<a name="l01174"></a>01174   
<a name="l01175"></a>01175       <span class="comment">! Fill in right quarter of matrix.</span>
<a name="l01176"></a>01176   
<a name="l01177"></a>01177       <span class="keyword">do</span> i = l + 2, 2 * l + 1
<a nam