function R=CalcSSP2(MEG);
R=[];
PCA=[];
for x=1:1000;
  ssp=sprintf('PCA_%04d',x);
  if exists(ssp);PCA=[PCA,eval(ssp)];
  else break;
  end;
end;
if isempty(PCA);return;end;
n=size(PCA,1);
SSP=eye(n,n);
for x=1:size(PCA,2)
    SSP=SSP-PCA(:,x)*PCA(:,x)';
end;
if     size(MEG,1)==n;R=SSP*MEG;
elseif size(MEG,2)==n;R=MEG*SSP;
end;
endfunction;


