time_even=0:2:22;
time_odd=1:2:23;
temp_even= [9.2, 9.0, 6.8, 5.9, 7.2, 9.2, 10.9, 10.6, 10.3, 8.2, 8.8, 7.5];
temp_odd=  [9.6, 7.7, 6.0, 4.3, 8.9, 10.2, 10.6, 10.7, 8.6, 8.0, 7.1, 7.5]; 
pol=polyfit(time_even,temp_even,1);
temp_int_odd=polyval(pol,time_odd);
MSE=sum((temp_int_odd-temp_odd).^2)/length(temp_odd);
temp=zeros(1,24);
temp((0:2:22)+1)=temp_even;
temp((1:2:23)+1)=temp_odd;
temp_int=zeros(1,24);
temp_int((0:2:22)+1)=temp_even;
temp_int((1:2:23)+1)=temp_int_odd;
plot(0:23,temp,'b-',0:23,temp_int,'r',0:23,polyval(pol,0:23),'k')
