Подсчитать количество точек множества, находящихся внутри этого треугольника

CLS
RANDOMIZE TIMER
SCREEN 12
10 INPUT "Kilkist tochok n="; n
IF n < 0 THEN 10
IF n <> INT(n) THEN 10
DIM x(1 TO n)
DIM y(1 TO n)
Xmin = 0
Xmax = 1000
Ymin = 0
Ymax = 1000
VIEW (20, 20)-(450, 450), , 2
WINDOW (Xmin, Ymin)-(Xmax, Ymax)
FOR i = 1 TO n
    x(i) = INT(Xmin + RND * (Xmax - Xmin))
    y(i) = INT(Ymin + RND * (Ymax - Ymin))
    PSET (x(i), y(i)), 3
NEXT
 
smax = 0
FOR i = 1 TO n
    FOR j = 1 TO n
        FOR k = 1 TO n
            s = ABS((x(i) - x(j)) * (y(i) + y(j)) + (x(j) - x(k)) * (y(j) + y(k)) + (x(k) - x(i)) * (y(k) + y(i))) / 2
            IF smax < s THEN
                smax = s
                X1 = x(i): Y1 = y(i)
                X2 = x(j): Y2 = y(j)
                X3 = x(k): Y3 = y(k)
            END IF
        NEXT k
    NEXT j
NEXT i
PRINT "Smax="; smax / 2
 
LINE (X1, Y1)-(X2, Y2), 12
LINE (X2, Y2)-(X3, Y3), 14
LINE (X3, Y3)-(X1, Y1), 10
 
k = 0
FOR i = 1 TO n
    k = k + s(X1, Y1, X2, Y2, X3, Y3, x(i), y(i))
NEXT i
 
PRINT "KOLICHESTVO TOCHEK ="; k - 2
 
END

CLS
RANDOMIZE TIMER
SCREEN 12
10 INPUT "Kilkist tochok n="; n
IF n < 0 THEN 10
IF n <> INT(n) THEN 10
DIM x(1 TO n)
DIM y(1 TO n)
Xmin = 0
Xmax = 1000
Ymin = 0
Ymax = 1000
VIEW (20, 20)-(450, 450), , 2
WINDOW (Xmin, Ymin)-(Xmax, Ymax)
FOR i = 1 TO n
    x(i) = INT(Xmin + RND * (Xmax - Xmin))
    y(i) = INT(Ymin + RND * (Ymax - Ymin))
    PSET (x(i), y(i)), 3
NEXT
 
smax = 0
a = 0
FOR i = 1 TO n
    FOR j = 1 TO n
        FOR k = 1 TO n
            FOR h = 1 TO n
 
                a1 = SQR((x(i) - x(j)) ^ 2 + (y(i) - y(j)) ^ 2)
                b1 = SQR((x(j) - x(k)) ^ 2 + (y(j) - y(k)) ^ 2)
                c1 = SQR((x(k) - x(i)) ^ 2 + (y(k) - y(i)) ^ 2)
                p12 = a1 + b1 + c1
                p1 = p12 / 2
                s1 = SQR(p1 * (p1 - a1) * (p1 - b1) * (p1 - c1))
                IF smax < s1 THEN
                    smax = s1
                    X1 = x(i): Y1 = y(i)
                    X2 = x(j): Y2 = y(j)
                    X3 = x(k): Y3 = y(k)
 
 
                    a2 = SQR((x(i) - x(j)) ^ 2 + (y(i) - y(j)) ^ 2)
                    b2 = SQR((x(j) - x(h)) ^ 2 + (y(j) - y(h)) ^ 2)
                    c2 = SQR((x(h) - x(i)) ^ 2 + (y(h) - y(i)) ^ 2)
                    p22 = a2 + b2 + c2
                    p2 = p22 / 2
                    s2 = SQR(p2 * (p2 - a2) * (p2 - b2) * (p2 - c2))
 
                    a3 = SQR((x(i) - x(h)) ^ 2 + (y(i) - y(h)) ^ 2)
                    b3 = SQR((x(h) - x(k)) ^ 2 + (y(h) - y(k)) ^ 2)
                    c3 = SQR((x(k) - x(i)) ^ 2 + (y(k) - y(i)) ^ 2)
                    p32 = a3 + b3 + c3
                    p3 = p32 / 2
                    s3 = SQR(p3 * (p3 - a3) * (p3 - b3) * (p3 - c3))
 
                    a4 = SQR((x(h) - x(j)) ^ 2 + (y(h) - y(j)) ^ 2)
                    b4 = SQR((x(j) - x(k)) ^ 2 + (y(j) - y(k)) ^ 2)
                    c4 = SQR((x(k) - x(h)) ^ 2 + (y(k) - y(h)) ^ 2)
                    p42 = a4 + b4 + c4
                    p4 = p42 / 2
                    s4 = SQR(p4 * (p4 - a4) * (p4 - b4) * (p4 - c4))
 
                    X4 = x(h): Y4 = y(h)
                    IF ABS(s2 + s3 + s4) = s1 THEN a = a + 1
 
                END IF
            NEXT h
        NEXT k
    NEXT j
NEXT i
 
PRINT "KILKIS TOCHOK ="; a
 
PRINT "Smax="; smax / 2
 
LINE (X1, Y1)-(X2, Y2), 12
LINE (X2, Y2)-(X3, Y3), 14
LINE (X3, Y3)-(X1, Y1), 10
 
END