Fitting multiple peaks procedure

#pragma rtGlobals=1     // Use modern global access method.
#include  <Multi-peak fitting 1.4>

Function Gauss3(w,x) : FitFunc
    Wave w
    Variable x

    //CurveFitDialog/ These comments were created by the Curve Fitting dialog. Altering them will
    //CurveFitDialog/ make the function less convenient to work with in the Curve Fitting dialog.
    //CurveFitDialog/ Equation:
    //CurveFitDialog/ f(x) = y0+A1*exp(((x-X1)/B1)^2)+A2exp(((x-X2)/B2)^2)+A3*exp(((x-X3)/B3)^2)
    //CurveFitDialog/ End of Equation
    //CurveFitDialog/ Independent Variables 1
    //CurveFitDialog/ x
    //CurveFitDialog/ Coefficients 7
    //CurveFitDialog/ w[0] = y0
    //CurveFitDialog/ w[1] = A1
    //CurveFitDialog/ w[2] = B1
    //CurveFitDialog/ w[3] = X1
    //CurveFitDialog/ w[4] = A2
    //CurveFitDialog/ w[5] = B2
    //CurveFitDialog/ w[6] = X2
    //CurveFitDialog/ w[7] = A3
    //CurveFitDialog/ w[8] = B3
    //CurveFitDialog/ w[9] = X3

    return w[0]+w[1]*exp(((x-w[3])/w[2])^2)+w[4]*exp(((x-w[6])/w[5])^2)+w[7]*exp(((x-w[9])/w[8])^2)

End

Function Run()

make/O/N = (10) coef

coef={25147,4683.1,144.85,40.332,983.84,261.7,59.678,1258.7,77.517,130.99}

FuncFit Gauss3 coef test /X = Xwave


End