We study the weak approximation error of a skew diffusion with bounded measurable drift and Hölder diffusion coefficient by an Euler-type scheme, which consists of iteratively simulating skew Brownian motions with constant drift. We first establish two sided Gaussian bounds for the density of this approximation scheme. Then, a bound for the difference between the densities of the skew diffusion and its Euler approximation is obtained. Notably, the weak approximation error is shown to be of order $h^{\eta/2}$, where $h$ is the time step of the scheme, $\eta$ being the Hölder exponent of the diffusion coefficient.
Add the publication’s full text or supplementary notes here. You can use rich formatting such as including code, math, and images.