Ekstremum lokalne w wykresie funkcji

syms x y %deklarujemy zmienne
f=x^3+y^3-9*x*y;
fx=diff(f,x)
fy=diff(f,y)%obliczamy pochodne cząstkowe
[X,Y]=solve(fx,fy)
fxx=diff(fx,x)%tworzymy hesjan i funkcje g
fyx=diff(fx,y)
fxy=diff(fy,x)
fyy=diff(fy,y)
H=[fxx fyx;fyx fyy]
g=det(H)
subs(g,{x,y},{X(1),Y(1)})
subs(g,{x,y},{X(2),Y(2)})
subs(fxx,{x,y},{X(2),Y(2)})
subs(f,{x,y},{X(2),Y(2)})