x=find(dtmp&gt;=min(t)&amp;dtmp&lt;=max(t));
0037      <span class="keyword">if</span> ~isempty(indx); dtmp=dtmp(indx);
0038      <span class="keyword">end</span>;
0039   <span class="keyword">end</span>;
0040   Nsp(ch)=length(dtmp);
0041   Msp(ch)=Nsp(ch)/length(t);
0042   <span class="keyword">if</span> Msp(ch)~=0;
0043       data_proj=interp1(t',tapers,dtmp);
0044       exponential=exp(-i*w'*(dtmp-t(1))');
0045       J(:,:,ch)=exponential*data_proj-H*Msp(ch);
0046   <span class="keyword">else</