x=0:0.1:1;
y=x;
[X,Y]=meshgrid(x,y);
Z=X+Y*(1-0.5*X)*(1+X).^2-0.5*X.^2-0.5*Y.^2;
>> mesh(X,Y,Z)
看你采用什么优化算法了
