m=[0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2];
T=[0.25, 0.37, 0.46, 0.52, 0.57, 0.65, 0.69, 0.78, 0.76, 0.80];
lm = log10(m);
lT = log10(T);
p = polyfit(lm,lT,1)
plot(lm,polyval(p,lm), 'r');
hold on;
plot(lm,lT,'k*');
legend ('numerical fit', 'Data', 'location', 'southeast');
xlabel('log(m)', 'FontSize', 16)
ylabel('log(T)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)