popt, pcov = curve_fit(gauss_function, x, y, p0=[1, mean, sigma])
plot(x, gauss_function(x, *popt))
plot(x, y, 'ok')
гауссовская подгонка с scipy.optimize.curve_fit в python с неправильными результатами