th=sym('th'); % or syms th phi
phi= sym('phi');
R=rotx(phi)*rotz(th)

% if th = 0 deg

th = 0;
R1 =rotx(phi)*rotz(th)
