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

%options=optimset(options,'LineSearch','cubicpoly');
% 采用本行參數設置時，表示計算方法為三次樣條內插搜索法
% 否則為混合多項式內插搜索法
x0=[-1.9,2];
JAC='[-20*x(1), 10; -1, 0; banplot2(x)]';
f='[10*(x(2)-x(1)^2),(1-x(1))]';
options = optimset (options,'Jacobian','on');
disp('[x,resnorm,residual,exitflag,output]=lsqnonlin({f,JAC},x0,[],[], options);');
[x,resnorm,residual,exitflag,output]=lsqnonlin({f,JAC},x0,[],[], options);
fval = resnorm;