%  maxEig   Maximum Eigenvalue of a random matrix
%     Input parameters
%       sd - seed for uniform random generator
%       dim - size of the square matrix (should be odd)
%     Output value
%       maxe - maximum real eigenvalue
function maxe = maxEig(sd,dim)
  rng(sd);
  ev = eig( randn(dim) );
  maxe = max( ev(imag(ev)==0) )
end