Describir: Energy preserving/decaying time integration schemes for multibody systems dynamics