In this paper we obtain an explicit representation of the transition density of the one-dimensional skew Brownian motion with (a constant drift and) two semipermeable barriers. Moreover we propose a rejection method to simulate this density in an exact way.