  In this work we study static spherically symmetric solutions of effective
field equations related to local and nonlocal higher-derivative gravity models,
based on the effective delta source approximation. We discuss several
possibilities for the equations of state, and how they influence the solutions.
In particular, we present an equation of state such that the solutions match
the Newtonian-limit metric in both regimes of large and small $r$. A
significant part of the work is dedicated to the study of the regularity of the
solutions and the comparison with the linearized solutions. Explicit solutions
are presented for some cases of local and nonlocal models. The results obtained
here can be regarded as a possible link between the previous researches on
linearized higher-derivative gravity and the solutions of the nonlinear
complete field equations.
