Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement