Revise paste sanislav

syms a b c alfaA alfaB alfaC betaA betaB betaC gamaA gamaB gamaC
syms xA yA zA xB yB zB xC yC zC
ax=a*cos(alfaA);
ay=a*cos(betaA);
az=a*cos(gamaA);
bx=b*cos(alfaB);
by=b*cos(betaB);
bz=b*cos(gamaB);
cx=c*cos(alfaC);
cy=c*cos(betaC);
cz=c*cos(gamaC);
Rx = ax + bx + cx;
Ry = ay + by + cy;
Rz = az + bz + cz;
R = [Rx Ry Rz];
Mx = yA*az-zA*ay+yB*bz-zB*by+yC*cz-zC*cy;
My = zA*ax-xA*az+zB*bx-xB*bz+zC*cx-xC*cz;
Mz = xA*ay-yA*ax+xB*by-yB*bx+xC*cy-yC*cx;
Mo = [Mx My Mz];
disp(R)
disp(Mo)
an=input('a= ');
bn=input('b= ');
cn=input('c= ');
xAn=input('xA= ');
yAn=input('yA= ');
zAn=input('zA= ');
xBn=input('xB= ');
yBn=input('yB= ');
zBn=input('zB= ');
xCn=input('xC= ');
yCn=input('yC= ');
zCn=input('zC= ');
alfaAn=input('alfaA= ');
betaAn=input('betaA= ');
gamaAn=input('gamaA= ');
alfaBn=input('alfaB= ');
betaBn=input('betaB= ');
gamaBn=input('gamaB= ');
alfaCn=input('alfaC= ');
betaCn=input('betaC= ');
gamaCn=input('gamaC= ');
alfaAn=deg2rad(alfaAn);
betaAn=deg2rad(betaAn);
gamaAn=deg2rad(gamaAn);
alfaBn=deg2rad(alfaBn);
betaBn=deg2rad(betaBn);
gamaBn=deg2rad(gamaBn);
alfaCn=deg2rad(alfaCn);
betaCn=deg2rad(betaCn);
gamaCn=deg2rad(gamaCn);
R_n = eval(subs(R, [a b c alfaA alfaB alfaC betaA betaB betaC gamaA gamaB gamaC], [an bn cn alfaAn alfaBn alfaCn betaAn betaBn betaCn gamaAn gamaBn gamaCn]))
Mo_n = eval(subs(Mo, [a b c alfaA alfaB alfaC betaA betaB betaC gamaA gamaB gamaC xA xB xC yA yB yC zA zB zC],[an bn cn alfaAn alfaBn alfaCn betaAn betaBn betaCn gamaAn gamaBn gamaCn xAn xBn xCn yAn yBn yCn zAn zBn zCn]))