import math


def point_az(x1, x2, y1, y2):
    """
    Determines azimuth between two points (in cartesian coordinates).

    Parameters
    ----------
    x1 : float
        point 1 x
    x2 : float
        point 2 x
    y1 : float
        point 1 y
    y2 : float
        point 2 y

    Returns
    -------
    float
        Azimuth between points (in degrees)

    """

    azimuth1 = math.degrees(math.atan2((x2 - x1), (y2 - y1)))
    azimuth = (azimuth1 + 360) % 360
    return azimuth