code: def gauss_box_model(x, amplitude=1.0, mean=0.0, stddev=1.0, hpix=0.5):
      
      z = (x - mean) / stddev
      z2 = z + hpix / stddev
      z1 = z - hpix / stddev
      return amplitude * (norm.cdf(z2) - norm.cdf(z1))
docs: Integrate a Gaussian profile.