dd>
<dt>npts<span class="classifier">int</span></dt><dd><p>number of points of tapers</p>
</dd>
</dl>
<p><em>Notes</em></p>
<p>This is a slight variation of Gleen Ierly’s code. What was
mainly done, was to avoid use of common blocks, defining all
variables and performing the numerical integration inside
(previously done by function pssevf).</p>
<p>Exponential convergence rate for analytic functions!  Much faster 
than Romberg; competitive with Gauss integration, without awkward 
weights.</p>
<p>Integrates