..</span>
0045                       xz(i+1)-z(k),xz(i+1)-z(k+1),M(i,k),Cov(i,k));
0046         <span class="keyword">end</span>
0047             <span class="keyword">end</span>
0048         <span class="keyword">end</span>    
0049     <span class="keyword">else</span>
0050         <span class="keyword">for</span> k = 1:size(M,2),
0051             <span class="keyword">for</span> i = 1:size(M,1),
0052         <span class="keyword">if</span> xz(i)==0,
0053                   fprintf(fid,<span class="string"