  Novel multi-step predictor-corrector numerical schemes have been derived for
approximating decoupled forward-backward stochastic differential equations
(FBSDEs). The stability and high order rate of convergence of the schemes are
rigorously proved. We also present a sufficient and necessary condition for the
stability of the schemes. Numerical experiments are given to illustrate the
stability and convergence rates of the proposed methods.
