exo2h.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 2

question h)

> f:=(x,y)->piecewise(x>0,exp(y*ln(x)),NULL);

f := proc (x, y) options operator, arrow; piecewise(0 < x, exp(y*ln(x)), NULL) end proc

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

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

> Limit('f'(x,exp(-1/x)), x=0)=limit(f(x,-1/ln(x)),x=0);

Limit(f(x, exp(-1/x)), x = 0) = undefined

Les limites suivants les courbes y=0 et y=-1/ln(x)  différent, il n'y a donc pas de limite globale

Profil suivant y=0

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

[Plot]

> coupe_plan(f,[t,0],[0,1]); courbe t->f(t,0)

[Plot]

Profil suivant y=x^2-x

> coupe(f,[0,1],[-1,2],[t,-1/ln(t)],[0,0.5],-(y+1/ln(x)),-50,100);

[Plot]

> coupe_plan(f,[t,-1/ln(t)],[0,1]); courbe t->f(t,-1/ln(t))

[Plot]

Représentation

> coupe(f,[0,1],[-1,2],[t,-1/ln(t)],[0,0.5],1,-10,50);

[Plot]

> profilx:=theta->t;profily:=theta->(1-theta)*0+theta*(-1/ln(t));t='t':

profilx := proc (theta) options operator, arrow; t end proc

profily := proc (theta) options operator, arrow; (1-theta)*0-theta/ln(t) end proc

Voici l'évolution de la route tracée sur la surface lorsque le trcé plan passe de la droite y=x à la courbe y=x^3

> profil_courbe(f,profilx,profily,[0,0.5],[0,1],-5);

[Plot]

> F:=proc(theta)  coupe_plan(f,[profilx(theta),profily(theta)],[0,0.5]); end:

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

[Plot]

>