Implementation of the dual only potential reduction and primal dual method for the determination of the 2DM for the 1D Hubbard model. The Hessian for the solution of the Newton equations is explicitly constructed