The characteristic equation for the Mueller matrix of an arbitrary nondepolarizing optical system is obtained and solved. Polarization eigenstates for all types of nondepolarizing systems are studied. An efficient method for obtaining eigenvectors of matrices as functions of eigenvalues is described.