We show how the anticommutation relations for Fermi operators can be implemented with computer algebra using SymbolicC++. We describe applications to the Hubbard model. An important identity for Fermi operators is proved. Then, we test for higher order constants of motion for the Hubbard model. Finally, the matrix representation for the four point Hubbard model is calculated.