#define A 1
function  gauss_from_def1
loop
x[i] = p[0] * exp(-(t[i] - p[1])*(t[i] - p[1])*p[2]);
end function
/* comment line */
jacobian jacgauss_from_def1
double arg, expf;
loop
arg = t[i] - p[1];
expf = exp(-arg*arg*p[2]);
d0 = expf;
d1[i] = p[0]*2*arg*p[2]*expf;
d2 = p[0]*(-arg*arg)*expf;
end jacobian