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];
time_even=0:2:22;
time_odd=1:2:23;
T=polyfit(time_even,temp_even,1);
f=polyval(T,time_odd);
MSE=sum((f-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_e=zeros(1,24);
temp_e((0:2:22)+1)=temp_even;
temp_e((1:2:23)+1)=f;
plot(0:23,temp,'b',0:23,temp_e,'g',0:23,polyval(T,0:23))