This paper addresses an attraction–repulsion chemotaxis system governed by Neumann boundary conditions within a bounded domain (Formula presented.) that has a smooth boundary. The primary focus of the study is the chemotactic response of a species (cell population) to two competing signals. We establish the existence and uniqueness of a weak solution to the system by analyzing the solvability of an approximate problem and utilizing the Leray–Schauder fixed-point theorem. By deriving appropriate a priori estimates, we demonstrate that the solution of the approximate problem converges to a weak solution of the original system. Additionally, we conduct computational studies of the model using the finite element method. The accuracy of our numerical implementation is evaluated through error analysis and numerical convergence, followed by various numerical simulations in a two-dimensional domain to illustrate the dynamics of the system and validate the theoretical findings. © 2025 Elsevier B.V., All rights reserved.