We propose a numerical algorithm for backward stochastic differential equations based on time discretization and trigonometric wavelets. This method combines the effectiveness of Fourier-based methods and the simplicity of a wavelet-based formula, resulting in an algorithm that is both accurate and easy to implement. Furthermore, we mitigate the problem of errors near the computation boundaries by means of an antireflective boundary technique, giving an improved approximation. We test our algorithm with different numerical experiments.In this article, we shall adopt the quick variant of the SWIFT formula proposed in [13].The reason behind this is left for the error section, but this gives an approximation for E x p [v(t p+1 , X ∆ t p+1 )], which we see in the discrete-time approximation of FBSDE,