x = mu + sigma * qnorm( pnorm(a,mu,sigma) +
     runif(1)*(pnorm(b,mu,sigma) - pnorm(a,mu,sigma)) )