  We look at BPS systems involving two interacting Sine-Gordon like fields both
when one of them has a kink solution and the second one either a kink or an
antikink solution. The interaction between the two fields is controlled by a
parameter $\lambda$ which has to satisfy $| \lambda|< 2$. We then take these
solitonic static solutions (with solitons well localised) and construct from
them systems involving two solitons in each field (kinks and antikinks) and
then use them as initial conditions for their evolution in Lorentz covariant
versions of such models. This way we study their interactions and compare them
with similar interactions involving only one Sine-Gordon field. In particular,
we look at the behaviour of two static kinks in each field (which for one field
repel each other) and of a system involving kinks and anti-kinks (which for one
field attract each other) and look how their behaviour depends on the strength
of the interaction $\lambda$ between the two fields. Our simulations have led
us to look again at the static BPS solutions of systems involving more fields.
We have found that such ostensibly 'static' BPS solutions can exhibit small
motions due to the excitation of their zero modes. These excitations arise from
small unavoidable numerical errors (the overall translation is cancelled by the
conservation of momentum) but as systems of two or more fields have more that
one zero mode such motions can be generated and are extremely small. The energy
of our systems has been conserved to within $10^{-5}\%$.
