  A code to evolve boson stars in 3D is presented as the starting point for the
evolution of scalar field systems with arbitrary symmetries. It was possible to
reproduce the known results related to perturbations discovered with 1D
numerical codes in the past, which include evolution of stable and unstable
equilibrium configurations. In addition, the apparent and event horizons masses
of a collapsing boson star are shown for the first time. The present code is
expected to be useful at evolving possible sources of gravitational waves
related to scalar field objects and to handle toy models of systems perturbed
with scalar fields in 3D.
