x = [0;1;2;3;4;5;6]
y=[9.5;7.6;4.6;7.8;8.1;14.1;20.1]
A = [ x.^0, x.^1, x.^2]
c = pinv(A) * y
