In this work, we study the numerical approximation of a class singular fully coupled forward backward stochastic differential equations. These equations have degenerate component and non-smooth terminal condition. They are used, for example, in modelling carbon market [1] linked to scalar conservation law perturbed by diffusion. Classical FBSDEs methods fail capture correct entropy solution ass...