%loads the file temperatures.dat
%each col is temperature in a different place
%computes average temperature, plots

load temperatures.dat

moon_temps = temperatures(:,1) ;
chicago_temps = temperatures(:,2) ;
seattle_temps = temperatures(:,3) ;

number_of_temps=size(temperatures,1);

summed_moon_temp=0 ;
summed_chicago_temp=0 ;
summed_seattle_temp=0 ;

for j=1:number_of_temps
    summed_moon_temp=summed_moon_temp+moon_temps(j) ;
    summed_chicago_temp=summed_chicago_temp+chicago_temps(j) ;
    summed_seattle_temp=summed_seattle_temp+seattle_temps(j) ;
end

average_moon_temp=summed_moon_temp/number_of_temps
average_chicago_temp=summed_chicago_temp/number_of_temps
average_seattle_temp=summed_seattle_temp/number_of_temps


figure
plot(1:number_of_temps,moon_temps,'-','Linewidth',2) ; hold on 
plot(1:number_of_temps,chicago_temps,':','Linewidth',2) 
plot(1:number_of_temps,seattle_temps,'--','Linewidth',2) 

plot(1:number_of_temps,average_moon_temp*ones(1,number_of_temps),'-','Linewidth',1)  
plot(1:number_of_temps,average_chicago_temp*ones(1,number_of_temps),'-','Linewidth',1)  
plot(1:number_of_temps,average_seattle_temp*ones(1,number_of_temps),'-','Linewidth',1)  

legend('moon temps','chicago temps','seattle temps')
xlabel('day','FontSize',12) 
ylabel('temperature','FontSize',12)


save chicago.dat average_chicago_temp -ascii