def vecteur(xA, yA, xB, yB): return (xB - xA, yB - yA) def milieu(xA, yA, xB, yB): return ((xA + xB) / 2, (yA + yB) / 2) def recherche(): (xAB, yAB) = vecteur(5, 10, 9, -2) (xK, yK) = milieu(5, 10, 9, -2) (xAC, yAC) = vecteur(5, 10, 17, 4) (xJ, yJ) = milieu(5, 10, 17, 4) (xBC, yBC) = vecteur(9, -2, 17, 4) (xI, yI) = milieu(9, -2, 17, 4) for xD in range(5, 18): for yD in range(-2, 10): (xID, yID) = vecteur(xI, yI, xD, yD) (xJD, yJD) = vecteur(xJ, yJ, xD, yD) (xKD, yKD) = vecteur(xK, yK, xD, yD) prod1 = xKD * xAB + yKD * yAB prod2 = xJD * xAC + yJD * yAC prod3 = xID * xBC + yID * yBC if prod1 == 0 and prod2 == 0 and prod3 == 0: return (xD, yD)