mimmax.wxmx - Calcula e classifica os pontos críticos de f
(%i1) | kill(all); |
(%i1) | f:x^3+3·x·y^2−15·x−12·y; /· introduza aqui a função ·/ |
(%i4) | fx:diff(f,x)$ fy:diff(f,y)$ display(fx,fy)$ |
(%i5) | sol:algsys([fx=0, fy=0], [x,y]); |
(%i6) | n:length(sol); |
(%i7) |
for i:1 thru n do( P[i]:subst(sol[i],[x,y]), display(P[i]) ); |
(%i9) | define(fxx(x,y),diff(fx,x)); define(D(x,y),diff(fx,y)^2−diff(fx,x)·diff(fy,y)); |
(%i11) |
for i:1 thru n do( D[i]:D(P[i][1],P[i][2]),fxx[i]:fxx(P[i][1],P[i][2]), disp("−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−"), if D[i]>0 then display([P[i],D[i]]) else display([P[i],D[i],fxx[i]]), if D[i]=0 then disp("indeterminado") else if D[i]>0 then print(P[i],"é ponto de sela") else if fxx[i]>0 then print(P[i], "é ponto de mínimo local") else print(P[i],"é ponto de máximo local") )$ disp("−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−")$ |