def sun_rot_elements_at_epoch(T, d):
    """Calculate rotational elements for Sun.

    Parameters
    ----------
    T: float
        Interval from the standard epoch, in Julian centuries i.e. 36525 days.
    d: float
        Interval in days from the standard epoch.

    Returns
    -------
    ra, dec, W: tuple (float)
        Right ascension and declination of north pole, and angle of the prime meridian.

    """
    ra = 286.13
    dec = 63.87
    W = 84.176 + 14.1844000 * d

    return ra, dec, W