Abstract: The mechanical properties and service life of materials may be significantly influenced by the interaction between inhomogeneities and dislocations in engineering structures. Previous analytical studies on inhomogeneity-dislocation interactions have been concerned primarily on some special inhomogeneity shapes (e.g. circle and ellipse). On the other hand, the involved singularity issue in dislocation studies is challenging, even intractable for commercial finite element software. Using the numerical equivalent inclusion method (NEIM) in conjunction with the Fast Fourier Transforms technique, this work presents an effective computational scheme for evaluating the interaction energy between an edge dislocation and inhomogeneities. The proposed computational method may successfully circumvent the numerical singularity. The results of norm analyses on relative errors demonstrate that the stress disturbance field caused by the impurity has a great influence on the final solutions, especially when the dislocation is located in the neighborhood of the inhomogeneity. The proposed method in this work shows excellent numerical convergence and stability, and appears to be convenient and efficient for handling arbitrarily shaped inhomogeneities interacting with an edge dislocation.