  The interaction between aphids, ants and ladybirds has been investigated from
an ecological point of view since many decades, while there are no attempts to
describe it from a mathematical point of view. This paper introduces a new
mathematical model to describe the within-season population dynamics in an
ecological patch of a system composed by aphids, ants and ladybirds, through a
set of four differential equations. The proposed model is based on the
Kindlmann and Dixon set of differential equations, focused on the prediction of
the aphids-ladybirds population densities, that share a prey-predator
relationship. The population of ants, in mutualistic relationship with aphids
and in interspecific competition with ladybirds, is described according to the
Holland and De Angelis mathematical model, in which the authors faced the
problem of mutualistic interactions in general terms. The set of differential
equations proposed here is discretized by means the Nonstandard Finite
Difference scheme, successfully applied by Gabbriellini to the mutualistic
model. The constructed finite-difference scheme is positivity-preserving and
characterized by four nonhyperbolic steady-states, as highlighted by the
phase-space and time-series analyses. Particular attention is dedicated to the
steady-state most interesting from an ecological point of view, whose
asymptotic stability is demonstrated via the Centre Manifold Theory. The model
allows to numerically confirm that mutualistic relationship effectively
influences the population dynamic, by increasing the peaks of the aphids and
ants population densities. Nonetheless, it is showed that the asymptotical
populations of aphids and ladybirds collapse for any initial condition, unlike
that of ants that, after the peak, settle on a constant asymptotic value.
