code: def draw_normal_initial(self):
          
          means, scale = self.get_means_and_scales_from_q()
          return np.random.normal(means,scale,size=[self.sims,means.shape[0]]).T
docs: Draw parameters from a mean-field normal family