).*cosinefunc;
0101     <span class="comment">%         subplot(221); plot(s(:,1));</span>
0102     <span class="comment">%         subplot(222); plot(s(:,10));</span>
0103     <span class="comment">%         subplot(223); plot(s(:,100));</span>
0104     <span class="comment">%         subplot(224); plot(s(:,120));</span>
0105     <span class="comment">%         pause</span>
0106             s=trapz(f,s,1)';
0107             ss=s.*conj(s);
0108 <span class="comment">%         plot(tau,s)</span>
010