In this paper, we present an implicit finite difference method for the numerical solution of the Black–Scholes model of American put options without dividend payments. We combine the proposed numerical method by using a front-fixing approach where the option price and the early exercise boundary are computed simultaneously. We study the consistency and prove the stability of the implicit method by fixing the values of the free boundary and of its first derivative. We improve the accuracy of the computed solution via a mesh refinement based on Richardson’s extrapolation. Comparisons with some proposed methods for the American options problem are carried out to validate the obtained numerical results and to show the efficiency of the proposed numerical method. Finally, by using an a posteriori error estimator, we find a suitable computational grid requiring that the computed solution verifies a prefixed error tolerance.