Критерий степени подобия треугольников

REM
REM  Даны два треугольника с углами 
REM       (A1, B1, C1) и (A2, B2, C2)
REM  причем A1 >= A2, B1 >= B2, C1>= C2
REM  D - Степень подобия этих треугольников 
REM        D = |A1 - A2| + |B1 - B2| + |C1 - C2|
REM        При D = 0 треугольники подобны!!
REM 
 
DECLARE FUNCTION p! (x1!, y1!, x2!, y2!)
DECLARE FUNCTION u! (a!, b!, c!)
 
CLS
SCREEN 11
RANDOMIZE TIMER
CONST c1 = .4636476
CONST c2 = 1.107149
CONST c3 = 1.570796
DIM x(1 TO 40)
DIM y(1 TO 40)
 
FOR i = 1 TO 40
   x(i) = 400 * RND + 220
   y(i) = 400 * RND
   PSET (x(i), y(i))
NEXT i
 
LINE (20, 40)-(20, 400)
LINE -(200, 400)
LINE -(20, 40)
 
x3 = 300
y3 = 400
D = 7
 
FOR i = 1 TO 39
FOR j = i + 1 TO 40
   a = p(x(i), y(i), x(j), y(j))
   b = p(x(j), y(j), x3, y3)
   c = p(x3, y3, x(i), y(i))
   ua = u(a, b, c)
   ub = u(b, c, a)
   uc = u(c, a, b)
 
   IF ua < ub THEN SWAP ua, ub
   IF ua < uc THEN SWAP ua, uc
   IF ub < uc THEN SWAP ub, uc
 
   dd = ABS(ua - c3) + ABS(ub - c2) + ABS(uc - c1)
 
   IF D > dd THEN
      D = dd
      x1 = x(i)
      y1 = y(i)
      x2 = x(j)
      y2 = y(j)
   END IF
NEXT j, i
 
LINE (x1, y1)-(x2, y2)
LINE -(x3, y3)
LINE -(x1, y1)
 
END
 
FUNCTION p (x1, y1, x2, y2)
   p = SQR((x1 - x2) ^ 2 + (y1 - y2) ^ 2)
END FUNCTION
 
FUNCTION u (a, b, c)
   z = (b ^ 2 + c ^ 2 - a ^ 2) / (2 * b * c)
   IF 1 - z ^ 2 <= 0 THEN
      u = 1.570796
   ELSE
      u = 1.570796 - ATN(z / SQR(1 - z ^ 2))
   END IF
END FUNCTION
 
REM Автор предлагает вам удивительную по красоте программу
REM Эта программа - Откровение богов, к коим автор причисляет и себя
REM Вы - Счастливейшие из живущих под Солнцем.
REM Вам выпал шанс видеть Великое Творение Разума!
REM