Wiener systems consist of a linear dynamic system whose output is measured through a static non-linearity. In this paper we study the identification of single-input single-output Wiener systems with finite impulse response dynamics and polynomial output non-linearities. Using multi-index notation, we solve a least squares problem to estimate products of the coefficients of the non-linearity and the impulse response of the linear system. We then consider four methods for extracting the coefficients of the non-linearity and impulse response: direct algebraic solution, singular value decomposition, multi-dimensional singular value decomposition and prediction error optimization.