Premium
Numerical solution of systems of partial integral differential equations with application to pricing options
Author(s) -
Yousuf Muhammad
Publication year - 2018
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.22244
Subject(s) - mathematics , stability (learning theory) , convergence (economics) , partial derivative , numerical analysis , numerical stability , jump diffusion , partial differential equation , consistency (knowledge bases) , penalty method , jump , mathematical optimization , mathematical analysis , computer science , physics , geometry , quantum mechanics , machine learning , economics , economic growth
We introduce and analyze a strongly stable numerical method designed to yield good performance under challenging conditions of irregular or mismatched initial data for solving systems of coupled partial integral differential equations (PIDEs). Spatial derivatives are approximated using second order central difference approximations by treating the mixed derivative terms in a special way. The integral operators are approximated using one and two–dimensional trapezoidal rule on an equidistant grid. Computational complexity of the method for solving large systems of PIDEs is discussed. A detailed treatment for the consistency, stability, and convergence of the proposed method is provided. Two asset American option under regime–switching with jump–diffusion model when solved using a penalty term, leads to a system of two dimensional PIDEs with mixed derivatives. This model involves double probability density function which brings more challenges to the numerical solution in already a complicated partial integral differential equation. The complexity of the dense jump probability generator, the nonlinear penalty term and the regime–switching terms are treated efficiently, while maintaining the stability and convergence of the method. The impact of the jump intensity and other parameters is shown in the graphs. Numerical experiments are performed to demonstrated efficiency, accuracy, and reliability of the proposed approach.