  We propose a Constraint Molecular Dynamics model for Fermionic system. In
this approach the equations of motion of wave packets for the nuclear many-body
problem are solved by imposing that the one-body occupation probability
$\bar{f}(r,p,t)$ can assume only values less or equal to 1. This condition
reflects the Fermionic nature of the studied systems and it is implemented with
a fast algorithm which allows also the study of the heaviest colliding system.
The parameters of the model have been chosen to reproduce the average binding
energy and radii of nuclei in the mass region $A=30\sim 208$. Some comparison
to data is given.
