  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.
