This example shows how the __float128 from gcc libquadmath can be used with odeint.
Examplary implementation of the method of Heun.