def precesion(date):  # pragma: no cover
    """Precession as a rotation matrix
    """
    zeta, theta, z = np.deg2rad(_precesion(date))
    return rot3(zeta) @ rot2(-theta) @ rot3(z)
Precession as a rotation matrix