#p1(1, 1) p2(2, 2) p3(3, 5) p4(10, 21) p5(42, 22)
points = [1, 1;
2, 2;
3, 5;
10, 21;
42, 22];
disp(points);
ht = size(points)(1);
wd = size(points)(2);
for i = 1:ht
X(i, 1) = 1;
X(i, 2) = points(i, 1);
Y(i) = points(i, 2);
endfor
W = zeros(2, 1);
disp(X);
disp(Y);
disp(W);
hold on
for i=1:ht
plot(points(i, 1), points(i, 2), 'ro');
endfor
plot(points(1:5, 1), points(1:5, 2), 'b');