code: def make_arr(mu1, mu2, std1=1, std2=1, p=0.5):
      res = np.where(
          np.random.uniform(0, 1, n) > p,
          np.random.normal(mu1, std1, n),
          np.random.normal(mu2, std2, n),
      )
      return np.expand_dims(res, 1)
  
imports: from pathlib import Path
  import numpy
  import matplotlib.pylab
  import seaborn
  from sklego.naive_bayes import GaussianMixtureNB
label_desc: Use image node_id 1 for calling a global function with example usage: make_arr(mu1, mu2, std1, std2, p) and returns: np