On Wed, 17 Nov 2004, Norbert Nemec wrote: > def gap(A): > w = sort(real(linalg.eigwals(A))) > return min(w[1:]-w[:-1])/2 This would give incorrect result, wouldn't it? Try A with all eigenvalues being imaginary, for example. Pearu