We propose a simulation algorithm for non-causal vector autoregressive moving average (VARMA) processes. The algorithm is based on the Jordan canonical form of the companion matrix in the state space representation. We illustrate its performance for a non-causal VARMA(2, 2) process. (C) 2015 Elsevier B.V. All rights reserved.