提高odeint的精确度

Question

我使用odeint在Python 3中解决了四个耦合的ODE，得到N1 [x]，Ntau [x]，Nmu [x]和Ne [x]。我在x轴上解决它们在0.1到8之间，最后我只对Ntau [8] + Nmu [x] + Ne [x]的值感兴趣。

不幸的是，其中一些功能的顺序为+/- 10000，而总和通常为1e-14。这意味着在x = 8处评估的总和由数值误差支配。

我的实施相当复杂，涉及许多不相关的细节，所以我不会在这里分享（除非有要求），而是要求提供一般性建议。