n = 1; count = 0; x=[0:1:6]; y=[2, 2.1592, 3.1697, 5.4332, 9.1411, 14.406, 21.303]; for i = 1:n:6 count = count+1; f1 = x(i:count+1); f2 = y(i:count+1); m = polyfit(f1,f2,1); slope(:,count) = m(1); intercept(:,count) = m(2); end slope intercept