1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
| function [x] = normaletangent(u,v)
a=5;
thetarad=u*(pi/180);
phirad=v*(pi/180);
syms thetarad phirad
Mx=a*(2*cos (thetarad) + cos (2*thetarad));
My=a*cos (phirad)*(2*sin (thetarad)+ sin (2*thetarad));
Mz=a*sin (phirad)*(2*sin (thetarad)+ sin (2*thetarad));
Mx1=simplify(diff(Mx,thetarad));
My1=simplify(diff(My,thetarad));
Mz1=simplify(diff(Mz,thetarad));
Mx2=simplify(diff(Mx,phirad));
My2=simplify(diff(My,phirad));
Mz2=simplify(diff(Mz,phirad));
x=subs(Mx,u,v)+subs(Mx1,u,v)+subs(Mx2,u,v);
y=subs(My,u,v)+subs(My1,u,v)+subs(Mz2,u,v);
z=subs(My,u,v)+subs(My1,u,v)+subs(Mz2,u,v);
disp('Les équations paramétriques du plan tangent sont :');
disp('x =');disp(x);
disp('y =');disp(y);
disp('z =');disp(z);
end |
Partager