fun='(x(1)^2-2*x(1))*exp(-x(1)^2-x(1)*x(2)-x(2)^2)';
[xopt,fopt,exitflag,output]=fminsearch(fun,[1,1])
f='(x^2-2*x)*exp(-x^2-x*y-y^2)';
ezmeshc(f);
xlabel('\bf x');ylabel('\bf y');zlabel('\bf f(x,y)');
function z=fun1(x)
z=(x(1)^2-2*x(1))*exp(-x(1)^2-x(2)^2-x(1)*x(2))
end;
[z0 fval exitflag]=fminsearch('fun1',[-10 10]); %搜索范围自己设