D = readmatrix('Example.xlsx');
x = D(:,1);
y = D(:,2);
p = polyfit(x, y, 3);
yfit = polyval(p, x);
dp = polyder(p);
dyfit = polyval(dp,x(1));
bint = (yfit(1) - dyfit*x(1));
xinitslope = [x(1)-x(75), x(1)+x(75)];
yinitslope = dyfit*xinitslope + bint;
figure
plot(x, y)
hold on
plot(x, yfit, 'LineWidth',1.5)
plot(xinitslope, yinitslope, 'LineWidth',3)
hold off
grid
legend('Data', '3\circ Polynomial Fit', 'Initial Slope', 'Location','N')
text(xinitslope(1), yinitslope(1), sprintf('\\leftarrow Initial Slope = %7.2f',dyfit), 'HorizontalAlignment','left', 'VerticalAlignment','middle')
Best Answer