Составить программу для расчета траектории движения целей

SCREEN

PSET (0 + shiftX, 0 + shiftY), 15 ' centr koordinat gor oloskosti

SCREEN 
PSET (0 + shiftX, 0 + shiftY), 15 ' centr koordinat gor oloskosti  ...

OPEN "data.txt" FOR INPUT AS #1
i = 0 ...

OPEN """c:\my path\data.txt""" FOR INPUT AS #1

DIM azim(i) AS STRING, Vgor(i) AS STRING, Vver(i) AS STRING, kurs(i) AS STRING, vis(i) AS STRING, daln(i) AS STRING
DIM N AS INTEGER, Dr AS INTEGER
 
CLS
 
OPEN "data.txt" FOR INPUT AS #1
i = 0
DO UNTIL EOF(1)
    REDIM azim(i) AS STRING, Vgor(i) AS STRING, Vver(i) AS STRING, kurs(i) AS STRING, vis(i) AS STRING, daln(i) AS STRING
    INPUT #1, azim(i), Vgor(i), Vver(i), kurs(i), vis(i), daln(i)
    i = i + 1
LOOP
CLOSE
N = i - 1
 
'nazna4enie veli4ini massiva end
 
OPEN "data.txt" FOR INPUT AS #1
i = 0
DO UNTIL EOF(1)
    INPUT #1, azim(i), Vgor(i), Vver(i), kurs(i), vis(i), daln(i)
    PRINT azim(i), Vgor(i), Vver(i), kurs(i), vis(i), daln(i)
    i = i + 1
LOOP
CLOSE
 
'vvod dannyh v massiv end
 
PRINT "--------------", N
 
DIM x(N), y(N), z(N) AS LONG
DIM VgorX(N), VgorY(N) AS LONG
 
FOR i = 0 TO N
    x(i) = VAL(daln(i)) * SIN(VAL(azim(i)) * (3.1415926535# / 180))
    y(i) = VAL(daln(i)) * COS(VAL(azim(i)) * (3.1415926535# / 180)) * (-1)
    z(i) = VAL(vis(i))
    VgorX(i) = VAL(Vgor(i)) * SIN(VAL(kurs(i)) * (3.1415926535# / 180))
    VgorY(i) = VAL(Vgor(i)) * COS(VAL(kurs(i)) * (3.1415926535# / 180)) * (-1)
    PRINT VgorX(i), VgorY(i)
NEXT
 
DIM minDist AS DOUBLE
DIM RazDist(N, N) AS LONG
DIM first, second AS INTEGER
 
FOR i = 0 TO N
    FOR j = 0 TO N
        RazDist(i, j) = SQR((x(i) - x(j)) ^ 2 + (y(i) - y(j)) ^ 2)
    NEXT
NEXT
 
 
minDist = 9999999999#
FOR i = 0 TO N
    FOR j = 0 TO N
        IF RazDist(i, j) > 0 THEN
            IF minDist > RazDist(i, j) THEN
                minDist = RazDist(i, j)
                first = i
                second = j
            END IF
        END IF
    NEXT
NEXT
 
 
DIM shiftX, shiftY AS LONG
shiftX = 530
shiftY = 300
M = 300 / 200000
 
SCREEN 12
WIDTH 80, 60
WINDOW SCREEN (0, 0)-(799, 599)
 
PSET (0 + shiftX, 0 + shiftY), 15 ' centr koordinat gor oloskosti
CIRCLE (0 + shiftX, 0 + shiftY), 3, 15
LINE (270, 0)-(270, 800), 2
LINE (0 + shiftX, 0 + shiftY)-(shiftX, 50), 15
 
RAD = 10000
FOR i = 0 TO 14
    CIRCLE (0 + shiftX, 0 + shiftY), M * RAD, 8
    RAD = RAD + 10000
NEXT
 
ON ERROR GOTO Er
FOR i = 0 TO N
 
 
    PSET (x(i) * M + shiftX, y(i) * M + shiftY), 2
    CIRCLE (x(i) * M + shiftX, y(i) * M + shiftY), 5, 2
    LINE (x(i) * M + shiftX, y(i) * M + shiftY)-(x(i) * M + VgorX(i) + shiftX, y(i) * M + VgorY(i) + shiftY), 2
    Dr = -1
    LOCATE INT((y(i) * M + shiftY) / 10) - 1, INT((x(i) * M + shiftX) / 10) - 2
    IF Dr THEN PRINT i + 1
    CIRCLE (x(first) * M + shiftX, y(first) * M + shiftY), 7, 4
    CIRCLE (x(second) * M + shiftX, y(second) * M + shiftY), 7, 4
 
NEXT
ON ERROR GOTO 0
 
LOCATE 1, 1:
COLOR 2
FOR i = 0 TO N
    PRINT "Target:";
    PRINT i + 1
    PRINT "..........................."
    PRINT "BRG: ";
    PRINT azim(i);
    PRINT " SPD: ";
    PRINT Vgor(i)
    PRINT "CLB: ";
    PRINT Vver(i);
    PRINT " CRS: ";
    PRINT kurs(i)
    PRINT "HEIGHT: ";
    PRINT vis(i);
    PRINT " DST: ";
    PRINT daln(i)
    PRINT "---------------------------"
NEXT
 
COLOR 4
PRINT "Closest Targets"
PRINT "Target: ";
PRINT first + 1
PRINT "Target: ";
PRINT second + 1
 
 
 
SLEEP
END
 
Er: Dr = 0: RESUME NEXT

(VAL(kurs(j)) > 90 AND VAL(kurs(j)) <180)

(VAL(kurs(j)) < 90 OR VAL(kurs(j)) >180)

t = 6
maxdaln = 150001
FOR i = 0 TO t
    s(i) = 0
    FOR j = 0 TO N
    tx(i, j) = x(j) - (VgorX(j) * i)
    ty(i, j) = y(j) - (VgorY(j) * i)
    razdist(i, j) = SQR((x(j) - tx(i, j)) ^ 2 + (y(j) - ty(i, j)) ^ 2)
    IF (x(j) > 0 AND y(j) > 0) AND (VAL(kurs(j)) < 90 AND VAL(kurs(j)) > 0) THEN
    tdaln(j) = VAL(daln(j)) + razdist(i, j)
    ELSEIF (x(j) > 0 AND y(j) > 0 AND VAL(kurs(j)) > 90 OR VAL(kurs(j)) < 0) THEN
    tdaln(j) = VAL(daln(j)) - razdist(i, j)
    IF maxdaln > tdaln(j) THEN
    s(i) = s(i) + 1
    END IF
    END IF
    IF (x(j) > 0 AND y(j) < 0 AND VAL(kurs(j)) = 180) THEN
    tdaln(j) = VAL(daln(j)) + razdist(i, j)
    ELSEIF (x(j) > 0 AND y(j) < 0 AND (VAL(kurs(j)) > 180 OR 90 > VAL(kurs(j)))) THEN
    tdaln(j) = VAL(daln(j)) - razdist(i, j)
    IF maxdaln > tdaln(j) THEN
    s(i) = s(i) + 1
    END IF
    END IF
    IF x(j) < 0 AND y(j) < 0 AND 180 < VAL(kurs(j)) < 270 THEN
    tdaln(j) = VAL(daln(j)) + razdist(i, j)
    ELSEIF x(j) < 0 AND y(j) < 0 AND 180 > VAL(kurs(j)) > 270 THEN
    tdaln(j) = VAL(daln(j)) - razdist(i, j)
    IF maxdaln > tdaln(j) THEN
    s(i) = s(i) + 1
    END IF
    END IF
    IF x(j) < 0 AND y(j) > 0 AND 270 < VAL(kurs(j)) < 360 THEN
    tdaln(j) = VAL(daln(j)) + razdist(i, j)
    ELSEIF x(j) < 0 AND y(j) > 0 AND 270 > VAL(kurs(j)) > 360 THEN
    tdaln(j) = VAL(daln(j)) - razdist(i, j)
    IF maxdaln > tdaln(j) THEN
    s(i) = s(i) + 1
    END IF
    END IF
   
    IF s(i) > first THEN
    first = s(i)
    END IF
    NEXT
NEXT