T=[0.24 0.62 1.88 11.9 29.5 84]';
r=[0.39 0.72 1.52 5.2 9.5 19]';
lnT=log(T);
lnR=log(r);
A=[lnR ones(size(r))];
qTt0=A\lnT;
q=qTt0(1);
T0=exp(qTt0(2));
rr=(0.3:0.01:19.5)';
TT=T0*rr.^q;
plot(r,T,'*',rr,TT)