Within the framework of Einstein’s General Relativity we study strange quark stars assuming an interacting equation-of-state. Taking into account the presence of anisotropies in a sphere made of ultra dense matter, we employ the formalism based on the complexity factor. We integrate the structure equations numerically imposing the appropriate conditions both at the center and at the surface of the stars, thus obtaining interior solutions describing hydrostatic equilibrium. Making use of well-established criteria, we demonstrate that the solutions obtained here are well behaved and realistic. A comparison with another, more conventional approach, is made as well. Our numerical results are summarized in a number of figures.