code: def eul2m(angle3, angle2, angle1, axis3, axis2, axis1):
      
      angle3 = ctypes.c_double(angle3)
      angle2 = ctypes.c_double(angle2)
      angle1 = ctypes.c_double(angle1)
      axis3 = ctypes.c_int(axis3)
      axis2 = ctypes.c_int(axis2)
      axis1 = ctypes.c_int(axis1)
      r = stypes.emptyDoubleMatrix()
      libspice.eul2m_c(angle3, angle2, angle1, axis3, axis2, axis1, r)
      return stypes.cMatrixToNumpy(r)
docs: Construct a rotation matrix from a set of Euler angles.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/eul2m_c.html
  
      :param angle3: Rotation angle about third rotation axis (radians).
      :type angle3: float
      :param angle2: Rotation angle about second rotation axis (radians).
      :type angle2: float
      :param angle1: Rotation angle about first rotation axis (radians).
      :type angle1: float
      :param axis3: Axis number of third rotation axis.
      :type axis3: int
      :param axis2: Axis number of second rotation axis.
      :type axis2: int
      :param axis1: Axis number of first rotation axis.]
      :type axis1: int
      :return: Product of the 3 rotations.
      :rtype: 3x3-Element Array of floats