We present a new algorithm to discretize a decoupled forward‒backward stochastic differential equation driven by a pure jump Lévy process (FBSDEL for short). The method consists of two steps. In the first step we approximate the FBSDEL by a forward‒backward stochastic differential equation driven by a Brownian motion and Poisson process (FBSDEBP for short), in which we replace the small jumps by a Brownian motion. Then, we prove the convergence of the approximation when the size of small jumps ε goes to 0. In the second step we obtain the Lp-Hölder continuity of the solution of the FBSDEBP and we construct two numerical schemes for this FBSDEBP. Based on the Lp-Hölder estimate, we prove the convergence of the scheme when the number of time steps n goes to ∞. Combining these two steps leads to the proof of the convergence of numerical schemes to the solution of FBSDEs driven by pure jump Lévy processes.
