exo5c.mw

> restart;with(plots):with(plottools):with(linalg);

Warning, the name changecoords has been redefined

Warning, the assigned name arrow now has a global binding

Warning, the protected names norm and trace have been redefined and unprotected

[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...

>

Exercice 5

question c)

> f:=(x,y)->piecewise(y=0,1,tanh(x^2/y^2));

f := proc (x, y) options operator, arrow; piecewise(y = 0, 1, tanh(x^2/y^2)) end proc

>

Représentation

> plot3d(f(x,y),x=-2..2,y=-2..2,axes=normal);

[Plot]

>

Voici le profil de la section obtenu par une droite passant par l'origine et tournant au cours du temps

> profil_x:=theta->cos(theta)*t;profil_y:=theta->sin(theta)*t;t:='t':

profil_x := proc (theta) options operator, arrow; cos(theta)*t end proc

profil_y := proc (theta) options operator, arrow; sin(theta)*t end proc

> sect:=theta->cos(theta)*x+y*sin(theta);

sect := proc (theta) options operator, arrow; cos(theta)*x+y*sin(theta) end proc

> profil_courbe(f,profil_x,profil_y,[-5,5],[-Pi,Pi],0);

[Plot]

>

Le même profil vu sur un plan fixe (courbe du paramétrage)

> F:=proc(theta)  plot(f(profil_x(theta),profil_y(theta)),t=-5..5 ); end:

> animate(F,[theta],theta=-Pi..Pi, axes=normal,frames=41);

[Plot]

Voici les courpes obtenues par les sections correspondantes à ces droites tournantes

> profil_section(f,[-5,5],[-5,5],sect,[-Pi,Pi],-10);

[Plot]

>

>

En particulier les mimites diffèrent suivant y=x et y=0

profil selon y=0

> coupe(f,[-5,5],[-5,5],[t,0],[-5,5],y,-1,100);

[Plot]

> Limit('f'(x,0),x=0)=limit('f'(x,0),x=0);

Limit(f(x, 0), x = 0) = 1

profil selon y=x

> coupe(f,[-5,5],[-5,5],[t,t],[-5,5],y-x,-1,100);

[Plot]

> Limit('f'(x,x),x=0)=limit('f'(x,x),x=0);

Limit(f(x, x), x = 0) = tanh(1)

>