Quintic B-spline collocation approach for solving generalized BlackScholes equation governing option pricing
We develop a numerical algorithm for solving generalized BlackScholes partial differential equation, which arises in European option pricing. The method comprises the horizontal method of lines for time integration and ? ? -method, ??[1/2,1] ? ? [ 1 / 2 , 1 ] (?=1 ? = 1 corresponds to the back-ward Euler method and ?=1/2 ? = 1 / 2 corresponds to the CrankNicolson method) to discretize in temporal direction and the quintic B-spline collocation method in uniform spatial direction. The convergence analysis and stability of proposed method are discussed in detail, it is justifying that the approximate solution converges to the exact solution of orders O(k+h3) O ( k + h 3 ) for the back-ward Euler method and O(k2+h3) O ( k 2 + h 3 ) for the CrankNicolson method, where k k and h h are mesh sizes in the time and space directions, respectively. The proposed method is also shown to be unconditionally stable. This scheme applied on some test examples, the numerical results illustrate the efficiency of the method and confirm the theoretical behavior of the rates of convergence. Results shown by this method are found to be in good agreement with the known exact solutions. The produced results are also seen to be more accurate than some available results given in the literature.
quintic b-spline collocation generalized blackscholes equation option pricing convergence stability