L=[0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5];
M=[0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2];
T=[0.903, 1.011, 1.095, 1.193, 1.262, 1.338, 1.412, 1.488, 1.562, 1.632, 1.411, 1.424, 1.431, 1.428, 1.419, 1.430, 1.417, 1.422, 1.430, 1.432];
m = 0.49777;
n = 0.0034532;
x = L;
for l = 1:length(L)
	x(l) = (L(l)^m)*(M(l)^n);
end
p = polyfit(x,T,1)
plot(x,polyval(p,x), 'r');
hold on;
plot(x,T,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('L^m M^n', 'FontSize', 16)
ylabel('T', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)