We consider a system of non-linear eikonal equations in one space dimension that describes the evolution of interfaces moving with non-signed strongly coupled velocities. We have recently proven the global existence and uniqueness of viscosity solutions for this system, under a BV estimate. In this paper, we propose a semi-explicit scheme that satisfies the same BV estimate proven in the continuous case, at the discrete level, and we show that a certain linear interpolation of the discrete solution to the scheme converges to a viscosity solution of the main system considered. We also provide some numerical simulations in the case of dislocations dynamics