  The charge and spin dynamics of the two-dimensional Hubbard model in the
paramagnetic phase is first studied by means of the two-pole approximation
within the framework of the Composite Operator Method. The fully
self-consistent scheme requires: no decoupling, the fulfillment of both Pauli
principle and hydrodynamics constraints, the simultaneous solution of fermionic
and bosonic sectors and a very rich momentum dependence of the response
functions. The temperature and momentum dependencies, as well as the dependency
on the Coulomb repulsion strength and the filling, of the calculated charge and
spin susceptibilities and correlation functions are in very good agreement with
the numerical calculations present in the literature.
