In this paper we present a model approximation technique based on N-step-ahead affine representations obtained via Monte-Carlo integrations. The approach enables simultaneous linearization and model order reduction of nonlinear systems in the original state space thus allowing the application of linear MPC algorithms to nonlinear systems. The methodology is detailed through its application to benchmark model examples. (C) 2013 Elsevier Ltd. All rights reserved.