  We propose a general method for constructing quasi-exactly solvable
potentials with three analytic eigenstates. These potentials can be real or
complex functions but the spectrum is real. A comparison with other methods is
also performed.
