The fractional Black–Scholes equation has been widely studied by researchers in recent years. In this article, an efficient spectral collocation method based on fractional Pell functions is proposed for solving the time–fractional Black–Scholes equation. We introduce fractional Pell functions using the transformation x→xβ(β>0) on Pell polynomials, and we look for a solution of the model as a linear combination of these functions. Using operational matrices, we approximate the fractional derivative and other terms in a convenient form of the main equation. A system of algebraic equations is obtained by collocating resultant approximate equations. Convergence analysis of the numerical method has been investigated in Sobolev space. Finally, we have demonstrated the capability of the proposed method by considering numerical experiments in the form of tables and figures.