  We present an algorithm for computing the set of torsion points satisfying a
given system of multivariate polynomial equations. Its complexity is
quasilinear in the logarithm of the degree of the input equations and
exponential in their number of non zero terms and variables.
