The return map for a planar vector field with nilpotent linear part: a direct and explicit derivation

Using a direct approach the return map near a focus of a planar vector field with nilpotent linear part is found as a convergent power series which is a perturbation of the identity and 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.