<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def rotmat(m1, angle, iaxis):
 """ Rotmat applies a rotation of angle radians about axis iaxis to a matrix. This rotation is thought of as rotating the coordinate system. http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/rotmat_c.html :param m1: Matrix to be rotated. :type m1: 3x3-Element Array of floats :param angle: Angle of rotation (radians). :type angle: float :param iaxis: Axis of rotation X=1, Y=2, Z=3. :type iaxis: int :return: Resulting rotated matrix. :rtype: 3x3-Element Array of floats """

    m1 = stypes.toDoubleMatrix(m1)
    angle = ctypes.c_double(angle)
    iaxis = ctypes.c_int(iaxis)
    mout = stypes.emptyDoubleMatrix()
    libspice.rotmat_c(m1, angle, iaxis, mout)
    return stypes.cMatrixToNumpy(mout)