load 'MLOTT5y'
M = MLOTT5y;

mat1 = M (1:31,:);

years1 = mat1 (:,1);
T1 = mat1 (:,2);

bb1 = robustfit (years1,T1);



load 'MLOTT5y'
M = MLOTT5y;

mat2 = M (2:32,:);

years2 = mat2 (:,1);
T2 = mat2 (:,2);

bb2 = robustfit (years2,T2);


load 'MLOTT5y'
M = MLOTT5y;

mat3 = M (3:33,:);

years3 = mat3 (:,1);
T3 = mat3 (:,2);

bb3 = robustfit (years3,T3);


load 'MLOTT5y'
M = MLOTT5y;

mat4 = M (4:34,:);

years4 = mat4 (:,1);
T4 = mat4 (:,2);

bb4 = robustfit (years4,T4);



load 'MLOTT5y'
M = MLOTT5y;

mat5 = M (5:35,:);

years5 = mat5 (:,1);
T5 = mat5 (:,2);

bb5 = robustfit (years5,T5);


load 'MLOTT5y'
M = MLOTT5y;

mat6 = M (6:36,:);

years6 = mat6 (:,1);
T6 = mat6 (:,2);

bb6 = robustfit (years6,T6);


load 'MLOTT5y'
M = MLOTT5y;

mat7 = M (7:37,:);

years7 = mat7 (:,1);
T7 = mat7 (:,2);

bb7 = robustfit (years7,T7);



load 'MLOTT5y'
M = MLOTT5y;

mat8 = M (8:38,:);

years8 = mat8 (:,1);
T8 = mat8 (:,2);

bb8 = robustfit (years8,T8);


load 'MLOTT5y'
M = MLOTT5y;

mat9 = M (9:39,:);

years9 = mat9 (:,1);
T9 = mat9 (:,2);

bb9 = robustfit (years9,T9);


trendchangesplot = figure;
hold on
trenchangesplot = plot (years1,years1*bb1(2)+bb1(1),'blue','Linewidth',2);
hold on
grid on
trenchangesplot = plot (years2,years2*bb2(2)+bb2(1),'blue','Linewidth',2);
hold on
grid on
trenchangesplot = plot (years3,years3*bb3(2)+bb3(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years4,years4*bb4(2)+bb4(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years5,years5*bb5(2)+bb5(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years6,years6*bb6(2)+bb6(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years7,years7*bb7(2)+bb7(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years8,years8*bb8(2)+bb8(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years9,years9*bb9(2)+bb9(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('temperature (Celsius)','Fontsize',12);


load 'MLOTT5y'
MM = MLOTT5y;

load 'bb'
bb = bb;

years = MM (:,1);
T = MM (:,2);

trendchengesmean = figure;
hold on
trendchangesmean = plot (years, T, 'cyan','Linewidth',2);
hold on
trenschangesmean = plot (years, years*bb(2) + bb(1),'magenta','Linewidth',2);
hold on
grid on
trenchangesmean = plot (years1,years1*bb1(2)+bb1(1),'blue','Linewidth',2);
hold on
grid on
trenchangesmean = plot (years2,years2*bb2(2)+bb2(1),'blue','Linewidth',2);
grid on
trenchangesmean = plot (years3,years3*bb3(2)+bb3(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years4,years4*bb4(2)+bb4(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years5,years5*bb5(2)+bb5(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years6,years6*bb6(2)+bb6(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years7,years7*bb7(2)+bb7(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years8,years8*bb8(2)+bb8(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years9,years9*bb9(2)+bb9(1),'blue','Linewidth',2);
hold on
%trenchangesmean = plot (years10,years10*bb10(2)+bb10(1),'blue','Linewidth',2);

hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('temperature (Celsius)','Fontsize',12);
hold on
axis ([1967 2012 6 8]);
hold on
legend ('5 yr running average','0.04 yrs - 71.31 (C / yr)','30 yrs trends','location','Southeast');