Based on a complex variable boundary integral equation (CVBIE) suggested previously, this paper provides a numerical solution for the elastic inclusion problem using CVBIE. A dissimilar elastic inclusion is embedded in the infinite matrix. The original problem is decomposed into two problems. One is an interior boundary value problem (BVP) for the elastic inclusion, while the other is an exterior BVP for the matrix with notch. Both problems are connected by conventional boundary integral equations (BIEs) in complex variables. After performing discretization for the coupled BIEs, the inverse matrix technique is suggested to solve the relevant algebraic equations. Based on the properties of some integral operators, three ways for the inverse matrix technique are suggested. Several numerical examples are carried out to prove the efficiency of the suggested method.