Computational enhancement of discrete gradient method

We propose new numerical approach to non-conservative dynamical systems. Our method being of low order, enhances qualitative performance of standard discrete gradient algorithm, thank to new concept of a reservoir. Paper is of explanatory character, focusing on concrete non-linear physical systems. Superiority of our new method with respect to standard discrete gradient method is observed.