Dividend paying European stock options are modeled using a time-fractional Black–Scholes (tfBS) partial differential equation (PDE). The underlying fractional stochastic dynamics explored in this work are appropriate for capturing market fluctuations in which random fractional white noise has the potential to accurately estimate European put option premiums while providing a good numerical convergence. The aim of this paper is two fold: firstly, to construct a time-fractional (tfBS) PDE for pricing European options on continuous dividend paying stocks, and, secondly, to propose an implicit finite difference method for solving the constructed tfBS PDE. Through rigorous mathematical analysis it is established that the implicit finite difference scheme is unconditionally stable. To support these theoretical observations, two numerical examples are presented under the proposed fractional framework. Results indicate that the tfBS and its proposed numerical method are very effective mathematical tools for pricing European options.