By simulation, I observed that the efficient and best fit of the generated random variable is given by Algorithm 1(\ref{167436}). The shaded part is the density obtained from the random number generated and the curve in black covers original gamma distribution for that shape parameter. The vertical lines represent mean and median of the density obtained.