L=[0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65];
T=[0.903, 1.011, 1.095, 1.193, 1.262, 1.338, 1.412, 1.488, 1.562, 1.632];
lL = log(L);
lT = log(T);
p = polyfit(lL,lT,1)
plot(lL,polyval(p,lL), 'r');
hold on;
plot(lL,lT,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('log(L)', 'FontSize', 16)
ylabel('log(T)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)