Gaussian function:
y = a + b*exp(c*s^2)
| int n | number of data points. |
| double* x | x array (at least order+1 values). |
| double* y | y array (at least order+1 values). |
| double* coeff | array in which 3 coefficients are returned |
| (if NULL, | no coefficients returned). |