A Computational Scheme for the Interaction between an Edge Dislocation and an Arbitrarily Shaped Inhomogeneity via the Numerical Equivalent Inclusion Method