Most of the methods used to derive deterministic equations governing t
he evolution of linear waves in random media are based on the mean fie
ld approach. For a given linear system of equations with random coeffi
cients this approach results in an approximate deterministic equation
(or system) for the mean (averaged over the set of all realizations) f
ield. The coefficients of such equations are associated with the stati
stical moments of the random coefficients of the initial system. Somet
imes similar approaches are also applied to nonlinear problems. Howeve
r, recently, in a number of examples, it was shown that the mean field
approach in nonlinear problems may give the wrong results. The error
is related to the infinite growth of the root-mean-square phase fluctu
ations due to fluctuations of the wave velocity in an inhomogeneous me
dium. To overcome this effect the idea of eliminating the unbounded gr
owth of phase fluctuations by using an appropriate coordinate transfor
mation was proposed resulting in a method called the ''mean waveform m
ethod''. In the present paper we extend the idea of the mean waveform
method in order to develop a systematic approach which enables the con
struction of an approximate deterministic evolution equation for a giv
en quasi-hyperbolic and quasi-linear system of equations with weak non
linearity and stationary random coefficients. The applicability of the
proposed approach is demonstrated by means of an example.