code: def gauss(x,m_y,sigma):
      '''returns the gaussian with mean m_y and std. dev. sigma,
  calculated at the points of x.'''
  
      e_y = [np.exp((1.0/(2*float(sigma)**2)*-(n-m_y)**2)) for n in np.array(x)]
      y = [1.0/(float(sigma) * np.sqrt(2 * np.pi)) * e for e in e_y]
      
      return np.array(yf gauss(x,m_y,sigma):
      '''returns the gaussian with mean m_y and std. dev. sigma,
  calculated at the points of x.'''
  
      e_y = [np.exp((1.0/(2*float(sigma)**2)*-(n-m_y)**2)) for n in np.array(x)]
      y = [1.0/(float(sigma) * np.sqrt(2 * np.pi)) * e for e in e_y]
      
      return np.array(y)
docs: returns the gaussian with mean m_y and std. dev. sigma,
  calculated at the points of x.