Dynamical system theory is a widely used technique to analyze the asymptotic behavior of a cosmological model. In this method, the equations that describes the dynamics are written in terms of dimensionless variables to make up a set of autonomous first-order differential equations. Then, the asymptotic dynamics of the model are encoded in the fixed points of the autonomous set. Usually, these points are analytical expressions for the variables in terms of the parameters of the model, allowing us to constrain its parameter space. However, analytical treatment could be impossible in some cases. In this work, we give an example of a dark energy (DE) model with no analytical fixed points, a tachyon field coupled to a vector field in an anisotropic Bianchi I background, and propose a numerical description of the parameter space, which allows us to find the bifurcation curves of possible accelerated attractors of the system. This work could serve as a template for numerical analysis of highly complicated dynamical systems.