code: def gaussian(x, mu, sig):
      a = (x - mu) ** 2 / (2 * (sig**2))
      b = 1 / (sig * (gs.sqrt(2 * gs.pi)))
      f = b * gs.exp(-a)
      l2_norm = gs.sqrt(gs.trapz(f**2, x))
      f_sinf = f / l2_norm
  
      return gs.array([f_sinf])
  
file_path: geomstats/tests/tests_geomstats/test_geometry/data/functions.py
imports: import geomstats.backend
  from geomstats.test.data import TestData
  from .manifold import ManifoldTestData
  from .riemannian_metric import RiemannianMetricTestData
label_desc: Use image node_id 1 for calling a global function with example usage: gaussian(x, mu, sig) and returns: gs