f = [0.20, 0.39, 0.59, 0.78, 0.98, 1.18, 1.37, 1.57, 1.77, 1.96];
x = [0.017, 0.034, 0.051, 0.068, 0.085, 0.103, 0.118, 0.135, 0.154, 0.171];
lf = log(f);
lx = log(x);
p = polyfit(lx, lf,1)
plot(lx,polyval(p,lx));
hold on;
plot(lx,lf,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('log(x)', 'FontSize', 16)
ylabel('log(F)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)