Dynamical system theory is a widely used technique in the analysis of
cosmological models. Within this framework, the equations describing the
dynamics of a model are recast in terms of dimensionless variables, which
evolve according to a set of autonomous first-order differential equations. The
fixed points of this autonomous set encode the asymptotic evolution of the
model. Usually, these points can be written as analytical expressions for the
variables in terms of the parameters of the model, which allows a complete
characterization of the corresponding parameter space. However, a thoroughly
analytical treatment is impossible in some cases. In this work, we give an
example of a dark energy model, a scalar field coupled to a vector field in an
anisotropic background, where not all the fixed points can be analytically
found. Then, we put forward a general scheme that provides a numerical
description of the parameter space. This allows us to find interesting
accelerated attractors of the system with no analytical representation. This
work may serve as a template for the numerical analysis of highly complicated
dynamical systems.Comment: 13 pages, 13 figures, 1 table. Changes match the published versio