ТЭЗҮ e pansion
Martin Berz, in Advances in Imaging and Electron Physics, 1999. 5.4.1 Differentiating ODE Solvers. Intuitively, the most direct method for obtaining Taylor expansions for the flow of an ODE is to recognize that a numerical ODE solver describes a functional dependency between initial conditions and final conditions.