A method is derived for computing the invariant zeros of a general non-square matrix-valued function of a complex variable in terms of the parameters of the corresponding state-space system model. Based on this a method is developed for choosing a ‘ squaring down ’ output matrix to give a desired set of zero locations for a system whose states are all accessible. This technique is then extended to deal with the situation when all the states are not accessible.