code: def volume_multigauss(loglike, ndim):
      """compute volume enclosed at loglike threshold"""
      sqr_radius = -2 * loglike
      radius = sqr_radius**0.5 * 0.01
      if radius >= 0.5:
          # the volume is still touching the unit cube
          return np.nan
      if radius >= (0.2**2 * ndim) ** 0.5:
          # the two peaks are still touching each other
          return np.nan
  
      # compute volume of a n-sphere
      return nsphere_volume(radius, ndim)
  
imports: import numpy
  from math import gamma, pi, exp
label_desc: Use image node_id 17 for calling a global function with example usage: volume_multigauss(loglike, ndim) and returns: nsphere_volume, np, np