M=[0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2];
T=[1.411, 1.424, 1.431, 1.428, 1.419, 1.430, 1.417, 1.422, 1.430, 1.432];
lM = log(M);
lT = log(T);
p = polyfit(lM,lT,1)
plot(lM,polyval(p,lM), 'r');
hold on;
plot(lM,lT,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('log(M)', 'FontSize', 16)
ylabel('log(T)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)