We describe three different methods for generating quasi-exactly solva
ble potentials, for which a finite number of eigenstates are analytica
lly known. The three methods are respectively based on (i) a polynomia
l ansatz for wave functions; (ii) point canonical transformations; (ii
i) supersymmetric quantum mechanics. The methods are rather general an
d give considerably richer results than those available in the current
literature.