In this section we will apply the approach of [19,20] to derive the evolution equation for the multi-mode PDF via introducing a generating functional, performing a weak-nonlinearity expansion and statistical averaging aided by a new graphical technique. We are going to demonstrate the phase independence property. This will also prepare us to answer the question of the next section: to what extend the amplitudes are going to remain statistically independent?
