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