| > | 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
![]()
![]()
![]()
| > |
| > |
Exercice 7
| > | f:=(x,y)->piecewise(x=0 and y=0,0,(x*y)/(x^2+y^2)); |
| > | plot3d(f(x,y),x=-0.5..0.5,y=-5..5,axes=normal); |
![[Plot]](images/exo7_6.gif)
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': |
| > | sect:=theta->cos(theta)*x+y*sin(theta); |
| > | profil_courbe(f,profil_x,profil_y,[-5,5],[-Pi,Pi],0); |
![[Plot]](images/exo7_10.gif)
| > |
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]](images/exo7_11.gif)
Voici les courpes obtenues par les sections correspondantes à ces droites tournantes
| > | profil_section(f,[-0.5,0.5],[-5,5],sect,[-Pi,Pi],-10); |
![[Plot]](images/exo7_12.gif)
Ainsi f n'est pas continue en (0,0) comme nous le montre le profil se lon y=x qui ne passe pas par la valeur f(0,0)=0
| > | coupe(f,[-0.5,0.5],[-5,5],[t,t],[-0.7,0.7],x-y,-1,100); |
![[Plot]](images/exo7_13.gif)
Cependant par calcul des taux d'acroissements de t->(f(0+th)-f(0))/t
| > |
| > | 'f'(t*h1,t*h2)=simplify(f(t*h1,t*h2)); |
| > | Limit('f'(t*h1,t*h2)/t,t=0)=limit(f(t*h1,t*h2)/t,t=0); |
La dérivée dans la direction (h1,h2) n'existe que si h_1=0 ou h2=0, i.e. dans la direction i ou j
| > |