Abstract : We provide two numerical schemes for the sample paths simulation of skew diffusions with piecewise constant coefficients. Both of these schemes use the resolvent kernel. Precisely, the first scheme use only the resolvent kernel while the second use it as less as possible. Moreover, in the second, some sampled positions are exact.