syms fi ;
syms theta;
syms psi;
rM=[cos(fi)*cos(theta)*cos(psi)-sin(fi)*sin(psi) -cos(fi)*cos(theta)*sin(psi)-sin(fi)*cos(psi) cos(fi)*sin(theta) ;
sin(fi)*cos(theta)*cos(psi)+cos(fi)*sin(psi) -sin(fi)*cos(theta)*sin(psi)+cos(fi)*cos(psi) sin(fi)*sin(theta);
-sin(theta)*cos(psi) sin(theta)*cos(psi) cos(theta)];
transition=[x y z]';
zeros=sym(zeros(1,3));
one=sym(ones(1));
hM=[rM transition ;
zeros one];
hM=double(hM);
Best Answer