% copyright Thomas Mc Avoy 1998
% Routine to carry out SVD analysis
function c=singu(z)

[mm,n]=size(z);
for j=1:mm;
for m=1:mm;
	c(j,m)=z(j,:)*z(m,:)'/(norm(z(j,:))*norm(z(m,:)));
end;
end;
c=c.*c;
