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