The method of integration of dynamical systems admitting non-linear superpositions is applied to four-dimensional non-linear dynamical systems. All four-dimensional dynamical systems admitting non-linear superpositions with four-dimensional Vessiot-Guldberg-Lie algebras are classified into 160 standard forms. The integration method is described and illustrated.