clear;

% set function handle
f_handle = @rosenbrock2d;

% Set up tolerance
atol = 1e-11;

% set starting point
x0 = [-1;0];

% call the quasi_newton_BFGS function
x = quasi_newton_BFGS(f_handle, x0, atol);

% print result of BFGS method
fprintf('Lösung: x=(%.6f,%.6f)\n',x(1),x(2))


