code: def loglike_multigauss(x):
      """two-peaked gaussian problem"""
      a = -0.5 * (((x - 0.4) / 0.01) ** 2).sum()
      b = -0.5 * (((x - 0.6) / 0.01) ** 2).sum()
      return np.logaddexp(a, b)
  
imports: import numpy
  from math import gamma, pi, exp
label_desc: Use image node_id 15 for calling a global function with example usage: loglike_multigauss(x) and returns: np