/*
	r_exp <varname>

	generates expotentially distributed random variable with mean 2
*/

program define r_exp
version 2.1
	gen %_1 = -2*log(uniform())
end