“…There is a considerable literature that reviews numerical simulations of the option price when markets are subject to transaction costs using PDEs, we quote [4,19,22], where several nonlinear Black-Scholes models in the presence of transaction costs have been studied for European and American options, [30] presented a positivity preserving scheme for the Barles' and Soner's model and they studied its convergence. [21] used a grid stretching technique on non-uniform meshes for the computation and [24,25] reviewed nonlinear parabolic equations arising from markets under transaction costs. However, there are only few works for the computation of nonlinear Black-Scholes equations depending on more than one underlying asset, among them we cite [29], they gave a numerical analysis for two dimensional nonlinear Black-Scholes equation for Spread option in incomplete markets and [3] where the authors gave a theoretical and numerical study of multi-asset nonlinear Black-Scholes equation arising from models with general transaction costs.…”