title: How can I add a kurtosis term in this Gaussian bell-curve formula?
body: <p>I am working with a formula to generate a curve in Python. I am tuning its parameters so that I can eyeball its matching over an underlying curve (already plotted).</p>
  
  <p>What I need is to adjust its Kurtosis, which currently is not a parameter of the formula:</p>
  
  <pre><code>def gaussian(x, peak_x, peak_y, sigma):
      return numpy.exp(-numpy.power(x - peak_x, 2.) / (2 * numpy.power(sigma, 2.))) * peak_y
  </code></pre>
  
  <p>I would need to expand the function to this signature:</p>
  
  <pre><code>def gaussian(x, peak_x, peak_y, sigma, KURTOSIS)
  </code></pre>
  
  <p>But I don't know where and how to change the formula.</p>