  The return map for planar vector fields with nilpotent linear part (having a
center or a focus and under an assumption generically satisfied) is found as a
convergent power series whose terms can be calculated iteratively. The first
nontrivial coefficient is the value of an Abelian integral, and the following
ones are explicitly given as iterated integrals built with algebraic functions.
