SAS provides such a random function.

Weibull Distribution
x = RAND('WEIBULL',a,b)
Arguments
x
is an observation from the distribution with the following probability density function:

Range:x ≥ 0
a
is a numeric shape parameter.
Range:a > 0
b
is a numeric scale parameter.
Range:b > 0