h = [0.5, 0.48, 0.45, 0.44, 0.42, 0.39, 0.37, 0.36, 0.34, 0.32, 0.31, 0.29, 0.27, 0.26, 0.24, 0.23, 0.21, 0.2, 0.18, 0.17, 0.16, 0.15, 0.13, 0.12, 0.11];
t = [0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 210, 220, 230, 240];
p = polyfit(t, h,2)
plot(t,polyval(p,t));
hold on;
plot(t, h,'k*');
h = legend ('numerical fit', 'Data', 'location', 'northeast');
set (h, "fontsize", 16);
xlabel('t (s)', 'FontSize', 16)
ylabel('h (m)', 'FontSize', 16)
axis('tight')
set(gca,'FontSize',16)