is desingularized when the Hessian becomes singular. The matrix A is called the adjoint matrix for H. The full Hessian matrix can be computationally expensive at each step of the positions x(t). However, it can be updated . A first numerical step from a stationary point goes in direction f . The following steps then secure that the gradient maintains this direction . The solution curves are called Newton trajectories (NT).