# Path: skyfield/constants.py
# ASEC2RAD = 4.848136811095359935899141e-6
# 
# T0 = 2451545.0


from numpy import array, cos, sin
from .constants import ASEC2RAD, T0


def compute_precession(jd_tdb):
    """Return the rotation matrices for precessing to an array of epochs.

    `jd_tdb` - array of TDB Julian dates

    The array returned has the shape `(3, 3, n)` where `n` is the number
    of dates that have been provided as input.

    """
    eps0 = 84381.406

    # 't' is time in TDB centuries.


# Path: skyfield/constants.py
# ASEC2RAD = 4.848136811095359935899141e-6
# 
# T0 = 2451545.0


# Path: skyfield/precessionlib.py
from numpy import array, cos, sin
from .constants import ASEC2RAD, T0



def compute_precession(jd_tdb):
    """Return the rotation matrices for precessing to an array of epochs.

    `jd_tdb` - array of TDB Julian dates

    The array returned has the shape `(3, 3, n)` where `n` is the number
    of dates that have been provided as input.

    """
    eps0 = 84381.406

    # 't' is time in TDB centuries.


    t = (jd_tdb - T0) / 36525.0