Several approximate methods have been used to find approximate solutions of elliptic systems of first order equations. One common method is the Newton imbedding approach, i.e. the parameter extension method. In this article, we discuss approximate solutions to discontinuous Riemann-Hilbert boundary value problems, which have various applications in mechanics and physics. We first formulate the discontinuous Riemann-Hilbert problem for elliptic systems of first order complex equations in multiply connected domains and its modified well-posedness, then use the parameter extensional method to find approximate solutions to the modified boundary value problem for elliptic complex systems of first order equations, and then provide the error estimate of approximate solutions for the discontinuous boundary value problem.