code: def draw_variable_local(self, size): 
          
          return ss.norm.rvs(loc=self.mu0, scale=self.sigma0, size=size)
docs: Simulate from the Normal distribution using instance values
  
          Parameters
          ----------
          size : int
              How many simulations to perform
  
          Returns
          ----------
          np.ndarray of Normal random variable