On pseudorandom numbers from multivariate polynomial systems

We bound exponential sums along the orbits of essentially arbitrary multivariate polynomial dynamical systems, provided that the orbits are long enough. We use these bounds to derive nontrivial estimates on the discrepancy of pseudorandom vectors generated by such polynomial systems. We generalize several previous results and in particular suggest a new approach that eliminates the need to control the degree growth of the iterations of these polynomial systems, which has been an obstacle in all previous approaches.