an class="comment">% Compute principal components</span>
0030 <span class="keyword">if</span> nargout == 1
0031     <span class="keyword">for</span> f=1:N
0032         Gxxf=squeeze(Gxx(f,:,:));
0033         PC(f,:)=svd(Gxxf);
0034         A=PC(f,:);
0035     <span class="keyword">end</span>
0036 <span class="keyword">elseif</span> nargout &gt;= 2
0037     VCxx=zeros(N,R,R);
0038     <span class="keyword">for</span> f=1:N
0039         Gxxf=squeeze(Gxx(f,:,:));
0040         [U1,S,U2]=svd(Gxxf);
0041         P