

function X = weibull_RV(U,a, b)

X = (-log(U)/a)^(1/b)

%inverse of a Weibull