code: def make_data(n=10, mask_shape=(10, 10, 10), axis=0, r=3, signal=5):
      """
      Generate Gaussian noise in a cubic volume
      + cubic activations
      """
      mask = np.zeros(mask_shape, int)
      XYZ = np.array(np.where(mask == 0))
      p = XYZ.shape[1]
      data = np.random.randn(n, p)
      I = np.where(
          np.square(XYZ - XYZ.max(axis=1).reshape(-1, 1) / 2).sum(
              axis=0
          )
          <= r**2
      )[0]
      data[:, I] += signal
      vardata = np.random.randn(n, p) ** 2
      if axis == 1:
          data = data.T
          vardata = vardata.T
      return data, vardata, XYZ
  
label_desc: Use image node_id 1 for calling a global function with example usage: make_data(n, mask_shape, axis, r, signal) and returns: data,  vardata,  XYZ