  In the analytic-Bargmann representation associated with the harmonic
oscillator and spin coherent states, the wavefunction as entire complex
functions can be factorized in terms of their zeros in a unique way. The
Schr\"odinger equation of motion for the wavefunction is turned to a system of
equations for its zeros. The motion of these zeros as a non-linear flow of
points is studied and interpreted for linear and non-linear bosonic and spin
Hamiltonians. Attention is given to the study of the zeros of the
Jaynes-Cummings model and to its finite analoque. Numerical solutions are
derived and discussed.
