x=[0 1 2 3 4 5];
y=[10 25 51 66 97 118];
plot(x,y,'x')
axis([-1 6 0 130])
xlabel('x')
ylabel('y')
title('Plot for Least Squares Example')
hold on

x2=linspace(0,5);
plot(x2, 22.0286.*x2 + 6.0952, 'k-')