x = 1:0.1:2;
[~, n] = size(x);
y = x.^2 + 4*x + 7 + normrnd(0,2,1,n);
plot(x,y, '*')
poly_fit(x,y, 5)