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