[D,x] = cheb(18);       % cheb.m from Trefethen's "Spectral Methods in MATLAB"
  A = D(2:end,2:end); 
  opts.npts=50;
  opts.ax = [-60 10 -35 35];
  opts.levels = -6:0;
  eigtool(A,opts)