from cmath import * def Droites_orthogonales(a,b,c): z=(c-a)/(b-a) if phase(z)==pi/2 or phase(z)==-pi/2: return True else: return False