  We consider three different approaches to analyze the quantum mechanical
problems in multi-well potentials:
  i) the standard matrix diagonalization technique in the basis sets of
harmonic oscillator eigenfunctions or plain waves;
  ii) the spectral method, which allows to reconstruct the spectrum and
stationary functions based on the time-dependent solution of the Schr\"odinger
equation;
  iii) approximations with exact solutions obtained by the supersymmetric
quantum mechanics technique.
  The latter approach proves to be the most promising as it gives a unique
possibility to include the specific multi-well features of the problem directly
in the calculation procedure.
