code: def gaussian1(height, x0, y0, a, b, c):
      
  
      return lambda x, y: height * np.exp(-0.5* (a*(x-x0)**2 + b*(x-x0)*(y-y0) + c*(y-y0)**2))
docs: height - the amplitude of the gaussian
      x0, y0, - center of the gaussian
      a, b, c - ellipse parameters (coefficients in the quadratic form)