load 'MLOprec5y'
M = MLOprec5y;

mat1 = M (1:21,:);

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

bb1 = robustfit (years1,T1);



load 'MLOprec5y'
M = MLOprec5y;

mat2 = M (2:22,:);

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

bb2 = robustfit (years2,T2);


load 'MLOprec5y'
M = MLOprec5y;

mat3 = M (3:23,:);

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

bb3 = robustfit (years3,T3);


load 'MLOprec5y'
M = MLOprec5y;

mat4 = M (4:24,:);

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

bb4 = robustfit (years4,T4);



load 'MLOprec5y'
M = MLOprec5y;

mat5 = M (5:25,:);

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

bb5 = robustfit (years5,T5);


load 'MLOprec5y'
M = MLOprec5y;

mat6 = M (6:26,:);

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

bb6 = robustfit (years6,T6);


load 'MLOprec5y'
M = MLOprec5y;

mat7 = M (7:27,:);

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

bb7 = robustfit (years7,T7);



load 'MLOprec5y'
M = MLOprec5y;

mat8 = M (8:28,:);

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

bb8 = robustfit (years8,T8);


load 'MLOprec5y'
M = MLOprec5y;

mat9 = M (9:29,:);

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

bb9 = robustfit (years9,T9);


load 'MLOprec5y'
M = MLOprec5y;

mat10 = M (10:30,:);

years10 = mat10 (:,1);
T10 = mat10 (:,2);

bb10 = robustfit (years10,T10);



load 'MLOprec5y'
M = MLOprec5y;

mat11 = M (11:31,:);

years11 = mat11 (:,1);
T11 = mat11 (:,2);

bb11 = robustfit (years11,T11);



load 'MLOprec5y'
M = MLOprec5y;

mat12 = M (12:32,:);

years12 = mat12 (:,1);
T12 = mat12 (:,2);

bb12 = robustfit (years12,T12);



load 'MLOprec5y'
M = MLOprec5y;

mat13 = M (13:33,:);

years13 = mat13 (:,1);
T13 = mat13 (:,2);

bb13 = robustfit (years13,T13);


load 'MLOprec5y'
M = MLOprec5y;

mat14 = M (14:34,:);

years14 = mat14 (:,1);
T14 = mat14 (:,2);

bb14 = robustfit (years14,T14);


load 'MLOprec5y'
M = MLOprec5y;

mat15 = M (15:35,:);

years15 = mat15 (:,1);
T15 = mat15 (:,2);

bb15 = robustfit (years15,T15);


load 'MLOprec5y'
M = MLOprec5y;

mat16 = M (16:36,:);

years16 = mat16 (:,1);
T16 = mat16 (:,2);

bb16 = robustfit (years16,T16);


load 'MLOprec5y'
M = MLOprec5y;

mat17 = M (17:37,:);

years17 = mat17 (:,1);
T17 = mat17 (:,2);

bb17 = robustfit (years17,T17);



load 'MLOprec5y'
M = MLOprec5y;

mat18 = M (18:38,:);

years18 = mat18 (:,1);
T18 = mat18 (:,2);

bb18 = robustfit (years18,T18);



load 'MLOprec5y'
M = MLOprec5y;

mat19 = M (19:39,:);

years19 = mat19 (:,1);
T19 = mat19 (:,2);

bb19 = robustfit (years19,T19);


load 'MLOprec5y'
M = MLOprec5y;

mat20 = M (20:40,:);

years20 = mat20 (:,1);
T20 = mat20 (:,2);

bb20 = robustfit (years20,T20);



load 'MLOprec5y'
M = MLOprec5y;

mat21 = M (21:41,:);

years21 = mat21 (:,1);
T21 = mat21 (:,2);

bb21 = robustfit (years21,T21);



load 'MLOprec5y'
M = MLOprec5y;

mat22 = M (22:42,:);

years22 = mat22 (:,1);
T22 = mat22 (:,2);

bb22 = robustfit (years22,T22);


load 'MLOprec5y'
M = MLOprec5y;

mat23 = M (23:43,:);

years23 = mat23 (:,1);
T23 = mat23 (:,2);

bb23 = robustfit (years23,T23);


load 'MLOprec5y'
M = MLOprec5y;

mat24 = M (24:44,:);

years24 = mat24 (:,1);
T24 = mat24 (:,2);

bb24 = robustfit (years24,T24);

load 'MLOprec5y'
M = MLOprec5y;

mat25 = M (25:45,:);

years25 = mat25 (:,1);
T25 = mat25 (:,2);

bb25 = robustfit (years25,T25);

load 'MLOprec5y'
M = MLOprec5y;

mat26 = M (26:46,:);

years26 = mat26 (:,1);
T26 = mat26 (:,2);

bb26 = robustfit (years26,T26);

load 'MLOprec5y'
M = MLOprec5y;

mat27 = M (27:47,:);

years27 = mat27 (:,1);
T27 = mat27 (:,2);

bb27 = robustfit (years27,T27);


load 'MLOprec5y'
M = MLOprec5y;

mat28 = M (28:48,:);

years28 = mat28 (:,1);
T28 = mat28 (:,2);

bb28 = robustfit (years28,T28);


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
trenchangesplot = plot (years10,years10*bb10(2)+bb10(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years11,years11*bb11(2)+bb11(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years12,years12*bb12(2)+bb12(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years13,years13*bb13(2)+bb13(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years14,years14*bb14(2)+bb14(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years15,years15*bb15(2)+bb15(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years16,years16*bb16(2)+bb16(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years17,years17*bb17(2)+bb17(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years18,years18*bb18(2)+bb18(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years19,years19*bb19(2)+bb19(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years20,years20*bb20(2)+bb20(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years21,years21*bb21(2)+bb21(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years22,years22*bb22(2)+bb22(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years23,years23*bb23(2)+bb23(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years24,years24*bb24(2)+bb24(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years25,years25*bb25(2)+bb25(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years26,years26*bb26(2)+bb26(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years27,years27*bb27(2)+bb27(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years28,years28*bb28(2)+bb28(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('precipitation (mm)','Fontsize',12);


load 'MLOprec5y'
MM = MLOprec5y;

load 'bbMLOprec'
bb = bbMLOprec;

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);
hold on
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
trenchangesmean = plot (years11,years11*bb11(2)+bb11(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years12,years12*bb12(2)+bb12(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years13,years13*bb13(2)+bb13(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years14,years14*bb14(2)+bb14(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years15,years15*bb15(2)+bb15(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years16,years16*bb16(2)+bb16(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years17,years17*bb17(2)+bb17(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years18,years18*bb18(2)+bb18(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years19,years19*bb19(2)+bb19(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years20,years20*bb20(2)+bb20(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years21,years21*bb21(2)+bb21(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years22,years22*bb22(2)+bb22(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years23,years23*bb23(2)+bb23(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years24,years24*bb24(2)+bb24(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years25,years25*bb25(2)+bb25(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years26,years26*bb26(2)+bb26(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years27,years27*bb27(2)+bb27(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years28,years28*bb28(2)+bb28(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('precipitation (mm)','Fontsize',12);
hold on
axis ([1957 2012 0 1000]);
hold on
legend ('5 yr running average','3.45 yrs + 7.29 e3 (mm / yr)','20 yrs trends');