In this paper we discuss time integrators for nonlinear differential equations. In recent years, splitting approaches have become an important tool for reducing the computational time needed to solve differential equations. Moreover, nonlinearity is a challenge to splitting schemes, while one has to extend the exp-functions in terms of a nonlinear Magnus expansion. Here we discuss a novel extension of the so-called multiproduct expansion methods, which is used to improve the standard Strang splitting schemes as to their nonlinearity. We present an extension of linear splitting schemes and concentrate on nonlinear systems of differential equations and generalise in this respect the recent MPE method; see (Chin and Geiser, 2011). Some first numerical examples, of rigid body problems, are given as benchmarks.