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