%Creates all 39 Q matrices, named qXXXX.

p = who('sc*');
for k = 1:length(p),
	seq = char(p(k));
	seq = seq(3:6);
	A = eval(char(p(k)));
	[vA, vpA] = eig(A);
	[vpA, s] = sort(diag(vpA)); %sort eigenvalues
	vA = vA(:,s);               %sort eigenvectors the same way
	vAA = reorder(vA);
	eval(['q' seq '=vAA;']);
end;
clear vA vpA vAA seq p A s k