An integral equation method for the numerical solution of a Dirichlet problem for second-order elliptic equations with variable coefficients