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)
def precesion ( date ) : zeta , theta , z = np . deg2rad ( _precesion ( date ) ) return rot3 ( zeta ) @ rot2 ( - theta ) @ rot3 ( z )
Precession as a rotation matrix