Variational method based on implicit integration scheme to simulate dynamics of quantum many-body systems. Developed in the context of TP-IV project in the Computational Quantum Science Lab @ EPFL led by prof. Giuseppe Carleo.