DATA [...
];

year = DATA(:,1);
ICE = DATA(:,2);
plot(year,ICE)
X = [year, year.^0];
A = inv(X'*X)*X'*ICE

A =

   1.0e+02 *

  -0.000829642589118
   1.719082964351848