  We present $\Lambda$CDM $N$-body cosmological simulations in the framework of
a static general scalar-tensor theory of gravity. Due to the influence of the
non-minimally coupled scalar field, the gravitational potential is modified by
a Yukawa type term, yielding a new structure formation dynamics. We present
some preliminary results and, in particular, we compute the density and
velocity profiles of the most massive group.
