A fractional Black–Scholes (B–S) model of order \(\gamma \) is obtained when the change is subjected to price under the basic fractal communication system. In this work, we have presented the exact solution of time fractional B–S equations determining European options using a new method known as the Sumudu transform iterative method. The solution is given in a convergent series form demanding no additive assumptions or discretization. Four examples are solved using the proposed method and the Mittag-Leffler function appears in the closed form solutions of these examples. Numerical results prove the accuracy and easy to use approach of the method for solving time fractional B–S equation. At the same time, numerical experiments indicate that the time fractional B–S equation is more competent with the actual financial market scenarios. The prime focus of this study is to minimize the difficulty of computational work related to other conventional methods. Moreover, it is an efficient technique to solve general time fractional order partial differential equations arising in many fields of engineering, economics, physics and numerous other disciplines.
