How to write a function to fit data to a sum of N Gaussian-like peaks without explicitly defining the expression for every possible N?

Question

I am trying to fit a progression of Gaussian peaks to a spectral lineshape. The progression is a summation of N evenly spaced Gaussian peaks. When code