  We develop a consistent perturbation theory in quantum fluctuations around
the classical evolution of a system of interacting bosons. The zero order
approximation gives the classical Gross-Pitaevskii equations. In the next order
we recover the truncated Wigner approximation, where the evolution is still
classical but the initial conditions are distributed according to the Wigner
transform of the initial density matrix. Further corrections can be
characterized as quantum scattering events, which appear in the form of a
nonlinear response of the observable to an infinitesimal displacement of the
field along its classical evolution. At the end of the paper we give a few
numerical examples to test the formalism.
