function E = particlebox(N)

  [D,x] = cheb(N);
  D = 2*D;
  L = D^2;  
  E = eig(L(2:end-1,2:end-1));
  E = -sort(E, 'descend');
  E = E(1:4);
