%Model tx=qq1*(t0-tx)-qq2*(tx-tr)
tt=1:10000;                         % time
tr=20; kk=1/25;                     % RT and slope
t0=kk*tt+tr;                        % base temperature
qq1=10; qq2=6;
tx=(qq1*t0+qq2*tr)/(1+qq1-qq2);
tx0=(qq1*t0)/(1+qq1);
plot(tt,t0);
hold on
plot(tt,tx,'m');
plot(tt,tx0,'g');
title('Linear is not a good model!');
%
%plot(tt,tx,'g');
%qq2=2;
%tx=(qq1*t0+qq2*tr)/(1+qq1-qq2);
%plot(tt,tx,'m');
%qq2=3;
%tx=(qq1*t0+qq2*tr)/(1+qq1-qq2);
%plot(tt,tx,'m');
%figure
%plot(tt,t0);
%tx=(qq1*t0+qq2*tr)/(1+qq1-qq2);
%hold on
