clc, clear all, close all
I = [-10 4 -10 4];
f = @(x, y) 2*log(x.^2 + y.^2 - 1) + x + 2*y;
fsurf(f, I)
H = [1/2, -1; -1 -1]
eig(H)