def nutation_matrix(epoch):
 (epsa, dpsi, deps) = nutation_components2000B(epoch.jd)
 return matrix_product(rotation_matrix((- (epsa + deps)), u'x', False), rotation_matrix((- dpsi), u'z', False), rotation_matrix(epsa, u'x', False))

def nutation_matrix ( epoch ) : ( epsa , dpsi , deps ) = nutation_components2000B ( epoch . jd ) return matrix_product ( rotation_matrix ( ( - ( epsa + deps ) ) , u'x' , False ) , rotation_matrix ( ( - dpsi ) , u'z' , False ) , rotation_matrix ( epsa , u'x' , False ) )
nutation matrix generated from nutation components .