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

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