x = [-2, -1, 0, 1, 2];
polyfit (x, x.^2+x+1)
