глюкануло
матрица Гессе
(8/x1^3 0 1/x1^2)
( 0 8/x2^3 1/x2^2)
( 1/x1^2 1/x2^2 0 )
получил как
( d2f/dx1^2 d2f/dx1dx2 d2f/dx1dk )
( d2f/dx2dx1 d2f/dx2^2 d2f/dx2dk )
( d2f/dkdx1 d2f/dkdx2 d2f/dk^2 )
а потом просто подставил k= -4