  A non-Abelian theory of fermions interacting with gauge bosons, the
constrained system, is studied. The equations of motion for a singular system
are obtained as total differential equations in many variables. The
integrability conditions are investigated and the set of equations of motion is
integrable. The Senjanovic and the canonical methods are used to quantize the
system, and the integration is taken over the canonical phase space
coordinates.
