import math


def normal(x1: float, y1: float, x2: float, y2: float) -> tuple[float, float]:
    """
    Compute the normal vector given two points
    """
    phi = math.atan2(y2 - y1, x2 - x1) + math.pi/2
    return (math.cos(phi), math.sin(phi))