Averaging principle for systems of RDEs with polynomial nonlinearities perturbed by multiplicative noise ∗