Backward Stochastic Differential Equations (BSDEs) provide a systematic way to obtain Feynman-Kac formulas for linear as well as nonlinear partial differential equations (PDEs) of parabolic and elliptic type, and the numerical approximation of their solutions thus provide Monte-Carlo methods for PDEs. BSDEs are also used to describe the solution of path-dependent stochastic control problems, and they further arise in many areas of mathematical finance.

In this talk, I will discuss the numerical approximation of BSDEs when the nonlinear driver is not Lipschitz, but instead has polynomial growth and satisfies a monotonicity condition. The time-discretization is a crucial step, as it determines whether the full numerical scheme is stable or not. Unlike for Lipschitz driver, while the implicit Bouchard-Touzi-Zhang scheme is stable, the explicit one is not and explodes in general. I will then present a number of remedies that allow to recover a stable scheme, while benefiting from the reduced computational cost of an explicit scheme. I will also discuss the issue of numerical stability and the qualitative correctness which is enjoyed by both the implicit scheme and the modified explicit schemes. Finally, I will discuss the approximation of the expectations involved in the full numerical scheme, and their analysis when using a quasi-Monte Carlo method.