NAME: lineintersection
PURPOSE: Calcule les coordonnees de l''intersection de 2 droites ou
d''une serie de 2 droites
CATEGORY:petit truc qui peut etre utile (sans boucles, ca va de soit!)
CALLING SEQUENCE: point=lineintersection(abc1, abc2)
INPUTS: abc est un tableau de dimensions 3,
nbre_de_couple_de_droites, dont chaque ligne contient les 3
parametres a, b, c de l''equation de droite du type ax+by+c=0
KEYWORD PARAMETERS:
/FLOAT: pour retourner l''output sous forme de tableau de
reel plutot que de vecteur decomplexes (par defaut)
OUTPUTS:2 possibilites:
1) par defaut: c''est une vecteur de complexe dont chaque
element est les coordonnees du point d''intersection d''un
couple de droites.
2) si /FLOAT est active, c'est un tableau de reels de
dimensiones 2, nbre_de_couple_de_droites dont chaque
ligne est les coordonnees du point d''intersection d''un
couple de droites
COMMON BLOCKS:
SIDE EFFECTS:si les deux droites dont paralleles, on retourne les
coordonnes (!values.f_nan,!values.f_nan)
RESTRICTIONS:attention a la precision de la machine qui fait que
les coordonnees calcules de verifient peut-etre pas exactement les
equations du couple de droites.
EXAMPLE:
IDL> abc1=linearequation(complex(1,2),[3,4])
IDL> abc2=linearequation(complex(1,2),[8,15])
IDL> print, lineintersection(abc1, abc2)
( 1.00000, 2.00000)
IDL> print, lineintersection(abc1, abc2,/float)
1.00000 2.00000