N = 50;
  r = -i./[1:N]; c = [-i pi i./[1:N-2]];
  T = toeplitz(c,r);
  opts.npts = 50;
  opts.ax = [-4 7 -6 3];
  opts.levels = [-5:-1];
  eigtool(T, opts)