L = [0.93, 0.99, 1.06, 1.03, 1, 0.9, 0.95, 0.96, 1.05, 0.95, 0.9922, 1.034, 0.98, 0.94, 0.97, 0.93, 0.94, 0.79, 0.94, 0.89,	0.91, 0.91, 0.84, 0.84, 0.91, 0.91, 0.95, 0.99, 0.96, 0.91, 0.91, 0.86, 0.96, 0.94, 0.90, 0.96, 0.89, 1.04, 0.99, 0.99];
d = [0.71, 0.71, 0.77, 0.77, 0.67, 0.73, 0.67, 0.77, 0.81, 0.71, 0.67, 0.67, 0.59, 0.56, 0.63, 0.61, 0.69, 0.73, 0.76, 0.61, 0.76, 0.67, 0.60, 0.66, 0.76, 0.68, 0.64, 0.76, 0.74, 0.68, 0.66, 0.62, 0.71, 0.73, 0.68, 0.64, 0.69, 0.74, 0.68, 0.78];
v = [1.19, 1.33, 1.27, 1.18, 1.30, 1.44, 1.24, 1.37, 1.40, 1.33, 1.16, 1.18, 1.12, 1.07, 1.15, 1.12, 1.36, 1.51, 1.36, 1.03, 1.70, 1.40, 1.07, 1.14, 1.61, 1.33, 1.20, 1.36, 1.50, 1.28, 1.26, 1.20, 1.26, 1.22, 1.40, 1.05, 1.21, 1.28, 1.40, 1.38];
lL = log(L);
lv = log(v);
p = polyfit(lL,lv,1)
plot(lL,polyval(p,lL), 'r');
hold on;
plot(lL,lv,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('log(L)', 'FontSize', 16)
ylabel('log(v walk)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)