#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