This study thoroughly investigates the efficiency of advanced numerical extrapolation methods aimed at enhancing the convergence of vector sequences in the realm of mathematical finance. Our focus lies in the application of polynomial extrapolation techniques to calculate finite difference solutions for the Black-Scholes (BS) equation–an indispensable model in options pricing. The performance of our algorithms undergoes rigorous evaluation through a comprehensive analysis involving both simulated and real-world data. Notably, our experiments uncover that a stochastic scheme, incorporating two extrapolation strategies and a random relaxation parameter, outperforms other proposed methods, excelling in both convergence and stability metrics. Our findings underscore the potential of this numerical extrapolation method to enhance the efficiency of financial calculations, particularly in the realm of option pricing. This innovation holds promise for refining financial models and addressing specific challenges within the field of mathematical programming, providing effective solutions to the primary computational bottlenecks commonly encountered in financial decision-making.