For a regime-switching model with a finite number of regimes and double phase-type jumps, Jiang and Pistorius (2008) derived matrix equations with real parameters for the Wiener–Hopf factorization. The Laplace transform of the first passage time distribution is expressed in terms of the solution of the matrix equations. In this paper we provide an iterative algorithm for solving the matrix equations of Jiang and Pistorius (2008) with complex parameters. This makes it possible to obtain numeric values of the Laplace transform with complex parameters for the first passage time distribution. The Laplace transform with complex parameters can be inverted by numerical inversion algorithms such as the Euler method. As an application, we compute the prices of defaultable bonds under a structural model with regime switching and double phase-type jumps.
