This paper presents an implementation of the integral equations with the generalized Neumann kernel to solve numerically the uniquely and the non-uniquely solvable Riemann-Hilbert problems in Jordan regions with smooth boundaries. The non-uniquely solvable problems are made uniquely solvable by requiring their solution to satisfy additional constraints. Two type of constraints are presented. Various test numerical examples are presented. The computational efficiency appears significantly excellent.