daily_temp=xlsread('G:\phd\hydroinformatics\assgn2\Temp_Even_23.5_78.5.xls');
avg(1)=mean(daily_temp(1:60,4));
avg(2)=mean(daily_temp(61:119,4));
avg(3)=mean(daily_temp(120:178,4));
avg(4)=mean(daily_temp(179:237,4));
avg(5)=mean(daily_temp(238:297,4));
avg(6)=mean(daily_temp(298:356,4));
avg(7)=mean(daily_temp(357:415,4));
avg(8)=mean(daily_temp(416:474,4));
avg(9)=mean(daily_temp(475:534,4));
avg(10)=mean(daily_temp(535:593,4));
avg(11)=mean(daily_temp(594:652,4));
avg(12)=mean(daily_temp(653:711,4));
avg(13)=mean(daily_temp(712:771,4));
avg(14)=mean(daily_temp(772:830,4));
avg(15)=mean(daily_temp(831:889,4));
avg(16)=mean(daily_temp(890:948,4));
avg(17)=mean(daily_temp(949:1008,4));
avg(18)=mean(daily_temp(1009:1067,4));
avg(19)=mean(daily_temp(1068:1126,4));
avg(20)=mean(daily_temp(1127:1185,4));
avg(21)=mean(daily_temp(1186:1245,4));
avg(22)=mean(daily_temp(1246:1304,4));
avg(23)=mean(daily_temp(1305:1363,4));
avg(24)=mean(daily_temp(1364:1422,4));
avg(25)=mean(daily_temp(1423:1482,4));
avg(26)=mean(daily_temp(1482:1542,4));
avg(27)=mean(daily_temp(1543:1600,4));
avg(28)=mean(daily_temp(1601:1659,4));
avg(29)=mean(daily_temp(1660:1719,4));
avg(30)=mean(daily_temp(1720:1778,4));
avg(31)=mean(daily_temp(1779:1837,4));
avg(32)=mean(daily_temp(1838:1896,4));
avg(33)=mean(daily_temp(1897:1956,4));
avg(34)=mean(daily_temp(1957:2015,4));
Year=1980:2013;
scatter(Year,avg);
xlabel('Year'); ylabel('Annual Average Temperature');
linear_coeff = polyval(p,Year,1);
hold on
plot(Year,linear_coeff)
hold off
