Let (P, Q) be a C 1 vector field defined in a open subset U ⊂ R2 . We call a null divergence factor a C 1 solution V (x, y) of the equation P ∂V + Q ∂V = ∂P + ∂Q V . In previous works ∂x ∂y ∂x ∂y it has been shown that this function plays a fundamental role in the problem of the center and in the determination of the limit cycles. In this paper we show how to construct systems with a given null divergence factor. The method presented in this paper is a generalization of the classical Darboux method to generate integrable systems.