Second order polynomial regression function
[f,df] = regpoly2(S)
m*n matrix with design sites
[1 S S(:,1)*S S(:,2)S(:,2:n) ... S(:,n)^2]
Jacobian at the first point (first row in S)
ycollet@freesurf.fr