This numerical toolbox can describe any multimode gaussian state, find its time evolution following open/closed dynamics arising from a set of quantum Langevin and Lyapunov equations and retrieve information from the resulting states, e.g., symplectic eigenvalues, logarithmic negativity, quantum fidelity between two states, ...