  In this paper we propose a generalized numerical scheme for backward
stochastic differential equations(BSDEs). The scheme is based on approximation
of derivatives via Lagrange interpolation. By changing the distribution of
sample points used for interpolation, one can get various numerical schemes
with different stability and convergence order. We present a condition for the
distribution of sample points to guarantee the convergence of the scheme.
