code: def sample_ellipsoid(p0, covmat, size=1):
      """
      Produce an ellipsoid of walkers around an initial parameter value,
      according to a covariance matrix.
  
      :param p0: The initial parameter value.
      :param covmat:
          The covariance matrix.  Must be symmetric-positive definite or
          it will raise the exception numpy.linalg.LinAlgError
      :param size: The number of samples to produce.
  
      """
      return np.random.multivariate_normal(
          np.atleast_1d(p0), np.atleast_2d(covmat), size=size
      )
  
imports: import warnings
  from functools import wraps
  import numpy
label_desc: Use image node_id 4 for calling a global function with example usage: sample_ellipsoid(p0, covmat, size) and returns: np