As discussed here: https://github.com/easyscience/dynamics-lib/discussions/145
As discussed here: #145