options=optimset('LargeScale','off');
options=optimset(options,'MaxFunEvals', 300);

options=optimset(options,'HessUpdate','steepdesc');
%options=optimset(options,'LineSearch','cubicpoly');
% 采用本行參數設置時，表示計算方法為三次樣條內插搜索法
% 否則為混合多項式內插搜索法
% 以下代碼同上
x0=[-1.9,2]
GRAD=[100*(4*x(1)^3-4*x(1)*x(2))+2*x(1)-2; …
100*(2*x(2)-2*x(1)^2); banplot2(x)]';
f='100*(x(2)-x(1)^2)^2+(1-x(1))^2';
options = optimset(options,'gradobj','on');
disp('[x,fval,exitflag,output] = fminunc({f,GRAD},x,options);');
[x,fval,exitflag,output] = fminunc({f,GRAD},x0,options);