returns rotation array along z-axis y-axis -z-axis
theta (float) – rotation angle
phi (float) – rotation angle
res – rotation array
(3,3)-np.ndarray