points = [
1 1; #I
1 5;
3 1;
4 4;
-1 1; #II
-1 5;
-4 1;
-5 5;
-1 -1; #III
-1 -3;
-4 -1;
-4 -2;
1 -1; #IV
4 -1;
1 -5;
3 -5;
];
ht = size(points)(1);
wd = size(points)(2);
H1 = zeros(wd + 1, ht);
for i = 1: size(H1)(2)
H1(3, i) = 1;
endfor
for i = 1: (size(H1)(1)) - 1
for j = 1: size(H1)(2)
H1(i, j) = points(i, j)
endfor
endfor
#hold on;
#X = -10 : 0.01 : 10;
#Y = -10 : 0.01 : 10;
#plot(X, 0, 'k');
#plot(0, Y, 'k');
#for i = 1 : size(points)(1)
#if points(i, 1) > 0 && points(i, 2) > 0
#plot(points(i, 1), points(i, 2), 'r*');
#elseif points(i, 1) > 0 && points(i, 2) < 0
#plot(points(i, 1), points(i, 2), 'b*');
#elseif points(i, 1) < 0 && points(i, 2) < 0
#plot(points(i, 1), points(i, 2), 'g*');
#else
#plot(points(i, 1), points(i, 2), 'm*');
#endif
#endfor
#hold off;