在求解 ODE 时，我想要 40 位或更多小数位的工作精度。谁能告诉我是否可以使用 gsl ODE 求解器（也许使用 libquadmath）获得更高精度的结果？从http://www.gnu.org/software/gsl/manual/gsl-ref.html#Driver的 gsl_odeiv2 文档中可以看出，它们只需要双打，但真的是这样吗？
如果没有，有人可以提出替代方案吗？'odeint' 之类的东西是否支持高于双精度？我必须从头开始编写自己的吗？
（作为参考，我的具体问题是：Writing Runge-Kutta ODE solver using gsl，现在由于 Daniel Fischer 修复了错误）