%start_octave_ferdavel 
%
%  nelder-mead á bls, 64
%
function y = f(x);
y= -100*(x(2)-x(1)^2)^2-(1-x(1))^2;
endfunction;
%
% Til að fá lausn eins og í bók þurfti ég að gera
%
x=[1;2];
n=length(x);
%
% og búa svo til simplexið
%
S = [x x(:,ones(n,1))+diag(0.05*x+0.00025*(abs(x)<eps))];

x=neldmead('f',x,S);
