="l00096" name="l00096"></a><span class="lineno">   96</span><span class="keywordtype">        real</span>(r8), <span class="keywordtype">parameter</span> :: dt = 1.0e-2</div>
<div class="line"><a id="l00097" name="l00097"></a><span class="lineno">   97</span>        </div>
<div class="line"><a id="l00098" name="l00098"></a><span class="lineno">   98</span>        f(x) = (1+x**2/real(dof, r8))**(-(dof+1)/2.0_r8)*gamma((dof+1)/2.0_r8)/sqrt(3.14*dof)/gamma(dof/2.0_r8)</div>
<div class="line"><a id="l00099" na