def BetaPrime(alpha, beta, tag=None):
assert (<EOL>alpha > <NUM_LIT:0> and beta > <NUM_LIT:0><EOL>), '<STR_LIT>'<EOL>x = Beta(alpha, beta, tag)<EOL>return x / (<NUM_LIT:1> - x)<EOL>
A BetaPrime random variate

Parameters
----------
alpha : scalar
    The first shape parameter
beta : scalar
    The second shape parameter