1
function d=projaff(x,y,pt)
2
// projaff - projection d'un point sur une droite
6
// x: vecteur des abscisses des 2 points definissant la droite
7
// y: vecteur des ordonnes des 2 points definissant la droite
8
// pt: coordonnes du point a projeter [abscisses ordonnees]
10
//l'equation de la droite est :
11
//(y2-y1)*(x-x1)-(x2-x1)*(y-y1)=0
12
//celle de la normale passant par le point:
13
//(x2-x1)*(x-xp)+(y2-y1)*(y-yp)=0
19
dx=x(2)-x(1);dy=y(2)-y(1)
21
d=m\ [-x(2)*y(1)+x(1)*y(2);m(2,:)*[pt(1);pt(2)]]