In this work, we are concerned with (R, S) – conjugate solutions to coupled Sylvester complex matrix equations with conjugate of two unknowns. When the considered two matrix equations are consistent, it is demonstrated that the solutions can be obtained by utilizing this iterative algorithm for any initial arbitrary [Formula: see text] – conjugate matrices [Formula: see text]. A necessary and sufficient condition is established to guarantee that the proposed method converges to the [Formula: see text] – conjugate solutions. Finally, two numerical examples are provided to demonstrate the efficiency of the described iterative technique.