syms a b c real
syms xA yA zA xB yB zB xC yC zC real
syms alfaA betaA gamaA alfaB betaB gamaB alfaC betaC gamaC real
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);
av = [ax ay az].';
bv = [bx by bz].';
cv = [cx cy cz].';
OA = [xA yA zA].';
OC = [xC yC zC].';
E1 = cross(OA,bv) + cross(av,OC);
E2 = cross(av,bv)+cv;
display(E1)
display(E2)
an = input('Modul a: ');
bn = input('Modul b: ');
cn = input('Modul 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 = deg2rad(input('alfaA: '));
betaAn = deg2rad(input('betaA: '));
gamaAn = deg2rad(input('gamaA: '));
alfaBn = deg2rad(input('alfaB: '));
betaBn = deg2rad(input('betaB: '));
gamaBn = deg2rad(input('gamaB: '));
alfaCn = deg2rad(input('alfaC: '));
betaCn = deg2rad(input('betaC: '));
gamaCn = deg2rad(input('gamaC: '));
E1n = subs(E1, [a b c xA yA zA xB yB zB alfaA betaA gamaA alfaB betaB gamaB], [an bn cn xAn yAn zAn xBn yBn zBn alfaAn betaAn gamaAn alfaBn betaBn gamaBn]);
E1n = eval(E1n);
display(E1n)
E2n = subs(E2, [a b c alfaA betaA gamaA alfaB betaB gamaB alfaC betaC gamaC], [an bn cn alfaAn betaAn gamaAn alfaBn betaBn gamaBn alfaCn betaCn gamaCn]);
E2n = eval(E2n);
display(E2n)