docs: def nutation(date, eop_correction=True, terms=106):  # pragma: no cover
      """Nutation as a rotation matrix
      """
      epsilon_bar, delta_psi, delta_eps = np.deg2rad(_nutation(date, eop_correction, terms))
      epsilon = epsilon_bar + delta_eps
  
      return rot1(-epsilon_bar) @ rot3(delta_psi) @ rot1(epsilon)