这是我的
[x,y]=meshgrid(linspace(-2,1,30),linspace(-7,1,30));
z=y.^3/9+3*x.^2.*y+9.*x.^2+y.^2+x.*y+9;
surf(x,y,z)
% 查看点的范围
contour(x,y,z)
% 极小值
[x,f,h]= fminsearch(@(x)(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9),[0 0])
% 极大值
[x,f,h]=fminsearch(@(x)-(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9),[0 -5])