y = < 75.995   91.972  105.711  123.203   ...
     131.669  150.697  179.323  203.212>'
t = < 1900:10:1970 >'
t = (t - 1940*ones(t))/40;   <t y>
n = 8;  A(:,1) = ones(t);  for j = 2:n, A(:,j) = t .* A(:,j-1);
A
c = A\y
