I derive the equation of motion in molecular dynamics for doing full l
attice QCD simulations with clover quarks. The even-odd preconditionin
g technique, expected to significantly reduce the computational effort
, is further developed for the simulations.