# Ejercicio 11.13

kill(all)$

## Apartado a)

f1(x,y):= x^2+x*y+y^2-6*x+3;

define ( Dxf1(x,y) , diff(f1(x,y),x) );
define ( Dyf1(x,y) , diff(f1(x,y),y) );

Puntos estacionarios

Sist01: [Dxf1(x,y)=0 , Dyf1(x,y)=0]$
solve(Sist01);

Ahora estudiaremos de que tipo de punto se trata

'Hf1(x,y)=Hf1xy:hessian( f1(x,y) , [x,y]);

'Hf1(P)=Hf1P:subst([x=4,y=-2],Hf1xy);

detHf1P=determinant(Hf1P);

Como que el determinante es positivo y el primer menor es igual a 2 y, por lo tanto, también positivo, la forma cuadrática asociada a HfP es definida positiva, es decir,
f(x,y)>f(P) para todo(x,y).
Por lo tanto, el extremo es un mínimo local.

wxplot3d(f1(x,y), [x,3,5], [y,-3,-1], [legend, "f1(x,y)"] ,
 [gnuplot_preamble, "set pm3d at s; unset surf; unset colorbox"])$

## Apartado b)

f2(x,y):= 4*x+2*y-x^2+x*y-y^2;

define ( Dxf2(x,y) , diff(f2(x,y),x) );
define ( Dyf2(x,y) , diff(f2(x,y),y) );

Puntos estacionarios

Sist01: [Dxf2(x,y)=0 , Dyf2(x,y)=0]$
solve(Sist01);

Ahora estudiaremos de que tipo de punto se trata

'Hf2(x,y)=Hf2xy:hessian(f2(x,y) , [x,y]);

'Hf2(P)=Hf2P:subst([x=10/3,y=8/3],Hf2xy);

detHf2P=determinant(Hf2P);

Como que el determinante es positivo y el primer menor es igual a -2 y, por lo tanto, negativo, la forma cuadrática asociada a HfP es definida negativa, es decir,
g(x,y)<g(P) para todo(x,y).
Por lo tanto, el extremo es un máximo local.

wxplot3d(f2(x,y), [x, 0,7], [y,-5,5], [legend, "f2(x,y)"] ,
 [gnuplot_preamble, "set pm3d at s; unset surf; unset colorbox"])$