function H = df2(x)
% gibt Hessematrix der Himmelblaufunktion aus
%
% Eingabe:
%   x...2-dimensionaler Vektor
%
% Ausgabe:
%   H...2x2 Matrix
%
% Die Funktion besitzt Nullstellen bei
% (3, 2)
% (-2.8051, 3.1313)
% (-3.7793, -3.2832)
% (3.5844, -1.8481)
% (-0.2708,-0.9230)
%
% Markus Grasmair,
% Wien, Mai 2011

X = x(1);
Y = x(2);

H = eye(2);
H(1,1) = 12*X.^2 + 4*Y-42;
H(1,2) = 4*X+4*Y;
H(2,1) = H(1,2);
H(2,2) = 12*Y.^2 + 4*X - 26;
