Numerical Solution of Fractional Black-Scholes Equation by Using the Multivariate Padé Approximation

In this study, a new application of multivariate Padé approximation method has been used for solving European vanilla call option pricing problem. Padé polynomials have occurred for the fractional Black-Scholes equation, according to the relations of "smaller than", or "greater than", between stock price and exercise price of the option. Using these polynomials, we have applied the multivariate Padé approximation method to our fractional equation and we have calculated numerical solutions of fractional Black-Scholes equation for both of two situations. The obtained results show that the multivariate Padé approximation is a very quick and accurate method for fractional Black-Scholes equation. The fractional derivative is understood in the Caputo sense.