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

%options=optimset(options,'LineSearch','cubicpoly');
% 采用本行參數(shù)設(shè)置時，表示計算方法為三次樣條內(nèi)插搜索法
% 否則為混合多項式內(nèi)插搜索法

%以下部分在思路一的各算法程序中相同，后面程序中不再重復(fù)
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);