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