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