plot made from octave with:

  T = [....];
  X = 1:1915;
  P = polyfit(X,T,1);
  M = P(1); C = P(2);
  Y = P(1)*X + C;
  plot(X,T,X,Y)
  ylim([-10,40])