We introduce a method for computation of hitting time distribution of an increasing Levy 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. (C) 2016 Elsevier B.V. All rights reserved.