Random numbers obeying beta distribution are utilized for simulation analyses of various stochastic models. However, the procedure for generating beta random numbers is limited depending on the values of the parameters. Therefore, in this paper, a unified procedure for generating beta random numbers independent of the values of the parameters is investigated using Wilson-Hilferty approximation and Stuart's proposition.