t=[0.040, 0.080, 0.120, 0.160, 0.200, 0.240, 0.280];
y=[0.009, 0.036, 0.085, 0.158, 0.251, 0.365, 0.501];
lt = log10(t);
ly = log10(y);
p = polyfit(lt,ly,1)
plot(lt,polyval(p,lt), 'r');
hold on;
plot(lt,ly,'k*');
h = legend ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
xlabel('log(t)', 'FontSize', 16)
ylabel('log(y)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)