Most recent qualitative models for financial assets assume that the dynamics of underlying equity prices follows a jump or Lévy process. It has been evident that some most intricate characteristics of such dynamics can be captured by CGMY and KoBoL procedures. The prices of financial derivatives with such models satisfy fractional partial differential equations or partial integro-differential equations. This study focuses at aforementioned fractional equations and discretizes them via a monotone Crank-Nicolson procedure. A spatial extrapolation strategy is introduced to ensure an overall second-order accuracy in approximations. Preconditioned conjugate gradient normal residual methods are incorporated for solving resulted linear systems. Numerical examples are given to illustrate the accuracy and efficiency of the novel computational approaches implemented.