AbstractIn this work, in order to obtain higher-order schemes for solving forward backward stochastic differential equations, we propose a new multi-step scheme by adopting the high-order multi-step method in Zhao et al. (SIAM J. Sci. Comput., 36(4): A1731-A1751, 2014) with the combination technique. Two reference ordinary differential equations containing the conditional expectations and their derivatives are derived from the backward component. These derivatives are approximated by using the finite difference methods with multi-step combinations. The resulting scheme is a semi-discretization in the temporal direction involving conditional expectations, which are solved by using the Gaussian quadrature rules and polynomial interpolations on the spatial grids. Our new proposed multi-step scheme allows for higher convergence rate up to ninth order, and are more efficient. Finally, we provide a numerical illustration of the convergence of the proposed method.