Text this: Adaptive integration algorithm for stiff ordinary differential equations