def rotation_matrix_2D(th):
s = np.sin(th)<EOL>c = np.cos(th)<EOL>return np.array([[c, -s], [s, c]])<EOL>
rotation_matrix_2D(th) yields a 2D numpy rotation matrix th radians counter-clockwise about the
origin.