function SSP=CalcSSP(PCA);
// PCA=[PCA_0001,...,PCA_000n]
n=size(PCA,1);
SSP=eye(n,n);
for n=1:size(PCA,2)
    SSP=SSP-PCA(:,n)*PCA(:,n)';
end;
endfunction;

