%
%
%
function [y,dy]=f(x);
y=-100*(x(2)-x(1)^2)^2-(1-x(1))^2;
dy=[2*(1-x(1))+400*(x(2)-x(1)^2)*x(1);-200*(x(2)-x(1)^2)];
endfunction;

 x=[1;2]
%
%  Hér þurfti ég að bæta við
%
k = size(x,1);
A = zeros(0,0);
%
%  og gera svo
% 
x=qnewton('f',x,A)

% ég fékk rétt svar með því að 
%

optset('qnewton','searchmeth',2);

%
%  Einnig má prófa
%
optset('qnewton','maxit',5000);
optset('qnewton',searchmeth',3);


