code: def draw(self):
          
          devs = normal(size=self._sig.shape[1])
          p = inner(self._sig, devs) + self._mu
          self._set_stochastics(p)
docs: N.draw()
  
          Sets all N's stochastics to random values drawn from
          the normal approximation to the posterior.