R function: rgamma

rgamma(5, shape=4, rate=0.333)
# [1] 18.416880  5.357108  9.588073 15.704966 13.772501