In this paper, we prove that the weak error between a stochastic differential equation with nonlinearity in the sense of McKean given by moments and its approximation by the Euler discretization with time-step h of a system of N interacting particles is O(1/N + h). We provide numerical experiments confirming this behaviour and showing that it extends to more general mean-field interaction and study the efficiency of the antithetic sampling technique on the same examples.