Нахождение нормали к ребру прямоугольника

[dr.Chas]

Здравствуйте. Задача реализовать алгоритм Кируса-Бека. Это я сделал, он работает, но не могу найти нормаль к каждому ребру прямоугольника, не силён в геометрии.

на возникший вопрос о том как же я реализовал алгоритм и узнал что он работает? Сразу отвечаю, к готовому прямоугольнику взял вычисленную нормаль.

И так дан прямоугольник:


Из справочника:

Как рассчитывается нормаль? Для начала возьмем три точки (a,b и c). По этим точкам необходимо вычислить вектора p (из точки a в b) и q (из a в с):
p = b - a
q = c - a

Ну беру я три точки a,b,c
px:=bx-ax;
py:=by-by;
qx:=cx-ax;
qy:=cy-ay;

Получив два этих вектора вычисляем нормаль как перекрестное произведение векторов p и q с нормализацией результата. Это можно представить в псевдокоде как:
n = Cross(p, q);
Normalize(n);

Вычисляю перекрёстное произведение векторов:
I := qx*py-px*qy;
По-моему здесь я делаю что то не правильно? Ну а дальше нормализация.

Для нормализации результата необходимо разделить каждый параметр вектора на его длину, которая вычисляется как:
length = Sqrt(x2+y2+z2);

Я в принципе сам бы сидел и мучался, может быть что то и "намучил", но времени не так много осталось.

[dr.Chas]

Что-то мне ни кто не хочет помогать. Не ужели ни кто не знает?

Вот переписал с си, на делфи.

prmx_[1]:=0;  prmx_[2]:=0; prmx_[3]:=7; prmx_[4]:=7;
prmy_[1]:=0;  prmy_[2]:=5; prmy_[3]:=5; prmy_[4]:=0;
{а если в prmx и prmy зададим размеры реального прямоугольника начерченного на форме то не  работает}

//wnormx, wnormy  - массив с нормалями
i:=0;  
vox:=prmx_[4]; voy:=prmy_[4];   //последние координаты прямоугольника
j:=1;
   repeat
    inc(i);
    if i>3 then j:=0;
    vix:=prmx_[i]; viy:=prmy_[i];
    vnx:=viy-voy;  vny:=vox-vix;
    vox:=vix; voy:=viy;
    wnormx[i]:=vnx; wnormy[i]:=vny;
   until j=0;

Работает но с не большим но. Если подставлять координаты т.е. (1,2) (декартовой системы координат). То всё работает. А если подставляю значения из делфи, то не работает. Почему?