This algorithm calculates the mutually conjugate directions of search
with respect to the Hessian matrix of
directly from the function evaluation and the gradient evaluation, but
without the direct evaluation of the Hessian of the function .
The quantity
is added to the gradient at each iteration, and when is a quadratic form
(positive definite), this results in a set of mutually conjugate vectors.