Numerical computation of hitting time distributions of increasing Lévy processes

We introduce a method for computation of hitting time distribution of an increasing Lévy process using the inverse Fourier transform and the Hilbert transform under the assumption that the characteristic function of the process is given. Its efficiency is demonstrated by the Kolmogorov-Smirnov test.