https://keroseneheaterreview.blogspot.com/2023/10/solving-ode-numerically-with-python.html