def Rx_matrix(theta):
return np.array([<EOL>[<NUM_LIT:1>, <NUM_LIT:0>, <NUM_LIT:0>],<EOL>[<NUM_LIT:0>, np.cos(theta), -np.sin(theta)],<EOL>[<NUM_LIT:0>, np.sin(theta), np.cos(theta)]<EOL>])<EOL>
Rotation matrix around the X axis