Markov chain approximations to scale functions of Lévy processes
arXiv:1310.1737
Abstract
We introduce a general algorithm for the computation of the scale functions of a spectrally negative Lévy process , based on a natural weak approximation of via upwards skip-free continuous-time Markov chains with stationary independent increments. The algorithm consists of evaluating a finite linear recursion with its (nonnegative) coefficients given explicitly in terms of the Lévy triplet of . Thus it is easy to implement and numerically stable. Our main result establishes sharp rates of convergence of this algorithm providing an explicit link between the semimartingale characteristics of and its scale functions, not unlike the one-dimensional Itô diffusion setting, where scale functions are expressed in terms of certain integrals of the coefficients of the governing SDE.
46 pages, 4 figures