A method of integration of non-stationary dynamical systems admitting nonlinear superpositions is presented. The method does not require knowledge of symmetries of the differential equations under consideration. The integration procedure is based on classification of Vessiot-Guldberg-Lie algebras associated with nonlinear superpositions. It is shown that the systems associated with one-and two-dimensional Lie algebras can be integrated by quadrature upon introducing Lie's canonical variables. It is not necessary to know symmetries of a system in question in this approach. Two-dimensional non-stationary dynamical systems with three-dimensional Vessiot-Guldberg-Lie algebras are classified into thirteen standard forms. Ten of them are integrable by quadrature. The remaining three standard forms lead to the Riccati equations. Integration of perturbed dynamical systems possessing approximate nonlinear superposition is discussed.