This paper presents a numerical method for solving the Dirichlet problem with discontinuous coefficients in bounded multiply connected regions. The method is based on reducing the problem of solving the Dirichlet problem with discontinuous coefficients to a problem of solving Dirichlet problem with continuous coefficients. The Dirichlet problem with continuous coefficients is then solved using a combination of a uniquely solvable boundary integral equation with generalized Neumann kernel and the Fast Multipole Method. Numerical results and comparison are given to illustrate the efficiency of the suggested method. Received: March 5, 2014 c © 2014 Academic Publications, Ltd. url: www.acadpubl.eu Correspondence author 448 M. Aghaeiboorkheili, A.H.M. Murid AMS Subject Classification: 30E25, 31B10