code: def skewed_gaussian(
      x, amplitude=1.0, center=0.0, sigma=1.0, gamma=0.0
  ):
      """Return a Gaussian lineshape, skewed with error function.
  
      Equal to: gaussian(x, center, sigma)*(1+erf(beta*(x-center)))
  
      where ``beta = gamma/(sigma*sqrt(2))``
  
      with ``gamma < 0``: tail to low value of centroid
           ``gamma > 0``: tail to high value of centroid
  
      For more information, see:
      https://en.wikipedia.org/wiki/Skew_normal_distribution
  
      """
      asym = 1 + erf(gamma * (x - center) / max(tiny, (s2 * sigma)))
      return asym * gaussian(x, amplitude, center, sigma)
  
imports: from numpy import arctan, copysign, cos, exp, isclose, isnan, log, log1p, maximum, minimum, pi, real, sin, sqrt, where
  from scipy.special import betaln
  from scipy.special import erf, erfc
  from scipy.special import gamma
  from scipy.special import loggamma
  from scipy.special import wofz
label_desc: Use image node_id 19 for calling a global function with example usage: skewed_gaussian(x, amplitude, center, sigma, gamma) and returns: unknown