This paper applies the extrapolated homotopy method to solve large sparse linear systems. An optimal extrapolation parameter obtained in this paper makes the approximate solution convergent rapidly. Numerical results demonstrate that the homotopy method embedded with the optimal extrapolation parameter is more accurate than the direct and iterative numerical methods.