  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.
