We propose a constructive method to prove the integrability of a given physical Hamiltonian in one dimension, which amounts to looking for appropriate polynomial R-matrices, solutions of GYBE, whose first coefficient in the power expansion with respect to the spectral parameter is the Hamiltonian itself. The method is applied to the extended Hubbard Hamiltonian, in particular to the cases in which it exhibits so(4) or gl(2,1) symmetries. We show that in the latter case the R-matrices are at most polynomial of second degree, whose coefficients are nothing but the Hamiltonian, the identity and the permutation operator. In this way, all known completely integrable cases are recovered. Also, the method allows to recognize that the possible integrability of the most general gl(2,1) invariant Hamiltonian depends on the existence of a non-additive R-matrix.