This monograph primarily concerns the use of power method polynomials in the context of simulating univariate and multivariate nonnormal distributions with specified cumulants and correlation matrices.