N=50;

A = GOE(N);
B = GOE(N);

tt = linspace(0,1, 1000);

spec = zeros(N,length(tt));
for j = 1:length(tt)
    spec(:,j) = eig(sqrt(1-tt(j))*A + sqrt(tt(j))*B);
end;

plot(tt, spec, 'k-');
