er of frequencies</span>
0022 <span class="keyword">if</span> nfreq~=length(findx); error(<span class="string">'frequency information (last two arguments) inconsistent'</span>); <span class="keyword">end</span>;
0023 H=fft(tapers,nfft,1);  <span class="comment">% fft of tapers</span>
0024 H=H(findx,:); <span class="comment">% restrict fft of tapers to required frequencies</span>
0025 w=2*pi*f; <span class="comment">% angular frequencies at which ft is to be evaluated</span>
0026 Nsp=zeros(1,C); Msp=zer