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

plot(year,ICE,'b.-',year,X*A,'r')
xlabel('Year');
ylabel('Ice Area');