[WIP] Some error in the numerical integration of the problem.
run illustrative_example.py to see what is wrong