This work focuses on numerical methods for finding optimal investment, dividend payment, and capital injection policies to maximize the present value of the difference between the cumulative dividend payment and the possible capital injections. The surplus is modeled by a regime-switching jump diffusion process subject to both regular and singular controls. Using the dynamic programming principle, the value function is a solution of the coupled system of nonlinear integro-differential quasi-variational inequalities. In this paper, the state constraint of the impulsive control gives rise to a capital injection region with free boundary, which makes the problem even more difficult to analyze. Together with the regular control and regime-switching, the closed-form solutions are virtually impossible to obtain. We use Markov chain approximation techniques to construct a discrete-time controlled Markov chain to approximate the value function and optimal controls. Convergence of the approximation algorithms is proved. Examples are presented to illustrate the applicability of the numerical methods. (C) 2013 Elsevier Ltd. All rights reserved.