=nt_cov_lags(x,y,shifts)</a>
0002 <span class="comment">%[C,tw,m]=nt_cov_lags(x,y,shifts) - covariance of [x,y] with lags</span>
0003 <span class="comment">%</span>
0004 <span class="comment">%  C: covariance matrix (3D if length(lags)&gt;1)</span>
0005 <span class="comment">%  tw: total weight</span>
0006 <span class="comment">%  m: number of columns in x</span>
0007 <span class="comment">%</span>
0008 <span class="comment">%  x,y: data matrices</span>
0009 <span class="comment">%  shifts: positive shift m