This paper presents simulation formulae of two-sided truncated normal rando
m variables using a completion distribution and its two corresponding condi
tionals generated via a Gibbs sampler. This procedure extends formulae give
n by Robert and Casella for the one-sided case.