contour actif avec méthode de ballon

Bonjour,
je travaille sur le contour actif et je n'arrive pas à programmer cette fonction sous matlab .son principe est le suivant:
La fonction normale_aff renvoie à partir de trois points, la normale au deuxième point. Cette fonction est notamment utile pour calculer la force de ballon, qui est dirigée suivant la normale au contour. La normale au deuxième point peut avoir deux direction opposées, suivant que l'on veut une normale orientée vers l'intérieur ou l'extérieur. Nous avons donc décidé donner un sens de parcours du contour et l'ordre dans lequel on donne les trois points détermine ainsi la normale extérieure au contour.
Merci de m'aider