%test for steepest descent
y=3;
fun = inline('2*x(1)^2 + 136*x(2)^2 - 10*x(1)*x(2) ');
%fun = inline('exp(x(1)) * (4*x(1)^2 + 2*x(2)^2 + 4*x(1)*x(2) + 2*x(2) + 1)');
%fun = inline('(x(1)) * (5*x(1)^2 + 2*x(2)^2 + 4*x(1)*x(2) + 2*x(2) + 1)');
OPTIONS=[];
OPTIONS=optimset('HessUpdate','steepdesc');
%x0 = [-1 1];
%options = optimset('LargeScale','off');
[x, fval, exitflag, output] = fminunc(fun, x0, options);

