<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def rotoreflection(axis, angle, origin=(0, 0, 0)):
 """ Returns a roto-reflection symmetry operation Args: axis (3x1 array):
 Axis of rotation / mirror normal angle (float):
 Angle in degrees origin (3x1 array):
 Point left invariant by roto-reflection. Defaults to (0, 0, 0). Return: Roto-reflection operation """

        rot = SymmOp.from_origin_axis_angle(origin, axis, angle)
        refl = SymmOp.reflection(axis, origin)
        m = np.dot(rot.affine_matrix, refl.affine_matrix)
        return SymmOp(m)