<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def rhalf_normal(tau, size=None):
 """ Random half-normal variates. """


    return abs(np.random.normal(0, np.sqrt(1 / tau), size))