x = linspace(0,1,1001);

for i=1:10000
a = 1 + 0.1*randn();
b = 1 + 0.1*randn();
c = 1 + 0.1*randn();
y = (a+b.*x+c.*x.^2)';
A = [x.^0; x; x.^2]';
lse(i,:) = A \ y;
end