def triangle_normal(a, b, c):
normal = np.cross(a - c, b - c)<EOL>norm = np.linalg.norm(normal)<EOL>return normal/norm<EOL>
Return a vector orthogonal to the given triangle

       Arguments:
         a, b, c  --  three 3D numpy vectors